注册 登录
Simwe仿真论坛(forum.simwe.com),CAE/CAD/CAM/,FEA/FEM/有限元分析论坛---(手机验证注册) 返回首页

TBE_Legend的个人空间 http://home.simwe.com/?355196 [收藏] [复制] [分享] [RSS]

日志

一道DAE的maple解答和mmtc解答

已有 1052 次阅读2010-4-15 10:37 |个人分类:mmtc_maple|

Clear["Global`*"]
SetOptions[Plot, PlotRange -> Full, Frame -> True];
sys = {a'[t] == -(a[t]/10) + 10000 b[t] c[t],
   b'[t] == a[t]/10 - 10000 b[t] c[t] - 10000000 b[t]^2,
   a[t] + b[t] + c[t] == 1, a[0] == 1, b[0] == 0};
sol = First[NDSolve[sys, {a, b, c}, {t, 0, 40000}]];
Plot[a[t] /. sol, {t, 0, 40000}]
Plot[Evaluate[{a[t], b[t], c[t]} /. sol], {t, 0, 40000}]
 
 
 
 
 
 
 

路过

鸡蛋

鲜花

握手

雷人

评论 (0 个评论)

facelist doodle 涂鸦板

您需要登录后才可以评论 登录 | 注册

小黑屋|手机版|Archiver|联系我们|中国仿真互动网 ( 沪ICP备07510919号 ) 云安全

GMT+8, 2019-10-24 11:32 , Processed in 0.064702 second(s), 16 queries , Gzip On.

Powered by Discuz! X3.4 Licensed

© 2001-2017 Comsenz Inc.

返回顶部