tspan = [0 12];
% 在同一目錄下，計算方程數值解。給出時間區間和初始值
y0 = [-1 1 0];
[t,Y] = ode45('odet3',tspan,y0);
% 采用ode45算法
plot(t,Y(:,1),'-',t,Y(:,2),'-.',t,Y(:,3),'.')
% 繪制計算結果并標注，如圖3-15
legend('Y1','Y2','Y3')