F=ode45(@(t,y)func1(t,y),[0,10],[5;0;0;300]); plot(F.y(3,:),F.y(4,:)); xlabel('x轴'),ylabel('y轴') function dF=func1(t,F) %F=[dx,dy,x,y]; m=1; k=0.2; g=9.8; dF=[-k./m.*F(1);-g-k./m.*F(2);F(1);F(2)]; end