function examp1_4
[t,y]=ode45(@lorenzeq,[0,100],[0;0;1e-10]); % 求解方程
plot3(y(:,1), y(:,2), y(:,3))  % 繪制三維相軌跡
axis([10 42 -20 20 -20 28]);   % 手動設置坐標系位置
function xdot = lorenzeq(t,x)  % 下面的函數直接對應于給出的微分方程
xdot=[-8/3*x(1)+x(2)*x(3);     % 這三行對應于微分方程的三個方程
   -10*x(2)+10*x(3);
   -x(1)*x(2)+28*x(2)-x(3)];