%
%
%   Testen von GaussLU, Forward, Backward
%
%

% Dimension des LGS
n=100

% Matrix erzeugen
disp('Generiere Matrix und rechte Seite ...')
A=rand(n,n);
% A=hilb(n);

% Exakte Loesung des LGS [1,2, ... , n]
x_exakt=ones(n,1);
for i=1:n 
  x_exakt(i,1)=i; 
end;

% Rechte Seite bauen
b=A*x_exakt;

% LU-Zerlegung
disp('LU-Zerlegung ...')
A=gausslu(A);

% FORWARD
disp('FORWARD ...')
b=forward(A,b);

% BACKWARD
disp('BACKWARD ...')
b=backward(A,b);

% Fehlernorm
disp('Norm des Fehlers ...')
err=norm(b-x_exakt)




