| 查看: 1798 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | ||
liuhaizhen木虫 (初入文坛)
|
[求助]
求mathematica编程
|
|
|
dy(1)/dt=(a+b+c)*y(1)+g*y(2)+b*y(3) dy(2)/dt=-(a+b)*y(3) dy(3)/dt=(a-b-c)*y(3)-g*y(1)+2*g*y(2)*y(1) 其中a,b,c,g为常量 |
» 猜你喜欢
Cas 72-43-5需要30g,定制合成,能接单的留言
已经有8人回复
求助:我三月中下旬出站,青基依托单位怎么办?
已经有6人回复
北京211副教授,35岁,想重新出发,去国外做博后,怎么样?
已经有8人回复
磺酰氟产物,毕不了业了!
已经有5人回复
论文终于录用啦!满足毕业条件了
已经有25人回复
2026年机械制造与材料应用国际会议 (ICMMMA 2026)
已经有3人回复
自荐读博
已经有3人回复
不自信的我
已经有5人回复
投稿Elsevier的杂志(返修),总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐,对您同样有帮助:
【软件】Mathmatica 软件 7.0版
已经有68人回复
mathematica 编程求助
已经有9人回复
如果编程语言是汽车……
已经有16人回复
这个复杂式子如何计算出求导结果
已经有5人回复

6楼2011-11-17 22:02:13
mshwangg
至尊木虫 (正式写手)
- 程序强帖: 5
- 应助: 206 (大学生)
- 金币: 10702.8
- 散金: 100
- 红花: 19
- 帖子: 597
- 在线: 195.4小时
- 虫号: 576702
- 注册: 2008-06-21
- 专业: 物理学I
2楼2011-11-17 00:30:03
【答案】应助回帖
★ ★
余泽成(金币+2): 谢谢参与应助! 2011-11-19 14:28:38
liuhaizhen(金币+10): 感谢你的应助 2011-11-28 10:00:33
liuhaizhen(金币+5): 优秀 2011-12-07 10:16:41
余泽成(金币+2): 谢谢参与应助! 2011-11-19 14:28:38
liuhaizhen(金币+10): 感谢你的应助 2011-11-28 10:00:33
liuhaizhen(金币+5): 优秀 2011-12-07 10:16:41
|
下面是程序代码,语法肯定没错,只要你再给定初始条件,就可以得到结果。否则解出来的是通解。 eq1 = y1'[t] == (a + b + c)*y1[t] + g*y2[t] + b*y3[t] eq2 = y2'[t] == -(a + b)*y3[t] eq3 = y3'[t] == (a - b - c)*y3[t] - g*y1[t] + 2*g*y2[t]*y1[t] DSolve[{eq1, eq2, eq3}, {y1[t], y2[t], y3[t]}, t] |

3楼2011-11-17 09:46:25

4楼2011-11-17 09:48:37









回复此楼