Implementing a Pecora-Carroll Lorenz system into python?

Question

I have the following system which is two Lorenz systems coupled together:

def coupled_dynamics_Lorenz(x,t,sigma=10.0, beta=8.0/3.0, rho=28.0, k=0.0):
    xdot