レポート課題

No.5 での課題の他、

$$\frac{ d^2 y}{d x^2 } + \sin (y) =0$$ (2)

（大振幅の振子の問題;初期値は各自与えること） をオイラー法、修正オイラー法で扱ってみよ。刻み幅に対して誤差がどのように 収束するかも調べること。

単に数字を出力結果として出すだけでなく、gnuplot などでグラフの形にしてみ よ。調和振動子や、振り子の課題では、$x,y(x)$ のグラフだけでなく、 $y,dy/dx$ (位相空間)でも表示して、物理的意味の考察もしてみよ。振り子では、 大振幅で回転運動になることに注意。

余裕があれば、 4次の Runge-Kutta 法にも取り組んでみよ。

この他、各自物理的に興味ある常微分方程式も扱ってみよ。

レポート提出期限は 6月20日まで。