Rabinovich Fabrikant 平衡點

編輯

> alpha := 1.1;

                             1.1

> gamma := .87;

                             0.87

> > > > q1 := sqrt(1+sqrt(1-gamma*alpha*(1-3*gamma/(4*alpha))));

                    1.3347123182858183570

> q2 := sqrt(1-sqrt(1-gamma*alpha*(1-3*gamma/(4*alpha))));

                    0.46748585798513339859

> s1 := [q2, -alpha/q2, 1-(1-gamma/alpha)*q2^2];

      [0.46748585798513339859, -2.3530123557983251267, 
        0.95430463972208895291]

> s2 := [-q2, -alpha/q2, 1-(1-gamma/alpha)*q2^2];

      [-0.46748585798513339859, -2.3530123557983251267, 
        0.95430463972208895291]

> s3 := [q1, -alpha/q1, 1-(1-gamma/alpha)*q1^2];

      [1.3347123182858183570, -0.82414763460993508052, 
        0.62751354209609286530]

> s3 := [-q1, -alpha/q1, 1-(1-gamma/alpha)*q1^2];

      [-1.3347123182858183570, -0.82414763460993508052, 
        0.62751354209609286530]

>