clc
display('Egy komolyabb feladat, diagonális domináns A-val:')
n = 1000;
A = rand(n,n);
A = A + n*diag(ones(n,1)); % Így biztosan diagonális domináns az A
b = rand(n,1);

x0    = b;
k_max = n;
d_max = 10^(-16);

display('A megoldás közelítése Jacobi-iterációval, időméréssel:')
tic()
[x_koz,k,hiba1]= Jacobi_iteracio(A, b, x0, k_max, d_max);
Jacobi_ido = toc()
display('A behelyettesítés hibájának végtelen normája:')
r = A*x_koz -b;
behelyhiba_inf = norm(r, 'inf')
display('NEM ismerjük a pontos megoldást, nézzük a gyári eredménytől való eltérést:')
x_gyari = A\b;
gyhiba_inf = norm(x_koz-x_gyari, 'inf')

display('A megoldás közelítése a ciklussal fölírt Jacobi-iterációval, időméréssel:')
tic()
[x_koz,k,hiba1]= Jacobi_iteracio_ciklus(A, b, x0, k_max, d_max);
Jacobi_ido_ciklus = toc()                     %% a ciklikus változat LASSÚ!
display('A behelyettesítés hibájának végtelen normája:')
r = A*x_koz -b;
behelyhiba_inf = norm(r, 'inf')
display('NEM ismerjük a pontos megoldást, nézzük a gyári eredménytől való eltérést:')
x_gyari = A\b;
gyhiba_inf = norm(x_koz-x_gyari, 'inf')
display('Konvergál-e az iteráció? Elfogadható-e a kapott közelítés?')
