3X+1问题

%从任意一个正整数开始，重复对其进行下面的操作：
%如果这个数是偶数，把它除以 2 ；
%如果这个数是奇数，则把它扩大到原来的 3 倍后再加 1
%你会发现，序列最终总会变成 4, 2, 1, 4, 2, 1, … 的循环。
X=input('请输入一个数：');
count=0;
number=0;
while(1)
    if rem(X,2)==0
        X=X/2;
        count=count+1;
    else
        X=3*X+1;
        count=count+1;
    end
    if X==1
        number=number+1;
    end
    fprintf('%d  ',X);
    if rem(count,7)==0
        fprintf('\n')
    end
    if number==3
        break;
    end
end
fprintf('\n')