3楼: Originally posted by chembetsey at 2012-07-27 23:03:38
不知道你想做什么，你的第二个方程是一个微分方程，如果除了L之外都是常数的话，你可把它积分出来，写出Fortran的格式是：
log(L/L0) = (Alpha*dT-Sigma/(3*K))*t  或
L = L0 * exp((Alpha*dT-Sigma/(3*K))*t)
当 ...

5楼: Originally posted by chembetsey at 2012-07-28 11:42:26
写成子函数就可以了：
Real*8 Function Delta(Alpha, dT, Sigma, K, t)
real*8 Alpha, dT, Sigma, K, t
Delta = exp((Alpha*dT-Sigma/(3.D0*K))*t) -1.D0
End