help lugauss
a=[4 1 0;1 -1 5;2 2 -3]
b=lugauss(a)
u=triu(b)
tril(b)
tril(b,-1)
l=tril(b,-1)+eye(3)
quit
