import math g = 9.8 # 定数を宣言する print ("Velocity [m/s]?") V = float (input()) # 初速度の大きさを入力する print ("Angle [deg.]?") angle = float (input()) # 仰角を入力する theta = 3.14159 / 180.0 * angle x = 0.0 # 初期値を計算する z = 0.0 u = V * math.cos (theta) w = V * math.sin (theta) f = open ("output.txt", "w") # ファイルを開く f.write ("%9.3f %9.3f\n" % (x, z)) # 結果をファイルに出力する dxdt = u # 時間微分を計算する dzdt = w dudt = 0.0 dwdt = - g x = x + 0.01 * dxdt # 次の時刻の値を計算する z = z + 0.01 * dzdt u = u + 0.01 * dudt w = w + 0.01 * dwdt f.write ("%9.3f %9.3f\n" % (x, z)) # 結果をファイルに出力する f.close () # ファイルを閉じる