Use a LU decomposition to solve the system A.x = b
[xsol] = sp_lusolve(A,b])
the A real sparse matrix of the system to be solved.
the b real vector of the system to be solved.
the solution vector of the system.
Solve the system A.x = b via a LU decomposition.