function dx = LorenzDifEqn1(t,x); 
% Lorenz系統微分方程
% 方程如下: 
% dx = -σ*(x-y)
% dy = -x*z+r*x-y
% dz = x*y-b*z
dx_1 = -(8/3)*x(1)+x(2)*x(3);
dx_2 = -10*x(2)+10*x(3);
dx_3 = -x(2)*x(1)+28*x(2)-x(3);
dx = [dx_1;dx_2;dx_3];