Skip to content

lorenz63.data

dyn_lorenz63(t, x, sigma=10.0, rho=28.0, beta=8.0 / 3)

Lorenz-63 dynamical model.

Source code in contrib/lorenz63/data.py 
13
14
15
16
17
18
19
def dyn_lorenz63(t, x, sigma=10., rho=28., beta=8./3):
    """ Lorenz-63 dynamical model. """
    x_1 = sigma*(x[1]-x[0])
    x_2 = x[0]*(rho-x[2])-x[1]
    x_3 = x[0]*x[1] - beta*x[2]
    dx  = np.array([x_1, x_2, x_3])
    return dx