matlab求解微分方程

y’‘-0.01(y’)^2+2y1=sin(t),y(0)=0,y’(0)=1,0<t<5,做y的图
新建M函数
function dy=weifen2(t,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=0.01y(2)^2-2y(1)+sin(t);
输入命令
[t,y]=ode15s(‘weifen2’,[0,5],[0 1]);
plot(t,y)

1.绘制飞船轨迹图
新建M函数
function dy=weifen3(t,y)
dy=zeros(4,1);
dy(1)=y(3);
dy(2)=y(4);
dy(3)=2y(4)+y(1)-(1-1/82.45)(y(1)+ 1/82.45)/((y(1)+1/82.45)^2+y(2)2)^{(3/2)-(1/82.45)*(y(1)+1/82.45-1)/((y(1)+1-1/82.45)}2+y(2)²⁾(3/2);
dy(4)=-2y(3)+y(2)-(1-1/82.45)*y(2)^{2/((y(1)+1/82.45)}2+y(2)²⁾(3/2)-(1/82.45)*y(2)/((y(1)+1-1/82.45)^2+y(2)2)^(3/2);
输入命令
[t,y]=ode15s(‘weifen3’,[0,10],[1.2 0 0 -1]);
plot(t,y)