function dx=star(t,x) %建立函数G=1;m1=1;m2=1;m3=1
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);(-G*m2*(x(5)-x(1)))/((r12)^3)-(G*m3*(x(9)-x(1)))/(r13)^3; x(4);(-G*m2*(x(7)-x(3)))/((r12)^3)-(G*m3*(x(11)-x(3)))/(r13)^3; x(6);(-G*m1*(x(1)-x(5)))/((r12)^3)-(G*m3*(x(9)-x(5)))/(r23)^3; x(8);(-G*m1*(x(3)-x(7)))/((r12)^3)-(G*m3*(x(11)-x(7)))/(r23)^3; x(10);(-G*m1*(x(1)-x(9)))/((r13)^3)-(G*m2*(x(5)-x(9)))/(r23)^3; x(12);(-G*m1*(x(3)-x(11)))/((r13)^3)-(G*m2*(x(7)-x(11)))/(r23)^3;]end
figure('name','三体运动');
x0=[0.3;-0.466;0.243;-0.433;0.970;-0.466;-0.243;-0.433;0;0.932;0;0.866] %调用函数
[t,x]=ode45(@star,[-1000,1000],x0)
plot(x(:,1),x(:,3),x(:,5),x(:,7),x(:,9),x(:,11))
你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答
本次提问扣除的有问必答次数，已经为您补发到账户，我们后续会持续优化，扩大我们的服务范围，为您带来更好地服务。
点赞