function x=rwe(A, b)
% Rueckwaertseinsetzen nach Proposition 2.2, Seite 29
% loest A*x=b, wenn A eine nichtsingulaere 
% rechte obere Dreiecksmatrix ist 
[m n]=size(A); % Groesse von A holen
if m~=n 
    error('Matrix ist unsymmetrisch');
end
x=zeros(n,1);  % Loesungsvektor dimensionieren
for j=n:-1:1   % RUECKWAERTS rechnen!!!!
    x(j,1)=(b(j,1)-A(j,j+1:end)*x(j+1:end,1))/A(j,j); 
    % das ist die Formel (2.3) des Buches    
end
