解决方法二（数值解）-Matlab微分技术

解决方法二（数值解）：1.编写m文件eq1.m如下： function dy=eq1(x,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=1/5sqrt(1+y(1)^2)/(1-x); 2.设定x0=0，xf=0.9999，编写主程序ff6.m如下： x0=0，xf=0.9999 [x,y]=ode15s('eq1',[x0 xf],[0 0]); plot(x,y(:,1),'b.') hold on y=0:0.01:2; plot(1,y,'b')。结论：导弹大致命中（1，0.2）处的目标乙舰。将y1=y，y2=y1'，将方程（3）转化为一阶微分方程组。