Solution of overdetermined equation

Do your homework today , It is required to solve a system of overdetermined equations , I thought I could use it directly sympy Library linsolve() Just get the function done , The result cannot be obtained , I guess it's because of the error

import numpy as np
from sympy import *

A = np.array([[-0.75, 0.25, 0, 0], 
            [0.2, -0.8, 0.25, 0], 
            [0.12, 0.12, -0.8, 0.25], 
            [0.43, 0.43, 0.55, -0.25], 
            [1.00, 1.00, 1.00, 1.00]])
b = np.array([0, 0, 0, 0, 1.0])
A = np.column_stack((A, b))
A = Matrix(A)

A = A.rref()  // The augmented matrix is transformed into a row stepped matrix 

A = np.array(A[0])
print(A)

In fact, it can also be used np.linalg.lstsq() Method , But this method is an approximation , Can't judge the solution of the equations , But you can use **np.linalg.matrix_rank()** First judge the rank of the augmented matrix, and then use this method to solve

A = np.array([[-0.75, 0.25, 0, 0], 
            [0.2, -0.8, 0.25, 0], 
            [0.12, 0.12, -0.8, 0.25], 
            [0.43, 0.43, 0.55, -0.25], 
            [1.00, 1.00, 1.00, 1.00]])
b = np.array([[0], [0], [0], [0], [1.0]])
x = np.linalg.lstsq(A, b, rcond=None)[0]