[Fortran] 纯文本查看 复制代码 SUBROUTINE FUN(X,Y,A,B)
DOUBLE PRECISION X,Y,A,B
Double COMPLEX ss,ff,lamda
double COMPLEX pt,fenzi,fenmu
DOUBLE PRECISION a0,r
a0=10.0
r=5.0
ss=CMPLX(X,Y) !+(r-a0)/cp
lamda=sqrt(ss**2*(ss+2)/(1+ss))
fenzi=cosh(lamda*(r-a0))
fenmu=cosh(lamda*a0)
c pt=-0.04*(1/ss-exp(-6.25*ss)/ss)
pt=(54-58*exp(-12.5*ss)+4*exp(-125*ss))/(4500*ss**2)
WRITE(10,*)pt !NEW CODE
ff=pt*fenzi/fenmu
a= real(ff)
b= imag(ff)
! RETURN
END
不知道为什么,用Simply Fortran调试的时候,加上write就能正常运行了,就很奇怪( |