import numpy as np A=[[1,1,1,1],[1,1,1,1],[1,1,1,1],[1,1,1,1],[1,1,0,0]] b=[1,1,1,1,1] X=np.linalg.lstsq(A, b, rcond = -1) print (X[0]) print(A @ X[0])

import numpy as np A=[[1,1,1,1],[1,1,1,1],[1,1,1,1],[1,1,1,1],[1,1,0,0]] b=[1,1,1,1,1] X=np.linalg.lstsq(A, b, rcond = -1) print (X[0]) print(A @ X[0])