各位，帮忙看一下这个微分方程哪里错了

[复制链接]

白影2

2 主题	4 听众	63 积分

收听TA
发消息

群组: Matlab讨论组

电梯直达

1^#

发表于 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)]