"""
万有引力の計算
F = G * (m1 * m2) / r^2
万有引力 = 重力定数G × (物体1の質量 × 物体2の質量) ÷ (物体間の距離)の2乗
"""
# グラフを描くためにmatplotlibをインポート
import matplotlib.pyplot as plt
# グラフを描くための関数
def draw_graph(x, y):
plt.plot(x, y, marker="o")
plt.title("Gravity force")
plt.xlabel("Distance(km)")
plt.ylabel("Gravitational force(N)")
plt.show()
# 万有引力の計算をするための関数
def calculate_F():
r = range(6400,64001,5000)
F = [] # 一旦、万有引力Fを空リストで作っておく
G = 6.672 * (10**-11) # 重力定数
m1 = 5.972 * (10**24) # 地球の重さ
# 以下で物体の重さを入力させる
# floatに変換しているため、文字を入力するとエラー終了する
m2 = float(input("Please enter mass:"))
# 先ほどrangeで作成した距離のリストrを一つずつ取り出して、forceとして万有引力を計算
for distance in r:
force = G*(m1*m2)/((distance*1000)**2)
F.append(force) # 作成しておいたリストFにひとつずつ追加(append)していく
#グラフを描画
draw_graph(r, F)
# メイン関数
if __name__ == '__main__':
calculate_F()