Clean code.
Added LU factorization from rosetta code.