tspan = [0 12];
y0 = [0 1 1];
% 在同一目錄下，計算方程數(shù)值解。輸入時間區(qū)間和初始條件
[t,Y] = ode45('rigid',tspan,y0);
% 采用ode45算法求解方程，options為默認值
plot(t,Y(:,1),'-',t,Y(:,2),'-.',t,Y(:,3),'.')
% 繪制計算結(jié)果并標注，如圖3-14
legend('Y1','Y2','Y3')