计算不出结果

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

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

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

01program ex_Ci
02implicit none
03    real Ls,Ci,Vs,Rl,Rv,Ht,hli,st,H0,LV
04    real,intrinsic::log
05    read(*,*) Ls,Vs,Rl,Rv,Ht,hli
06    Ci = exp(-4.531+1.6562*(H0)+5.5496*(H0)**2-6.4695*(H0)**3+&
07        (-0.474675+0.079*(H0)-1.39*(H0)**2+1.3212*(H0)**3)*log(LV)&
08        +(-0.07291+0.088307*(H0)-0.49123*(H0)**2+0.43196*(H0)**3)*&
09        (log(LV))**2)
10    print *,"标准系数Ci:",Ci
11end program
12function H0(Ht,hli)
13implicit none
14real Ht,hli
15real H0
16H0=Ht-hli
17end function
18 
19function LV(Ls,Vs,Rl,Rv)
20implicit none
21real Ls,Vs,Rl,Rv
22real LV
23LV = (Ls/Vs)*sqrt(Rl/Rv)
24end function

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

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

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

怎么改呢？？

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

我猜测应该是这样：

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

01program ex_Ci
02 implicit none
03     real Ls,Ci,Vs,Rl,Rv,Ht,hli,st,H0,LV
04     real :: h0_f , lv_f
05     real,intrinsic::log
06     read(*,*) Ls,Vs,Rl,Rv,Ht,hli
07     h0 = h0_f( ht , hli )
08     Lv = Lv_f( Ls,Vs,Rl,Rv )
09     Ci = exp(-4.531+1.6562*(H0)+5.5496*(H0)**2-6.4695*(H0)**3+&
10         (-0.474675+0.079*(H0)-1.39*(H0)**2+1.3212*(H0)**3)*log(LV)&
11         +(-0.07291+0.088307*(H0)-0.49123*(H0)**2+0.43196*(H0)**3)*&
12         (log(LV))**2)
13     print *,"标准系数Ci:",Ci
14 end program
15 function H0_f(Ht,hli)
16 implicit none
17 real Ht,hli
18 real H0_f
19 H0_f=Ht-hli
20 end function H0_f
21 
22 function LV_f(Ls,Vs,Rl,Rv)
23 implicit none
24 real Ls,Vs,Rl,Rv
25 real LV_f
26 LV_f = (Ls/Vs)*sqrt(Rl/Rv)
27 end function LV_f

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

可以了谢谢楼主！！！

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

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

规矩勋章

元老勋章

水王勋章