1> edit fattlu.m 2> a = [4 1 0; 1 -1 5; 2 2 -3] a = 4 1 0 1 -1 5 2 2 -3 3> b = fattlu(a) a = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 2.00000 2.00000 -3.00000 a = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 2.00000 2.00000 -3.00000 a = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 0.50000 2.00000 -3.00000 a = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 0.50000 2.00000 -3.00000 a = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 0.50000 -1.00000 -28.00000 b = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 0.50000 -1.00000 -28.00000 4> b b = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 0.50000 -1.00000 -28.00000 5> b = fattlu(a) a = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 2.00000 2.00000 -3.00000 a = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 2.00000 2.00000 -3.00000 a = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 0.50000 2.00000 -3.00000 a = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 0.50000 2.00000 -3.00000 a = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 0.50000 -1.00000 -28.00000 b = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 0.50000 -1.00000 -28.00000 6> edit fattlu.m 7> b = fattlu(a) b = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 0.50000 -1.00000 -28.00000 8> b = fattlu(a) b = 4.00000 1.00000 0.00000 0.25000 -2.00000 5.00000 0.50000 -1.00000 -28.00000 9> b = fattlu(a) b = 4.00000 1.00000 0.00000 0.25000 -1.25000 5.00000 0.50000 -1.20000 3.00000 10> l = tril(b,-1) + eye(3) l = 1.00000 0.00000 0.00000 0.25000 1.00000 0.00000 0.50000 -1.20000 1.00000 11> u = triu(b) u = 4.00000 1.00000 0.00000 0.00000 -1.25000 5.00000 0.00000 0.00000 3.00000 12> l * u ans = 4 1 0 1 -1 5 2 2 -3 13> a a = 4 1 0 1 -1 5 2 2 -3 14> quit