EDUğ a1 = [1 -1 0]'
EDUğ a2 = [2 0 -2]'
EDUğ a3 = [3 -3 3]'
EDUğ A = [a1 a2 a3]
EDUğ A(2,:) = A(2,:) + A(1,:)
EDUğ A(3,:) = A(3,:) + A(2,:)
%
% now we get an upper triagular matrix
%
EDUğ cd a:\Tcodes
EDUğ [l,u]=slu([a1 a2 a3])
%
% one step to get LU decomposition
%
EDUğ help grams
EDUğ [Q, R] = grams([a1 a2 a3])