Help with LU factorization

[ettersi]

You can get the n \times n identity matrix using Matrix{eltype(A)}(I,(n,n)).