myxyz.m

来自「太阳—地月系平衡解,5个平衡解」· M 代码 · 共 24 行

function q=myxyz(p)
x=p(1);y=p(2);z=p(3);

%m2与m2质量比mu
global mu;
%mu=1/(328900.5+1);
%一天文单位r
r=1;

d=[x;y;z];
%M1:
s=[mu*r;0;0];
%M2:
e=[-(1-mu)*r;0;0];
R1=d-s;
R2=d-e;

r1=sqrt(R1(1)^2+R1(2)^2+R1(3)^2);
r2=sqrt(R2(1)^2+R2(2)^2+R2(3)^2);

q(1)=x-(1-mu)*(x-mu)/r1^3-mu*(x+1-mu)/r2^3;
q(2)=y*(1-(1-mu)/r1^3-mu/r2^3);
q(3)=-z*((1-mu)/r1^3+mu/r2^3);