Simple implemention of LU decomposition
Simple subroutines for Vec and Mat
Perhaps i should refer to BLAS and LAPACK later