function dr=star(t,r) %建立函数
g=10;
m1=0.3;
m2=0.03;
m3=0.03
r(1)=sqrt(x(1)^2+x(2)^2);
r(2)=sqrt(x(5)^2+x(6)^2);
r(3)=sqrt(x(9)^2+x(10)^2);
dr=[r(2);(-g*m2*(r(1)-r(3))/(r(1)-r(3))^3)-(g*m3*(r(5)-r(1))/(r(5)-r(1))^3);
r(4);(-g*m3*(r(3)-r(5))/(r(3)-r(5))^3)-(g*m1*(r(1)-r(3))/(r(1)-r(3))^3);
r(6);(-g*m1*(r(5)-r(1))/(r(5)-r(1))^3)-(g*m2*(r(3)-r(5))/(r(3)-r(5))^3)]
end
x0=[2;0;-2;0.2;0;-0.2;2;0;-2;-0.2;0;0.2] %调用函数
[t,y]=ode45(@star,[0,100],r0)
plot(y(:,1),y(:,3))
题目要求:
求解三体问题。设质量为m1=0.3,m2=m3=0.03.(a)画出初始条件(x1(0),y1(0))=(2,2), ,(x2(0),y2(0))=(0,0),,以及(x3(0),y3(0))=(-2,-2), 以及满足一定初始速度时的轨迹。