例4-7-1.m

这是一个MATLAB7.0基础与提高例题的所有源码

function sequence=collatz(n)
%柯雷茨（Collatz）猜想
%   对于任意的正整数n:
%   如果n是偶数则除以2
%   如果是奇数则乘3加1
%   重复计算，直至结果为1%  

sequence = n;
next_value = n;
while next_value > 1
    if rem(next_value,2)==0
        next_value = next_value/2;
    else
        next_value = 3*next_value+1;
    end
    sequence = [sequence, next_value];
end