% Ax = b lin. e.r. megoldasa Gauss-Seidel-iteracioval
%
function x = Gauss_Seidel_reltol(A,b,maxiter)
tol = 1.e-8; 	% eloirt hibahatar (10^-8)
relerr = inf;	 % nagy relativ hibat allitunk be kezdetben
S = tril( A )  % S = L + D 
T = S - A      % T = -U = S - A
x = zeros(size(b)); 	% kezdovektor  (x = 0)
it=0;
while (relerr > tol) & (it< maxiter) 	% maxiterig vagy az eloirt hiba elereseig iteraljon
it = it+1; 	% noveljuk az iteracioszamlalot
x_uj = S \ (T*x+b);
relerr = norm(x_uj-x,inf)/norm(x_uj,inf);
x=x_uj;
end
