Question

matlab解微分方程组并画图

𝑑𝜃/𝑑𝑧=−（1/𝑅）[n−(1−n)(2cos𝜃+ cos^3 𝜃)]
𝑑𝑅/𝑑𝑧=−1/tan𝜃
初始R=1，𝜃=90°，n=0.9，要求以R为横坐标，Z为纵坐标画图。
最好附上详细过程，谢谢。

Answer 1

http://zhidao.baidu.com/question/205690281.html

function dx=appollo(t,x)
mu=1/82.45;
mustar=1-mu;
r1=sqrt((x(1)+mu)^2+x(3)^2);
r2=sqrt((x(1)-mustar)^2+x(3)^2);
dx=[x(2)
2*x(4)+x(1)-mustar*(x(1)+mu)/r1^3-mu*(x(1)-mustar)/r2^3
x(4)
-2*x(2)+x(3)-mustar*x(3)/r1^3-mu*x(3)/r2^3];
--------------------------------------------------------------------------------------
x0=[1.2;0;0;-1.04935751];%x0(i)对应与xi的初值
options=odeset('reltol',1e-8);%该命令的另一种写法是options=odeset;options.reltol=1e-8;
tic
[t,y]=ode45(@appollo,[0,20],x0,options);%t是时间点，y的第i列对应xi的值，t和y的行数相同
toc
plot(y(:,1),y(:,3))%绘制x1和x3，也就是x和y的图形
title('Appollo卫星运动轨迹')
xlabel('X')
ylabel('Y')
请采纳。

追问

这也对不上啊