First draft over LU factorization

Need to be improved…
Lacks numerical experiments.