function Yd=DYDt2(t,Y,flag,G,ME)
% flag	按ODE文件格式規(guī)定，必須是第三輸入宗量。對它的賦值由ode45指令自動產(chǎn)生。
% 第4、5宗量是被傳遞的參數(shù)
switch flag
case ''           %按規(guī)定：這里必須是空串。在此為"真"時，完成以下導(dǎo)數(shù)計算。
	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