8楼: Originally posted by WF1987 at 2011-12-16 16:36:38
=ode45('appollo2',,);
plot(t,y(:,3));figure;plot(t,y(:,4))
上为调用程序
%……………………………………………………………………
function dx=appollo2(t,x)
mu=1/82.45*(x(2)<0)-1*(x(2)>0);  %这是关键 ...