例6-7(1).m

来自「上传最新Matlab7基础与提高,与各位爱好者分享」· M 代码 · 共 7 行

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;]       	%运动方程