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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 大家帮我看一个fortran的程序，我总是计算不出正确的结果
53/1返回列表
查看: 1617  |  回复: 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

至尊木虫 (职业作家)
★ ★
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，表示程序就此中止了，一般用在出现错误无法继续的情况下……
赞一下(2人) 回复此楼
21楼2011-12-12 15:53:06
 已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 先中青基、再中博后面上的情况也挺多吧？ +6 WKKWKK 2024-06-29 6/300 2024-06-29 13:45 by 好人与坏人
[访问学者] 蹲! 希望今天能有大红奖状 +15 326lhpqk 2024-06-28 19/950 2024-06-29 11:47 by wenxixiaoer
[论文投稿] 关于论文第一通讯作者的问题？ +11 winsaint 2024-06-26 11/550 2024-06-29 11:25 by vmuch520
[考博] 30岁 想读非全博，求建议 +7 八千里路云雨泥 2024-06-27 19/950 2024-06-29 11:24 by 八千里路云雨泥
[基金申请] 我劝基金委要善良+一个大胆的设想 +13 我是大王 2024-06-28 14/700 2024-06-29 09:32 by Higraduate
[基金申请] 恳请大家帮我看看如何提高 +10 drwuwenhan 2024-06-28 13/650 2024-06-28 22:14 by 好人与坏人
[基金申请] 有当年结题的在研国基对申请影响大不大 1+4 sunjc 2024-06-28 13/650 2024-06-28 17:09 by lijian7338
[硕博家园] 哭到鼻子不透气，要窒息 +3 pvrw0224 2024-06-28 3/150 2024-06-28 17:01 by 仲夏夜的星星
[基金申请] 博后面上今天下午会公布吗？大家有无消息？ +11 地理学1995 2024-06-24 16/800 2024-06-28 14:54 by 地理学1995
[考博] 25申请博士 +3 lxcqwer 2024-06-27 3/150 2024-06-28 12:28 by highxixi
[基金申请] 要持续整治滥发“帽子”、“牌子”之风 +9 babu2015 2024-06-25 9/450 2024-06-28 09:11 by newfuzzy1
[论文投稿] 论文润色 +3 自强不息a?a 2024-06-26 3/150 2024-06-28 01:23 by Cherry_1215
[论文投稿] 为何全世界都用JCR分区，咱们要弄个中科院分区？意义何在？ +6 意得辑_editage 2024-06-25 7/350 2024-06-27 20:59 by keyaner23
[基金申请] 基金评审时，评审专家还回去看申请人代表作之外的文章吗？ +12 孤单12站 2024-06-26 18/900 2024-06-27 18:30 by 孤单12站
[基金申请] 今年什么时候会评啊 +9 lancet0903 2024-06-24 9/450 2024-06-27 13:22 by yrzhu1984
[有机交流] 做什么表征可以检测塑料中碳氢氧的含量 5+3 pzr的sci之路 2024-06-25 7/350 2024-06-26 19:32 by wrgeng
[基金申请] 青年基金E02口青基去年几个函评专家？ +6 他山攻玉之石 2024-06-25 9/450 2024-06-26 15:09 by 他山攻玉之石
[基金申请] 说博后基金7月出的真打电话了吗？ +12 antonysole 2024-06-24 14/700 2024-06-24 13:39 by sizhouyi
[有机交流] 生成亚胺的反应怎么能进行完全 +3 1369836 2024-06-23 3/150 2024-06-23 18:44 by hwqMSE
[论文投稿] OSA期刊审稿逾期 +3 Thomas_Squid 2024-06-22 3/150 2024-06-23 15:20 by wspglt
信息提示
关闭
请填处理意见
关闭