時間、速度、角度を出力するつもりでやったのですがコンパイルした時の値が違います
おそらく運動方程式が間違っているのだと思いますが(t,s辺りが違い気がする)どこがどう違うか教えて下さい
real L
data g, L / 9.8, 1.0 /
write (＊,＊) ’ Δt (sec) ＝’
read (＊,＊) dt
＊initial
2 write (＊,＊) ’ θ0 (radian) ＝’
read (＊,＊) s0
s ＝s0
v ＝0
t ＝0
vL＝v
＊ repeat
1 a＝ −g＊sin(s)
v ＝a＊t
w＝v/L
s ＝s0+v＊t/L
t ＝t+dt
write(＊,＊)’t,v,s’
write (＊,’(3f8.4)’)t,v,s
if (vL.gt.0.0.and.v.lt.0.0) go to 2
vL＝v
go to 1
end