>> % Illustration of LU (see 77-83)
>> A=[3 2 -1;9 7 -5;-6 6 3];
>> b=[4 8 15]';
>> [l,u]=slu(A)
l =
     1     0     0
     3     1     0
    -2    10     1
u =
     3     2    -1
     0     1    -2
     0     0    21
>> l*u
ans =
     3     2    -1
     9     7    -5
    -6     6     3
>> b
b =
     4
     8
    15
>> x= l\b
ans =
     4
    -4
    63
>> u\x
ans =
     1
     2
     3