function b=BACKWARD(U,b,unit)
% function b=BACKWARD(U,b)
% R"uckw"artsseinsetzen b := U\b, Algorithmus aus der Vorlesung
%

% Dimension der Matrix
n=size(U,1);


% oberes Dreieckssystem 
% erster, d.h. n-ter Schritt
b(n)=b(n)/U(n,n); 

% und nun dei anderen
for i=n-1:-1:1 
  b(i)=(b(i)-U(i,i+1:n)*b(i+1:n))/U(i,i); 
  end; 


