計算機数学レポートの答え

計算機数学のソース（ルンゲクッタ法？による 微分方程式の解法）
手計算でといたところ、V(t)=-0.75*t*exp{(-r/2l)t}
のようです。
（注）タグでのタブの打ち方がわかりません。適当に調整してください。
ただしx0=0,y0=0(Vにあたる）、z0=0(dV/dtにあたる）,dx=0.000001、xend=0.0001

以下本文

real dx,xend,x,x0
real y0,y,z0,z
read (5,*)x0,y0,z0,dx,xend


x=x0
y=y0
z=z0

do 10 x=x0,xend,dx
write(6,600)x,y,0.000001*z
600 format(' ',f10.8)
y=y0+dx*z0
z=z0-dx*(10**5)*(2*z0+y0*10**5-10**5)
y0=y
z0=z
10 continue

stop
end