In [4]:
import numpy as np
import matplotlib.pyplot as plt
import scipy.linalg as sla
%matplotlib inline
In [16]:
A = np.random.rand(10,10)
In [17]:
P, L, U = sla.lu(A)
In [25]:
plt.spy(P.dot(L))
Out[25]:
<matplotlib.image.AxesImage at 0x1071dc3c8>