24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1977  |  回复: 20
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

princeado

金虫 (正式写手)

[求助] 大家帮我看一个fortran的程序,我总是计算不出正确的结果

最近刚学fortran,开始编小小程序试手,下面是我弄的一个程序,没有错误,但总是运行不出正确的结果,程序是想通过do循环对数组a进行赋值,用的是语句a(i)=jiecheng(i),jiecheng是我自定义的函数,求每个数对应的阶乘。所有程序如下,哪位高手帮我看看问题在哪里。
        program test01
implicit none
integer :: i
real :: a(5)
real,external :: jiecheng
     do i=1,5
       a(i)=jiecheng(i)   
           end do
             write(*,*) a
           stop
           end
                   function jiecheng(x)
                   implicit none
                   integer x
            real jiecheng
           integer :: ii=1
           integer :: m=1
                   do ii=1,x
                       m=m*ii
                       end do
                           jiecheng=m
                                                  return
                           stop
                           end
回复此楼

» 猜你喜欢

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

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

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖


princeado(金币+20): 谢了,问题解决了,我再好好研究下原因 2011-11-08 16:10:31
余泽成(金币+1): 谢谢参与应助! 2011-11-08 18:10:24
内容已删除
2楼2011-11-08 15:50:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)


余泽成(金币+1): 谢谢参与应助! 2011-11-08 18:10:46
fortran 对声明阶段赋值的变量自动赋予 SAVE 属性,这样,在这个子程的连续调过程,在下一次执行时,其初始值为一次调用结束时的值……

另外,子程序的 return 之后,不要用 stop,stop 只在应该结束程序时才使用……
5楼2011-11-08 16:18:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

★ ★
余泽成(金币+2): 谢谢参与应助! 2011-11-08 18:10:59
余泽成: 专家考核存档! 2011-11-08 18:11:11
引用回帖:
7楼: Originally posted by princeado at 2011-11-08 16:37:35:
还有就是,integer :: m
                   m=1
与 integer :: m=1 两种声明变量的方式的区别是什么啊?

区别在于
integer :: m
m=1
中的 m是一个局变变量,不具有 SAVE 属性,每次调用该子程时,m 都会被重新赋值。

而 integer :: m=1 中的 m 虽然是一个局部变量,但具有 SAVE 属性,这里的 1 只是主程序初始化时的取值,在随后的执行中,它的值会是上一次该子程序执行结束时的值。

至于什么是 SAVE 属性,去查一下标准吧,我还真讲不清楚。大概的意思是,该变量的最值取值会被保留下来,直到主程序结束。
9楼2011-11-08 16:53:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)


余泽成(金币+1): 鼓励讨论! 2011-11-13 21:03:21
引用回帖:
15楼: Originally posted by 长歌-废墟 at 2011-11-10 14:07:03:
建议使用call子程序的写法,这样基本不会出错。

他这个问题跟用 function 还是 subroutine 没有关系。
17楼2011-11-10 18:47:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)


余泽成(金币+1): 鼓励交流! 2011-11-13 21:03:41
引用回帖:
13楼: Originally posted by princeado at 2011-11-09 07:37:21:
晚上回去翻书也看到了啊,我这fortran初学者就遇到了这样的问题,昨天下午还给彭国伦发邮件请教去了,呵呵

我很好奇,他给你啥样的回复?

这个 save 属性是标准中要求的,所有的编译器都是一样的行为……
18楼2011-11-10 18:48:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

★ ★
jjdg(金币+1): 感谢参与 2011-12-12 22:29:15
jjdg(金币+1): 感谢应助 2011-12-13 00:50:03
引用回帖:
20楼: Originally posted by 1303haojian at 2011-12-12 14:35:40:
函数里面也要用stop吗   不是return   end就可以了吗

可以用 stop,表示程序就此中止了,一般用在出现错误无法继续的情况下……
21楼2011-12-12 15:53:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 princeado 的主题更新
信息提示
请填处理意见