%
% Loesung eines Gleichungssystmes Ux = y,
% wobei U eine obere Dreiecksmatrix ist
%
% Bemerkung: Es werden zwei Varianten gezeigt. Man braucht jedoch nur eine!
%
% Vorlesung "Einfuerung in die Numerik", Uni Basel, FS 11
%

function x = Usolve(U, y)

% Variante 1 (ohne geschachtelte for-Schleife)
n = length(y);
x = zeros(n, 1);

x(n) = ... ;

for i = n-1 : -1 : 1
    x(i) = ... ;   
end

% Variante 2 (mit geschachtelter for-Schleife)
% n = length(y);
% x = zeros(n, 1);
% 
% for i = n : -1 : 1
%     sum = 0;
%     for k = i+1 : n
%         sum =... ;
%     end
%     x(i) = 1/...*(... - sum);
% end

end
