I know that for: a+ bx1+ cx2=0, a+ 2bx1+ cx2=8 I can solve the equation using linalg:
g = np.array([[a,b], [a,2b]])
w = np.array([0,8])
solution = np.linalg.solve(g, w)
But how can I solve the following equations using np: a+ bx1+ cx2 + dx1^2+ 3ex2^2=0, a+ 2bx1+ 2cx2 + dx1^2+ ex2^2=8