%在三維空間演示一個由著名的勞倫茲非線性奇異方程 所描述的無序運動。先介紹一
%下勞倫茲方程，其形式如下：dy/dt=Ay
%方程中包含一個由函數y(t)確定的矢量和一個矩陣A。且A依賴于y，存在著下面的
%關系：A=[.8/3 ,0, y(2);0,.10,10;.y(2),28,.1]
%該程序的代碼如下：
A=[.8/3 0 .10;0 .10 10;10 28 .1];
y=[35;.10;.7];
h=2.01;
p=plot3(y(1,1),y(2,1),y(3,1),'*','erasemode','none');
axis([0 50 .25 25 .25 25]);
hold on;
for i=1:8000
    A(1,3)=y(2,1);
    A(3,1)=y(2,1);
    dita_y=h*A*y;
    y=y+dita_y;
set(p,'xdata',y(1,1),'ydata',y(2,1),'zdata',y(3,1));
    drawnow;
end