function [X]=stacionarna(A,b,P,maxiter,tol)

% Hermiteova
if ~isequal(A,A')
    error('Mora da vazi A*=A')
end

% Pozitivno definitna
[R,p] = chol(A);
if p~=0
    error('Matrica nije pozitivno definitna');
end

M=max(abs(eig(A)));
m=min(abs(eig(A)));

tau=2/(M+m);

X = P;
for k=1:maxiter   
     r = A*X - b;
     if norm(r,inf)<tol         
         return
     end
     X = X - tau*r;
end

norm(r,inf)
