import numpy as np coefficients = np.array([ [3.0, 1.0], [1.0, 2.0], ]) rhs = np.array([9.0, 8.0])