|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本帖最后由 maomao_hhu 于 2012-11-11 21:09 编辑
因为是非线性微分方程组的,所以老师说用deplot算一个six equations each with second order的微分方程组,
每个方程类似于:
de1 := diff(x1(t), t, t)+(1/2)*(sqrt((x3(t)-x1(t))^2+(y3(t)-y1(t))^2)-2)*(-2*x3(t)+2*x1(t))/sqrt((x3(t)-x1(t))^2+(y3(t)-y1(t))^2)+(1/2)*(sqrt((x2(t)-x1(t))^2+(y2(t)-y1(t))^2)-2)*(-2*x2(t)+2*x1(t))/sqrt((x2(t)-x1(t))^2+(y2(t)-y1(t))^2) = diff(x2(t), t)-2*(diff(x1(t), t))+diff(x3(t), t)
组成系统:
sys := {de1, de2, de3, de4, de5, de6}
最后输入命令:
DEplot([sys], [x1(t), x2(t), x3(t), y1(t), y2(t), y3(t)], t = 0 .. 100, [[x1(0) = 1, x2(0) = 3, x3(0) = 2, y1(0) = 1, y2(0) = 1, y3(0) = 2.73, (D(x1))(0) = 1, (D(x2))(0) = 0, (D(x3))(0) = 0, (D(y1))(0) = 0, (D(y2))(0) = 0, (D(y3))(0) = 0]])
结果出现:
Error, (in DEtools/DEplot) system must have same number of dependent variables as DE's.
我觉得变量没有问题啊,求各位高手指教
|
|