clc
display(' *** Diagonális domináns véletlen A mátrixokra néha egész gyors a konvergencia!')
%% Jacobi iteráció demó ismert megoldású rendszerekre
for n = 100:100:1000
  display(n)
  
  A = randi([0 10],n);          % Az elemek 0 és 10 közti egészek
  A = A + 10*diag(sum(A'));     % Így "nagyon erõsen" diagonális domináns az A mátrix
  xpontos=ones(n,1);            % Legyen a pontos megoldás a csupa 1-es vektor
  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 = n;                    % Legföljebb n iterációs lépést végzünk
  d_min = 10^(-15);
  
  [x_koz, k, hiba1] = Jacobi_iteracio(A,b,x0,k_max,d_min);
  k
  gyarihiba = norm(xpontos - xgyari)
  iteraciohiba = norm(xpontos - x_koz)
  Jacobi_pontosabb = iteraciohiba < gyarihiba   % MINDEN esetben ez teljesül!
endfor

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