dydt2.m

精通matlab

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