Easy-to-use UMFPACK interface

The following functions can be used for LU decompositions and solving equation systems:

spsolve(A, b)	Solve the sparse linear system Ax=b, where b may be a vector or a matrix.
splu(A)	Compute the LU decomposition of a sparse, square matrix.
UmfpackLU(A)	LU factorization of a sparse matrix.