function x=treennepiuuno(x1)
% x=TREENNEPIUUNO(x1)
% calcoloa la successione di collatz definita come
%          | x(n)/2 se x(n) e' pari 
%   x(n+1)= <
%          | 3x(n)+1 se x(n) e' dispari
% partendo da x(1)=x1 e fermandosi non appena x(n)=1

x(1)=x1;
n=1;
while x(n)~=1
    if rem(x(n), 2)==0
        x(n+1)=x(n)/2;
    else
        x(n+1)=3*x(n)+1;
    end
    n=n+1;
end

