function [V] = itprjr(MZ,SV,MXV,epp)
   
   V = SV;
   
   MZp = pinv(MZ);
   RK = 1;
   
   vh = MZp*V(:,RK);
   
   while (rank([V vh],epp) > RK) & (RK<MXV) 
      V = [V vh];
      RK = RK+1;
      vh = MZp*V(:,RK);
   end
end
