% Ax = b lin. e.r. megoldasa Gauss-Seidel-iteracioval
%
function x = Gauss_Seidel_it(A,b,maxiter)
S = tril( A )  % S = L + D 
T = S - A      % T = -U = S - A
x = zeros(size(b)); 	% kezdovektor  (x = 0)
for i = 1:maxiter, 	
x= S\(T*x+b);
end
