# Evgenii B. Rudnyi, http://MatrixProgramming.com

from numpy import *
from time import *
from scipy.linalg import *

dim = 2000
t = clock()
mat = random.rand(dim,dim)
print 'time for matrix', clock() - t
print mat.shape
print mat[:10,:10]
t = clock()
lu = lu_factor(mat)
print 'time for LU', clock() - t