C-N解偏微分方程的程序。解的是定步长的抛物型偏微分方程： du/dx - a * d2u/dx2 = 0 在程序中可以更改 a 的值以实现不同系数的解。在循环中改变 a 的值以实现变系数

Rossler equation A. Wolf, J. B. Swift, H. L. Swinney and J. A. Vastano, "Determining Lyapunov Exponents from a Time Series," Physica D, Vol. 16, pp. 285-317, 1985.-Rossler equation dx =-y- z dy = x+ a*y dz = b+ z*(x-c)