function X = pspust_crtanje2D(A,b,P,max1,tol)


[m,n]=size(A);
hold on
plot(P(1),P(2),'o');

for k=1:max1
    r=A*P-b;
    ind=mod(k-1,m)+1;
    d=zeros(1,n)';
    d(ind)=1;

    if(norm(r,inf)<tol)
        plot(X(1),X(2),'ro'); 
        disp('Broj potrebnih iteracija: ');
        disp(k);
        break;
    end
        
    L=-r'*d/((A*d)'*d);
       
    X=P+L*d;
    pause 
    plot([X(1),P(1)],[X(2),P(2)]);
    P=X;
end

hold off
