计算不出结果

中森天下 · 发表于 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

fcode · 发表于 2015-9-1 17:13:37

H0 和 LV 作为函数，你需要给予传递参数。

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

怎么改呢？？

fcode · 发表于 2015-9-1 17:40:45

我猜测应该是这样：

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

program ex_Ci
 implicit none
     real Ls,Ci,Vs,Rl,Rv,Ht,hli,st,H0,LV
     real :: h0_f , lv_f
     real,intrinsic::log
     read(*,*) Ls,Vs,Rl,Rv,Ht,hli
     h0 = h0_f( ht , hli )
     Lv = Lv_f( Ls,Vs,Rl,Rv )
     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_f(Ht,hli) 
 implicit none
 real Ht,hli
 real H0_f
 H0_f=Ht-hli
 end function H0_f

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

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

可以了谢谢楼主！！！

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

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

规矩勋章

元老勋章

水王勋章