求matlab高人指点，解一道数学题啊！急急急！如图！

Question

第三题。用matlab解答。周一晚上要交。请高人指点啊！

Answer 1

function zd503781491

% 求解微分方程
Y0 = [1.2; 0; 0; -1.04935371];
[t,Y] = ode45(@ode, [0 100], Y0);

% 绘制轨迹图
x = Y(:, 1);
y = Y(:, 3);
plot(x,y)
axis equal

function dY = ode(t, Y)
% 阿波罗飞船运动微分方程
mu = 1 / 82.45;
lambda = 1 - mu;

x = Y(1);
dx = Y(2);
y = Y(3);
dy = Y(4);

r1 = sqrt( (x+mu)^2 + y^2 );
r2 = sqrt( (x+lambda)^2 + y^2 );

dY = [
dx;
2*dy + x - lambda*(x+mu) / r1^3 - mu*(x-lambda) / r2^3;
dy;
-2*dx + y - lambda*y / r1^3 - mu*y / r2^3;
];

Answer 2

我就不解了，时间赶不上了！

Answer 3

buzhidao