lorenz-system

Runge Kutta constants diverging for Lorenz system? Question: I’m trying to solve the Lorenz system using the 4th order Runge Kutta method, where dx/dt=a*(y-x) dy/dt=x(b-z)-y dx/dt=x*y-c*z Since this system doesn’t depend explicity on time, it’s possibly to ignore that part in the iteration, so I just have dX=F(x,y,z) def func(x0): a=10 b=38.63 c=8/3 fx=a*(x0[1]-x0[0]) fy=x0[0]*(b-x0[2])-x0[1] …