% 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 = 1; 	% elkezdjuk szamolni a lepeseket
D = diag( diag(A) ); 	
T = D - A;	 
X(:,1) = zeros(size(b)); 	% kezdovektor  (x = 0)
while (relerr > tol) & (it< maxiter) 	% maxiterig vagy az eloirt hiba elereseig iteraljon
X(:,it+1) = D \ (T*X(:,it)+b);
relerr = norm(X(:,it+1)-X(:,it),inf)/norm(X(:,it+1),inf);
it = it+1; 	% noveljuk az iteracioszamlalot
end
x= X(:, it);