プログラムのエラー理由（fortran）を教えてください、お願いします

Question

前提・実現したいこと

CygwinとFORTRANとgnuplotで振り子の動き（tとθのグラフ）を観測したいのですが、エラーメッセージがds1～の行ででてきます。

発生している問題・エラーメッセージ

Unclassifiable statement at(1) ※1は表示されない

該当のソースコード

open(10,file='kadai.dat')
n=1000
pi=3.14159
l=1.0
g=9.807
t=0
tf=1
h=(tf-t)/n
s=pi/4
z=0

write(*,*)t,s,z
do i=1,n
t=t+h

ds1=h*g(t,s,z)
dz1=h*f(t,s,z)

ds2=h*g(t+h/2,s+1./2*ds1,z+1./2*dz1)
dz2=h*f(t+h/2,s+1./2*ds1,z+1./2*dz1)
ds3=h*g(t+h/2,s+1./2*ds2,z+1./2*dz2)
dz3=h*f(t+h/2,s+1./2*ds2,z+1./2*dz2)
ds4=h*g(t+h,s+ds3,z+dz3)
dz4=h*f(t+h,s+ds3,z+dz3)
s=s+(ds1+2*ds2+2*ds3+ds4)/6
z=z+(dz1+2*dz2+2*dz3+dz4)/6
write(*,*)t,s,z
end do
stop
end

function g(t,s,z)
g=z
return
end

function f(t,s,z)
f=-g*s/l
return
end

試したこと

うち間違いかなと思って何回か打ち直しましたが上手くいきませんでした。

補足情報（FW/ツールのバージョンなど）

初心者なもので、よくわかりません、すみません…

Answer 1

重力加速度gと関数gが被っていたためエラーが出たようです