factorization code snippets in matlab