function [U,L] = lurowvect(A)
%  Gaussian Elimination with pivoting and
%  back substitution on n x n A and rhs b
% 
[n,n] = size(A);
for k=1:n-1
  A(k+1:n,k) = A(k+1:n,k)/A(k,k); 
  for j=k+1:n
    A(j,k+1:n) = A(j,k+1:n) - A(j,k)*A(k,k+1:n);
  end
end
U = triu(A);
L = eye(n)+tril(A,-1);