added matrix subtraction/negation and MAPM to map a function across the elements of matrices
Added LU decomposition functions, including inverse of symmetric matrix