gusucode.com > 《精通matlab 7》包括配套光盘 > 配套光盘里的程序附件/Ex-22/Ex2204.m

    %例22-4：求描述某非刚性体的运动方程的微分方程 
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')