#NEWLINE# def diff_lorenz( value_array , sigma , beta , rho ) : #NEWLINE# #INDENT# 'The Lorenz attractor differential equation\n\n    :param value_array: 3d array containing the x,y, and z component values.\n    :param sigma: Constant attractor parameter\n    :param beta: FConstant attractor parameter\n    :param rho: Constant attractor parameter\n\n    :return: 3d array of the Lorenz system evaluated at `value_array`\n\n    ' #NEWLINE# diff_array = np . zeros ( 3 ) #NEWLINE# diff_array [ 0 ] = ( sigma * ( value_array [ 1 ] - value_array [ 0 ] ) ) #NEWLINE# diff_array [ 1 ] = ( ( value_array [ 0 ] * ( rho - value_array [ 2 ] ) ) - value_array [ 1 ] ) #NEWLINE# diff_array [ 2 ] = ( ( value_array [ 0 ] * value_array [ 1 ] ) - ( diff_array * value_array [ 2 ] ) ) #NEWLINE# return diff_array