计算不出结果

中森天下 · 发表于 2015-9-1 16:46:43

系统版本 win8.1 64位，编译环境microsoft visual studio2010 Intel.Visual.Fortran.Composer.XE.2011
是不是因为内部函数log需要另外定义一下？？求指点

[Fortran] 纯文本查看 复制代码

program ex_Ci
implicit none
    real Ls,Ci,Vs,Rl,Rv,Ht,hli,st,H0,LV
    real,intrinsic::log
    read(*,*) Ls,Vs,Rl,Rv,Ht,hli
    Ci = exp(-4.531+1.6562*(H0)+5.5496*(H0)**2-6.4695*(H0)**3+&
        (-0.474675+0.079*(H0)-1.39*(H0)**2+1.3212*(H0)**3)*log(LV)&
        +(-0.07291+0.088307*(H0)-0.49123*(H0)**2+0.43196*(H0)**3)*&
        (log(LV))**2)
    print *,"标准系数Ci:",Ci
end program
function H0(Ht,hli) 
implicit none
real Ht,hli
real H0
H0=Ht-hli
end function

function LV(Ls,Vs,Rl,Rv) 
implicit none
real Ls,Vs,Rl,Rv
real LV
LV = (Ls/Vs)*sqrt(Rl/Rv)
end function

中森天下 · 发表于 2015-9-1 17:35:47

怎么改呢？？

中森天下 · 发表于 2015-9-1 17:52:54

可以了谢谢楼主！！！

[数值问题] 计算不出结果