有人用xppaut或matlab画过Neimark–Sacker分叉图吗，微分方程组是如下：
    dxi/dt=beta1*(-x(i)+alpha/(1+(z(i))^n));
    dyi/dt=beta2*(-y(i)+alpha/(1+(x(i))^n));
    dzi/dt=beta3*(-z(i)+alpha/(1+(y(i))^n)+K*s(i)/(1+s(i)));
    dsi/dt=-kS0*s(i)+kS1*y(i)-Eta*(s(i)-Sext)
这里i=1,2,3,共有12各个方程组。