dydt2.m

来自「学matlab还有比这个更好的吗?我千兴万苦才找到的,这是一个学习matlab的」· M 代码 · 共 10 行

function Yd=DYDt2(t,Y,flag,G,ME)
% flag	按ODE文件格式规定，必须是第三输入宗量。对它的赋值由ode45指令自动产生。
% 第4、5宗量是被传递的参数
switch flag
case ''
	X=Y(1:2);V=Y(3:4);r=sqrt(sum(X.^2));Yd=[V;-G*ME*X/r^3];
otherwise
  error(['Unknown flag ''' flag '''.']);
end