Algorithms, Linear Equation.
belonging to U (transformed matrix) are irrelevant to the calculation of solutions. import numpy as np def gaussElim(a,b): n = len(b) # Elimination phase for k in range(0,n-1): for i in range(k+1,n): if a[i,k] != 0.0: #if not null define λ lam = a [i,k]/a[k,k] #we calculate the new row of the