apolloeq.m

来自「matlab的大量实例代码」· M 代码 · 共 8 行

function dx=apolloeq(t,x)
mu=1/82.45; mu1=1-mu; 
r1=sqrt((x(1)+mu)^2+x(3)^2); r2=sqrt((x(1)-mu1)^2+x(3)^2);
dx=[x(2); 
    2*x(4)+x(1)-mu1*(x(1)+mu)/r1^3-mu*(x(1)-mu1)/r2^3;
    x(4); 
    -2*x(2)+x(3)-mu1*x(3)/r1^3-mu*x(3)/r2^3];