clc
%%
disp('Gauss-Seidel iteráció demó')
disp('A megoldandó egyenletrendszer:');   % A diagonális domináns!
A = [10, 1,  2,  1;
      1, 10,  1,  0;
      1, 2, -5,  1;
      1, 1,  1, 10
    ]
b = [22; 24; -6; 46]
%%
disp('A pontos megoldása:');
xpontos = [1; 2; 3; 4]
x_gyari = A\b
%% 
disp('A megoldás közelítése Gauss-Seidel-iterációval:');
k_max = 10
[x_koz]= Gauss_Seidel_iteracio_egyszeru(A, b, k_max)

hiba = norm(xpontos - x_koz,inf)
disp('Konvergál-e az iteráció? Elfogadható-e a kapott közelítés?');
%% 
k_max = 100
[x_koz]= Gauss_Seidel_iteracio_egyszeru(A, b, k_max)
hiba = norm(xpontos - x_koz,inf)
disp('Konvergál-e az iteráció? Elfogadható-e a kapott közelítés?');