% Jacobi-iteracio hibakuszob beallitasaval
% Az X tomb i-edik oszlopaban lesz az i. iteracio
function x = Jacobi_it_reltol(A,b,maxiter)
tol = 1.e-8; 	% eloirt hibahatar (10^-8)
relerr = inf;	 % nagy relativ hibat allitunk be kezdetben
it = 0; 	% elkezdjuk szamolni a lepeseket
D = diag( diag(A) ); 	
T = D - A;	 
x = zeros(size(b)); 	% kezdovektor  (x = 0)
while (relerr > tol) & (it< maxiter) 	% maxiterig vagy az eloirt hiba elereseig iteraljon
it = it+1; 	% noveljuk az iteracioszamlalot
x_uj = D \ (T*x+b);
relerr = norm(x_uj-x,inf)/norm(x_uj,inf);
x=x_uj;
end
%x= X(:, it);
