harj32.m

% Funktio saa argumentin n, palauttaa muuttujan pal
function pal = harj32( n )
    x = n;
    while n > 1
        if mod(n,2) == 0 
            % Parillinen
            n = n/2;
        else
            % Pariton
            n = 3*n+1;
        end
        x(end+1) = n; % Tallenna uusi muuttujan n arvo vektoriin
    end

    disp( x )        % Tulosta vektori x ruudulle
    pal = length(x); % Aseta palautusarvo
end