这个程序运行不起来，是matlab用微分方程解三体运动

Question

function dr=star(t,x) %建立函数

G=6.67*10^(-11);

m1=0.3;

m2=0.03;

m3=0.03

r12=sqrt((x(5)^2-x(1)^2)^2+(x(7)^2-x(3)^2)^2);

r23=sqrt((x(9)^2-x(5)^2)^2+(x(11)^2-x(7)^2)^2);

r13=sqrt((x(9)^2-x(1)^2)^2+(x(11)^2-x(3)^2)^2);

dx=[x(2)+x(4);(-G*m2)/((r12)^2)-(G*m3)/(r13)^2;

    x(6)+x(8);(-G*m3)/((r23)^2)-(G*m1)/(r12)^2;

    x(10)+x(12);(-G*m1)/((r13)^2)-(G*m2)/(r23)^2]

end

x0=[2;0;-2;0.2;0;-0.2;2;0;-2;-0.2;0;0.2] %调用函数

[t,y]=ode45(@star,[0,1000],x0)

plot(y(:,1),y(:,3))

figure('name','三体运动');

axis equal

Answer 1

你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答

本次提问扣除的有问必答次数，已经为您补发到账户，我们后续会持续优化，扩大我们的服务范围，为您带来更好地服务。