znxpydfun.m

来自「研究阻尼斜抛运动的MATLAB程序」· M 代码 · 共 6 行

function ydot=znxpydfun(t,y,flag,b,p,m)
ydot=[y(2);
   -b/m*y(2)*(y(2).^2+y(4).^2)^(p/2);
   y(4);
   -9.8-b/m*y(4)*(y(2).^2+y(4).^2)^(p/2)];