LU Decomposition

function LUDecompose has been added which takes
a matrix A, and produces L, U and the permutation 
matrix satisfying PA = LU