www.pudn.com > matlab8.zip > DYDt2.m

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