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