Linear Algebra, Python, Numpy, Matrix Decomposition, Eigenvectors.
information about how the PLU decomposition works can be found here. import scipy import scipy.linalg a = np.array([[1,2,4],[3,8,14],[2,6,13]]) p,l,u = scipy.linalg.lu(a) print("Lower Matrix") print(l) print("\n") print("Upper Matrix") print(u) Lower Matrix [[ 1. 0. 0. ] [ 0.66666667 1. 0. ] [