[交流] 问一个关于函数调用的小问题

请教个问题，看看我这个函数，fortran77写的：
   FUNCTION HEAVI(PRIS)
         REAL*8 PRIS
   IF(PRIS.GT.0.0) THEN
      HEAVI=1.0
   ENDIF
   RETURN
C
      IF(PRIS.LT.0.0) THEN
      HEAVI=0.0
      ENDIF
   RETURN
      END

调用的时候出现下列情况：
HEAVI(-3)=1.00000000000000，
HEAVI(-3.0D0)=0.000000000000000E+000，
   HEAVI(-3.0)=1.00000000000000
HEAVI(3)=1.00000000000000，
HEAVI(3.0D0)=1.000000000000000E+000，
HEAVI(3.0)=1.00000000000000
不知道什么原因？

改成这样怎么就好了呢？

   Program ex005
   real HEAVI
   write(*,*) HEAVI(3.0d0)
   end

   FUNCTION HEAVI(PRIS)
   REAL PRIS
   Real HEAVI
   IF(PRIS.GT.0.0) THEN
      HEAVI=1.0
   else
      HEAVI=0.0
   ENDIF
   RETURN
   END
语言Fortran77，gfortran环境。
求解释。

回复此楼

» 猜你喜欢

拟解决的关键科学问题还要不要写已经有9人回复
最失望的一年已经有17人回复
为什么nbs上溴没有产物点出现呢已经有6人回复
求推荐博导已经有4人回复
存款400万可以在学校里躺平吗已经有34人回复
求助一下有机合成大神已经有4人回复
求推荐英文EI期刊已经有5人回复
26申博已经有3人回复
基金委咋了？2026年的指南还没有出来？已经有10人回复
疑惑？已经有5人回复

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-11-11 21:53:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

maomao1210

金虫 (正式写手)

程序强帖: 5
应助: 2 (幼儿园)
金币: 1432.8
帖子: 991
在线: 442.1小时
虫号: 253215

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
jjdg(金币+1): 感谢参与 2011-11-12 00:41:19
余泽成: 专家考核存档 2011-11-13 21:35:24

引用回帖:

1楼: Originally posted by guojunhang at 2011-11-11 21:53:59:
请教个问题，看看我这个函数，fortran77写的：
   FUNCTION HEAVI(PRIS)
         REAL*8 PRIS
   IF(PRIS.GT.0.0) THEN
      HEAVI=1.0
   ENDIF
   RETURN
C
      IF(PRIS. ...

貌似变量定义也有问题，第一个函数中用的是real*8,主程序给的有时候是real*8 ,有时候是real*4

[ Last edited by maomao1210 on 2011-11-11 at 23:00 ]