Résoudre un système linéaire avec numpy

https://fr.moonbooks.org/Articles/R%C3%A9soudre-un-syst%C3%A8me-d%C3%A9quations-lin%C3%A9aires-avec-Numpy/

import numpy as np
a = np.array([[2,3,3,1], [-4,-6,3,2], [-1,1,1,1], [-2,-1,1,1]]) #matrice
b = np.array([15,3,5,1]) #second membre
x = np.linalg.solve(a, b)