应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 大家帮我看一个fortran的程序，我总是计算不出正确的结果
51/1返回列表
查看: 1614  |  回复: 20
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2011-11-08 15:26:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

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

另外，子程序的 return 之后，不要用 stop，stop 只在应该结束程序时才使用……
赞一下(1人) 回复此楼
5楼2011-11-08 16:18:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 21 个回答

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★
princeado(金币+20): 谢了，问题解决了，我再好好研究下原因 2011-11-08 16:10:31
余泽成(金币+1): 谢谢参与应助！ 2011-11-08 18:10:24
把 jiecheng 中的下列语句
[cdoe]
     integer :: ii=1
     integer :: m=1
[/code]
修改成下面就可以了：
CODE:
      integer :: ii
      integer :: m
      m=1

至于为什么，可以参考这里的讨论
http://www.newsmth.net/nForum/article/NumComp/102903
赞一下(1人) 回复此楼
2楼2011-11-08 15:50:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hhucy

实习版主 (著名写手)

DOE锅炉工
引用回帖:
2楼: Originally posted by snoopyzhao at 2011-11-08 15:50:49:
把 jiecheng 中的下列语句
[cdoe]
     integer :: ii=1
     integer :: m=1
[/code]
修改成下面就可以了：
CODE:
      integer :: ii
      integer :: m
      m=1

至于为什么，可以参考 ...

好可怕的变量
回复此楼
人生那么多不确定，你怕什么
3楼2011-11-08 16:06:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

princeado

金虫 (正式写手)
引用回帖:
2楼: Originally posted by snoopyzhao at 2011-11-08 15:50:49:
把 jiecheng 中的下列语句
[cdoe]
     integer :: ii=1
     integer :: m=1
[/code]
修改成下面就可以了：
CODE:
      integer :: ii
      integer :: m
      m=1

至于为什么，可以参考 ...

原因不是很清楚，您能再帮我解释一下吗？
赞一下 回复此楼
4楼2011-11-08 16:18:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 21 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 博后网上办公系统维护 +13 取名字烦人 2024-06-25 21/1050 2024-06-26 19:41 by msjy
[基金申请] 今年什么时候会评啊 +8 lancet0903 2024-06-24 8/400 2024-06-26 19:11 by binfriedman
[硕博家园] 怎么带研究生？ +24 豆豆小小2 2024-06-20 42/2100 2024-06-26 16:13 by wangyu1986w
[基金申请] 博后面上今天下午会公布吗？大家有无消息？ +9 地理学1995 2024-06-24 13/650 2024-06-26 16:02 by kyukitu
[基金申请] 博后面上和特助今天出吗？ +41 逗您玩 2024-06-21 78/3900 2024-06-26 16:00 by HAPPY_0225
[基金申请] 要持续整治滥发“帽子”、“牌子”之风 +6 babu2015 2024-06-25 6/300 2024-06-26 14:52 by felicity6056
[硕博家园] 【45岁以上博士】柔性人才引进项目，有补贴 +9 Dreamsummit 2024-06-24 18/900 2024-06-26 12:28 by hujm159
[硕博家园] 考博英语 +5 198新用户 2024-06-25 10/500 2024-06-26 08:09 by 地狱伞兵
[硕博家园] 申博雅思求助 +6 刘嘉研究究 2024-06-24 7/350 2024-06-26 07:08 by 小李醒yy
[硕博家园] 夏至，要不要硕博联谊 +4 我是王小帅 2024-06-21 6/300 2024-06-25 14:25 by 庭前花开月儿明
[教师之家] 有没有今年的影响因子？ +6 jurkat.1640 2024-06-22 10/500 2024-06-24 22:48 by 潇湘之迷
[基金申请] 国自然资助比率是不是要下降了？？ +8 今晚推荐22 2024-06-21 12/600 2024-06-24 11:15 by Pickfoot
[有机交流] 三颈瓶加热 5+4 wzjwx888 2024-06-20 4/200 2024-06-24 10:11 by w87437985
[论文投稿] 论文提交二审还有三天就三个月了，连续问了编辑部几次 10+3 大王叫我来寻山� 2024-06-22 9/450 2024-06-24 08:50 by 大王叫我来寻山�
[基金申请] 国自然青年基金，1A4B能上会吗？青年和面上的上会标准是一样的吗？ +19 今晚推荐22 2024-06-20 32/1600 2024-06-23 23:17 by andywei1028
[论文投稿] OSA期刊审稿逾期 +3 Thomas_Squid 2024-06-22 3/150 2024-06-23 15:20 by wspglt
[基金申请] F03青年基金函评结果 +5 暨阳一只柴 2024-06-19 6/300 2024-06-23 14:30 by adsqsj
[有机交流] 求助 45+7 脂质纳米粒 2024-06-20 9/450 2024-06-23 07:52 by buhui7829
[基金申请] 教育部基金 +5 m1393 2024-06-21 5/250 2024-06-21 21:13 by odes
[基金申请] 听大佬说今年信息口本子数量大幅增加？ +8 wutzxt 2024-06-21 9/450 2024-06-21 19:58 by wutzxt
信息提示
关闭
请填处理意见
关闭