main.m

利用四阶runge-kutta法

clear; 
[t,y]=ode45('fun',[0,35],[0,680,5000,0]); 
X=y(:,1);
Vx=y(:,2);
Y=y(:,3);
Vy=y(:,4);
subplot(1,3,1),plot(X,Y),axis([0,inf,0,5000]),grid on,xlabel('X'),ylabel('Y'),title('弹道曲线');
subplot(1,3,2),plot(t,Vx),xlabel('t'),grid on,ylabel('Vx'),title('水平速度曲线');
subplot(1,3,3),plot(t,Vy),xlabel('t'),grid on,ylabel('Vy'),title('垂直速度曲线');