function [U,b] = matrix_U(n) 

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

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

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

for i=2:n
  for j=1:n-1
    U(i,j) = 0;
  end
end

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

U = c*U;

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




