- numpy - scipy

The Cat Computing Universe: PyLab

import numpy as np import scipy.linalg as la x = np.random.rand(3,3) print("Performing LU Decomposition of ", x) (P, L, U) = la.lu(x) print("P is ", P, "L is ", L, "U is ", U)