数学建模社区-数学中国

标题: 各位，帮忙看一下这个微分方程哪里错了 [打印本页]

作者: 白影2 时间: 2012-10-23 15:03
标题: 各位，帮忙看一下这个微分方程哪里错了
clear,clc
x0=[-6045;-3.457;-3490;6.618;2500;2.533];
options=odeset('reltol',1);
tic
[t,y]=ode45(@rigid,[-100000,100000],x0,options)
toc
plot3(y(:,1),y(:,2),y(:,3))

function dx=rigid(t,x)
u=398600;
r1=sqrt(x(1)^2+x(3)^2+x(5)^2);
dx=[x(2)
(u/r1^3)*x(1)
x(4)
(u/r1^3)*x(3)
x(6)
(u/r1^3)*x(5)]

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)