function mh = jorana(V,A,z,epp)

   nv = size(V,2);

   em = A*V-V*jorblk(z,nv);
   nem = norm(em);

   if nem > epp
      mh = jorana(V(:,1:nv-1),A,z,epp);
   else
      mh = V;
   end

