www.pudn.com > 200120237.rar > apolloeq.m

function dx=apolloeq(t,x) 
%APOLLOEQ 函数用来描述 Apollo 卫星的数学模型，它可以直接用于 
%  常微分方程的数值求解。 
 
%Designed by Prof D Xue (c) 2000 
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];