function [x] = usol(U,b);
%
% solve the system Ux=b, where U is assumed to be nonsingular upper
% triangular (the code looks only at the upper triangular
% elements).
%