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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】fortran程序未出错，但运行不出结果，中途中断！
71/1返回列表
查看: 1778  |  回复: 6
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2010-09-02 21:05:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakekill

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

LZQ19871

木虫 (小有名气)

没有陷入死循环

我可以确认没有陷入死循环，
    因为这个程序开始是可以运行，能出结果的，但我只把常微分方程组中某几项系数改变，循环没有变化。
赞一下 回复此楼
3楼2010-09-12 14:54:36
 已阅   回复此楼   关注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的回帖

ifmc1234

木虫 (著名写手)
引用回帖:
Originally posted by snoopyzhao at 2010-09-12 15:15:43:
给个能用 gfortran 编译的版本吗？呵呵……

你要在什么环境下运行呢？有代码可以自己编译的吧
赞一下 回复此楼
5楼2010-09-12 19:23:03
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)
★
余泽成(金币+1):辛苦了！ 2010-09-13 11:03:49
LZQ19871(金币+10): 2010-11-12 09:42:57
就是你的代码用 gfortran 编译时出错，呵呵……
赞一下(1人) 回复此楼
6楼2010-09-12 21:02:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

LZQ19871

木虫 (小有名气)

我用的是 CVF IVF 我也分不清楚

我用的是 CVF  IVF 我也分不清楚
  不好意思，我没有你说的编译器。
赞一下 回复此楼
7楼2010-09-13 09:18:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 LZQ19871 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭