function dX = qi_hyper(t,X)

 
global a b c d e f; % 變量不放入參數表中,通過命令窗口賦值

 
x=X(1); y=X(2); z=X(3);w=X(4);
% Y的四個列向量為相互正交的單位向量
Y = [X(5), X(9), X(13),X(17);
    X(6), X(10), X(14),X(18);
    X(7), X(11), X(15),X(19);
    X(8), X(12), X(16),X(20)];
% 輸出向量的初始化
dX = zeros(20,1);
 
% qi吸引子
dX(1) = a*(y-x)+y*z;
dX(2) = b*(x+y)-x*z;
dX(3) = -c*z-e*w+x*y;
dX(4) = -d*w+f*z+x*y;
 
% qi吸引子的Jacobi矩陣


Jaco = [-a  a+z y   0;
        b-z b  -x   0;
        y   x  -c  -e;
        y   x   f  -d];
         
dX(5:20) = Jaco*Y;