function pT = projekcija(A,b,T)
%pT = projekcija(A,b,T) vrne krajevni vektor pT 
%pravokotne projekcije rT na afin podprostor 
%dan s poddolocenim sistemom Ax=b.

pT = pinv(A)*(b - A*T) + T;