1> a = [4 1 0; 1 -1 5; 2 2 -3]
a =

   4   1   0
   1  -1   5
   2   2  -3



2> b = fattlu(a)
b =

   4.00000   1.00000   0.00000
   0.25000  -1.25000   5.00000
   0.50000  -1.20000   3.00000



3> u = triu(b)
u =

   4.00000   1.00000   0.00000
   0.00000  -1.25000   5.00000
   0.00000   0.00000   3.00000



4> l = tril(b)
l =

   4.00000   0.00000   0.00000
   0.25000  -1.25000   0.00000
   0.50000  -1.20000   3.00000



5> l = tril(b,-1)
l =

   0.00000   0.00000   0.00000
   0.25000   0.00000   0.00000
   0.50000  -1.20000   0.00000



6> 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



7> l * u
ans =

   4   1   0
   1  -1   5
   2   2  -3



8> a
a =

   4   1   0
   1  -1   5
   2   2  -3



9> c =[0 4; 1 -2]
c =

   0   4
   1  -2



10> d = fattlu(c)
d =

     0     4
   Inf  -Inf



11> c
c =

   0   4
   1  -2



12> % usiamo una matrice di permutazione p = [0 1 ; 1 0]
    % c * p fa scambiare le colonne di c
    % p * c fa scambiare le righe di c

12> e = c * [0 1; 1 0]
e =

   4   0
  -2   1



13> e = [0 1; 1 0] * c
e =

   1  -2
   0   4



14> c
c =

   0   4
   1  -2



15> fattlu(e)
ans =

   1  -2
   0   4



16> [l, u, p] = lu(c)
l =

   1   0
   0   1

u =

   1  -2
   0   4

p =

Permutation Matrix

   0   1
   1   0



17> quit