====== Math 753/853 LU decomp ====== ===== Gaussian elimination ==== % Reduce Ax=b problem to ===== LU factorization with out pivoting ====