chaoszhao_ly.m

we determin a chaos mode and alyalize the time application to find if its really chaos and will have

function dX =chaoszhao_ly(t,X)
global a1 a2 a3 b1 b2 b3 c1 c2 c3 M N K;
x=X(1); y=X(2); z=X(3);
Y = [X(4), X(7),  X(10);
     X(5), X(8),  X(11);
     X(6), X(9),  X(12)];

dX = zeros(12,1);

dX(1) = a1*x*(1-x/M)-a2*(y-K+z)-a3*x;
dX(2) = b1*(y-K)-b2*x*y-b3*z;
dX(3) = c1*z*(c2*x-c3);

Jaco = [a1*(1-2*x/M)-a3     -a2        -a2;
        -b2*y              b1-b2*x     -b3;
         c1*c2*z          0      c1*(c2*x-c3)];       
dX(4:12) = Jaco*Y;