例6-7(1).m

这是一个MATLAB7.0基础与提高例题的所有源码

function rdot=e31a(t,r)
c=0.01;
g=9.81;
m=1;                                          	%空气阻力系数及重力加速度（m/s^2）
vm=sqrt(r(3)^2+r(4)^2);                        	%速度大小
rdot=[r(3);r(4);-c*vm*r(3)/m;-c*vm*r(4)/m-g;]       	%运动方程