# LUDecomposition D¶

LUDecomp(A) computes a LU decomposition of A using the algorithm of Crout. LU contains both triangular matrices; Perm is the permutation used for partial pivoting and Pivots yields the used pivots.
LUInverse(A) computes the inverse of A using a LU decomposition.
LUSolve(LU, Perm, B) uses a previously computed LU decomposition to solve a linear system with right hand side B. LU and Perm are as given by LUDecomp.