function b=forward(L,b)
% b=forward(L,b)
% Vorw"artseinsetzen b := L\b, 
% Algorithmus aus der Vorlesung
%

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

% unteres Dreieckssystem 
% b(1)=b(1) ist sowieso klar 

% nun von 2 bis n
for i=2:n 
  b(i)=(b(i)-L(i,1:i-1)*b(1:i-1)); 
  end


