Algorithms, Linear Equation.
b=np.array([1.0,1.0,-5.0]) aOrig = a.copy() # save original matrix A bOrig = b.copy() #save original vector b x = gaussElim(a,b) #print A transformed for check print(a) print("x =\n",x) #det = np.prod(np.diagonal(a)) #determinant #print("\ndet =",det) #check result and numerical precision