Rubyのルンゲクッタ法がうまくいきません

質問

Rubyでルンゲクッタ法でy=exp(x**2)と比較しようという問題で、下のようなプログラムを組むとyの値がかなり大きくなってしまいます。どこが間違っているのでしょうか。教えてください。

#! ruby -Ks

#ルンゲックッタ
def df(a,b)
z=2*a*b
return z
end

n=10
x=Array.new(n)
y=Array.new(n)

x[0]=0.0
y[0]=1.0

h=1.0/n #任意で変更

fw=File.open("output2.txt","w")

for i in 0..n
k1=df(x[i],y[i])*h
k2=df(x[i]+h/2,y[i]+k1/2)*h
k3=df(x[i]+h/2,y[i]+k2/2)*h
k4=df(x[i]+h,y[i]+k3)
k=(k1+2*k2+2*k3+k4)/6

x[i+1]=x[i]+h
y[i+1]=y[i]+k
fw.print x[i]," ",y[i]," ","\n"
end

fw.close()

回答No.1

Tacosan · ベストアンサー

k4

お礼

いち早い指摘ありがとうございました。無事できました

回答No.2

MillenniuM

すでに答えは出てますが、私なりに考えたものを
http://pastebin.com/kgAPLFdd

お礼

ご呈示いただいたプログラムはより分かり良いものでした。ありがとうございます。