float x = 0, y = 0, z = 0, dx, dy, dz, dt = 0.001, t, a = 10, b = 28, c = 8/3.0;
for (t = 0; t<10000; t = t + dt)
{
dx = a * (y - x)*dt;
dy = (x*(b - z) - y)*dt;
dz = (x*y - c * z)*dt;
x = x + dx;
y = y + dy;
z = z + dz;
}
如何判断洛伦兹方程在所有不同初值x,y,z和参数a,b,c组合下,是混沌还是稳定?