clc
%% Jacobi iteráció demó ismert megoldású rendszerekre
display(' *** Diagonális domináns véletlen A mátrixokra néha egész gyors a konvergencia!')
for n = 100:100:1000
  display(n)
  
  A       = rand(n);              % Az elemek 0 és 1 közti véletlen lebegőpontos számok
  A       = A + 2*diag(sum(A'));  % Így biztosan diagonális domináns az A mátrix
  xpontos = rand(n,1);            % A véletlenül választott xpontos megoldás
  b       = A*xpontos;            % Az ehhez tartozó b vektor
  
  xgyari  = A\b;             % A beépített megoldó ezt adja
  
  x0      = b;               % A kezdővektor legyen a b
  k_max   = floor(n/2);      % Legföljebb n/2 iterációs lépést végzünk
  
  [x_koz] = Jacobi_iteracio_egyszeru(A,b,k_max);
  
  gyarihiba    = norm(xpontos - xgyari)
  iteraciohiba = norm(xpontos - x_koz)
  Jacobi_jobb  = iteraciohiba < gyarihiba   % esetleg n = 100 kivételével MINDIG!
endfor

%% Kísérletezzünk tovább a különböző paraméterekkel!