Factorisation LU

import pprint
import numpy
import scipy
import scipy.linalg   # SciPy Linear Algebra Library

A = numpy.array([ [7, 3, -1, 2], [3, 8, 1, -4], [-1, 1, 4, -1], [2, -4, -1, 6] ])
P, L, U = scipy.linalg.lu(A)

print ("A:")
pprint.pprint(A)

#print ("P:")
#pprint.pprint(P)

print ("L:")
pprint.pprint(L)

print ("U:")



pprint.pprint(U)