| 查看: 2899 | 回复: 2 | |||
[求助]
mathematica 使用NDSolve时出现的问题已有1人参与
|
| 在mathematica 9.0中，使用NDSolve 时，只是将边界条件中的0改为0.03后，为什么会出现图中的问题呢？ 在 x == 0.00116` 处，步长实际上为零；可能存在奇点或者刚性系统。怎么解决？？ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : mathematica.png
2014-04-12 17:22:45, 34.12 K
» 猜你喜欢
导师想让我从独立一作变成了共一第一
已经有9人回复
-
博士读完未来一定会好吗
已经有23人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
微分方程的求解
已经有9人回复- 请教mathematic问题？关于解微分方程组的和画图的一些问题。
已经有19人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kexiaolong: 金币+30 2014-12-11 09:22:01
kexiaolong: 金币+30 2014-12-11 09:22:01
|
0.03的问题我没仔细研究，但是从你的提问来看，你是因为0处失败了才转向了0.03吧？那么，对0处的失败，可以添加选项MaxSteps->Infinity来解决。不过，你确定你这方程没问题吗？： eqn = With[{i = 0}, Rationalize[{y''[x] - y[x]/0.00008585917334 + 1/30.02 Exp[38 + 32.633/(2.586 + 138 x + 25.655 x^0.732)] == 0, y' == 100/30.02 y, y[11.6 10^-4] == 0}, 0]]; sol = NDSolve[eqn, y, {x, 0, 400}, MaxSteps -> Infinity, WorkingPrecision -> 16]; y[100] /. sol (* {-9.6701192450*10^4701} *) 调了下精度发现冒似没什么改善……直觉上你最好再分析下你的方程。 |
2楼2014-07-03 18:51:18
3楼2014-07-03 18:52:49