clc
%% Jacobi iteráció demó ismert megoldású rendszerekre
display(' *** Véletlen A mátrixokkal már általában nem fog konvergálni az iteráció!')
for n = 10:10:100
  display(n)
  
  A       = rand(n);              % Az elemek 0 és 1 közti véletlen számok
  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

  k_max   = 10*n;                 % Legföljebb 10*n 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   % Ez  csak ritkán teljesül!
endfor

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