すいません。
自分でなんとか書いてみたんですが、合ってるかどうか見てもらえますか？
本当ならこれを学校のFortranでやってみて、確かめて、gnuplotでグラフ。
とまでやりたかったんですが、学校のが使えませんでした。（家も使えません。）
どなたか数値計算していただけない・・ですよね。。
もし計算結果がでたらエクセルでグラフを書けばなんとかできるかなと思っているんですが、
↓こんなんです