fortran ödevime bir yardım!

[sicilium]

merhaba arkadaşlar,

bana fortranda alt fonksiyon yapma syntaxını bir örnek de vererek anlatır mısınız?

daha doğrusu bir program yazdım ama parce error veriyor. Fonksiyonun return yaptığı yerde bunu söylüyor. Bu yüzden fonksiyon syntaxına ihtiyacım var.

Yardımlarınız için çok teşekkürler:)

program homework

real a,b,t1,d,tau,t2ilk,t2,cpilk,p1,p2
a=3.376
b=0.557*((10)**(-3))
d=-0.031*10**(-5)

t1=300
t2=0

tau=1.1
cpilk=a+(b*t1+(d/(tau**2*t1**2))*((tau+1)/2))*((tau-1)/log(tau))
t2ilk=exp((1/cpilk)*log(p2/p1))
t2=t2ilk

1 tau=t2/t1
t2ilk=exp((1/cp(tau))*log(p2/p1))
if(((t2ilk/t2)/t2ilk).ge.10**(-3)) goto 1

print*,"outlet temp=",t2ilk
stop
end program

function cp(tauu)
real tauu, cpig
cpig=a+(b*t1+(d/(tauu**2*t1**2))*((tauu+1)/2))*((tauu-1)/log(tauu))
cp=cpig
return
end function