24小时热门版块排行榜    

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

LZQ19871

木虫 (小有名气)

[交流] 【求助】fortran程序未出错,但运行不出结果,中途中断!

各位请教下
   我用“徐士良”,“积分一步的变步长龙格-库塔法”解常微分方程组,正常情况下解得t=0,p=1......t=399.99,p=0.50774,t=400,p=0.507614. 而当我改变常微方程组系数时,解到t=3.01,p=0.999047时,程序就好像卡在那了,运行很久也没反应,也没有显示“Press any key to continue”
     同样的参数,我用MATLAB“ ode45”求解可以得到结果。
请问下有可能是哪方面的原因?

   附件中位FORTRAN程序,能运行出结果,可以直接打开EXCEL,查看。
但把main.f90中的 rd=0.4*w0 改成 4*w0 运行一会就不运行了,打开EXCEL发现停在某个数值上了。一直等都没有结果,也没有任何出错信息。
    附件:http://good.gd/690578.htm

[ Last edited by LZQ19871 on 2010-9-3 at 08:39 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

LZQ19871(金币+5): 2010-09-15 09:42:37
给个能用 gfortran 编译的版本吗?呵呵……
4楼2010-09-12 15:15:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

hakekill

木虫 (小有名气)


余泽成(金币+1):谢谢参与应助! 2010-09-10 16:16:57
LZQ19871(金币+5): 2010-09-15 09:42:26
程序是不是陷入死循环了。你可以单步调试看看。
2楼2010-09-10 11:39:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

LZQ19871

木虫 (小有名气)

没有陷入死循环

我可以确认没有陷入死循环,
    因为这个程序开始是可以运行,能出结果的,但我只把常微分方程组中某几项系数改变,循环没有变化。
3楼2010-09-12 14:54:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ifmc1234

木虫 (著名写手)

引用回帖:
Originally posted by snoopyzhao at 2010-09-12 15:15:43:
给个能用 gfortran 编译的版本吗?呵呵……

你要在什么环境下运行呢?有代码可以自己编译的吧
5楼2010-09-12 19:23:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 8/400 2026-08-25 05:12 by BZKMTicpDhFj
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见