(* *) (* Numerical Mathematics and Computing, Fifth Edition *) (* Ward Cheney & David Kincaid *) (* Brooks/Cole Publ. Co. *) (* (c) 2003 *) (* *) (* File: lu_factor *) (* LU Factorization *) MatrixForm[a = {{6,-2,2,4}, {12,-8,6,10}, {3,-13,9,3}, {-6,4,1,-18}}] {lu, perm, cond} = LUDecomposition[a]; MatrixForm[lu] MatrixForm[perm]