>>435なんですが
１：０〜２０の区間を適当に分割
２：初期値入力
３：２階常微分方程式を一階常微分に直す
　　dy/dx=z dz/dx=f(x,y,z)=-sin(y)としました。
４：ルンゲクッタ法
　　y=y+(k1+2*k2+2*k3+k4)/6
k1=h*f(x,y)
k2=h*f(x+h/2,y+k1/2)
k3=h*f(x+h/2,y+k2/2)
k4=h*f(x+h,y+k3)
　　これの()内が３つあるバージョン。分割した回数繰り返す。
５：ファイルにデータ書き込み出力

こういう流れで作りたいんですがどこが間違ってるのか出力すると初期値を
繰り返したデータになってしまいます。２〜５あたりだけでもいいので
わかる方いたら教えてもらえませんか？
イメージとしては↓のサイトのようにやったんですがうまくいきませんでした
http://www.akita-nct.ac.jp/~yamamoto/lecture/2005/5E/differential_equations/text/html/node3.html