function [L,b] = matrix_L(n) 

% Dieses Programm erzeugt eine zuf"allige linke untere Dreiecksmatrix
% der Dimension n x n und einen Vektor der Dimension n
% Der Aufruf ist [L,b] = matrix_L(n)

c = 20*(rand(1)-0.5);

L = rand(n,n)-0.5*ones(n,n);

for i=1:n
  for j=i+1:n
    L(i,j) = 0;
  end
end

for i=1:n
  if L(i,i) == 0
    L(i,i) = 0.4711;
  end
end

L = c*L;

b = rand(n,1) - 0.5*ones(n,1);
b = c*b;




