function x=success2b(x1,N)
%calcola al piu' N termini della successione
%        / x(n)/2+1 se x(n) pari
%  x(n) <
%        \ 3*x(n)-1 se x(n) dispari
% partendo da x(1) finche' x(n)=14
%
x(1)=x1;
n=1;
while (x(n)~=14)
   if rem(x(n),2)==0
      x(n+1)=x(n)/2+1;
   else
      x(n+1)=3*x(n)-1;
   end
   if n>N
      break;
   end
   n=n+1;
end
