function [x,count] = Upper(U,y)
n = size(y);
x(:,1)=y;count = 0;
for j = n:-1:2
    for i = 1:j-1
        count = count +2;
        x(i) = x(i) - U(i,j)*x(j);
    end
end
    