function xdot = lorenzeq(t,x)
%LORENZEQ 本函數用來描述Lorenz方程。
%給定時間變量t和狀態變量x, 則由此函數可以求出導數xdot

%Designed by Prof D Xue (c) 2000
xdot=[-8/3*x(1)+x(2)*x(3); 
      -10*x(2)+10*x(3); 
      -x(1)*x(2)+28*x(2)-x(3)];