VBA implementation of LU decomposition with partial pivoting, forward and back substitution, matrix multiplication