def diff_lorenz(value_array, sigma, beta, rho):
diff_array = np.zeros(<NUM_LIT:3>)<EOL>diff_array[<NUM_LIT:0>] = sigma * (value_array[<NUM_LIT:1>]-value_array[<NUM_LIT:0>])<EOL>diff_array[<NUM_LIT:1>] = value_array[<NUM_LIT:0>] * (rho - value_array[<NUM_LIT:2>]) - value_array[<NUM_LIT:1>]<EOL>diff_array[<NUM_LIT:2>] = value_array[<NUM_LIT:0>] * value_array[<NUM_LIT:1>] - beta * value_array[<NUM_LIT:2>]<EOL>return diff_array<EOL>
The Lorenz attractor differential equation

    :param value_array: 3d array containing the x,y, and z component values.
    :param sigma: Constant attractor parameter
    :param beta: FConstant attractor parameter
    :param rho: Constant attractor parameter

    :return: 3d array of the Lorenz system evaluated at `value_array`