
function  f_matlab
zahl_max = input('zahl_max')
for zahl = 1 : zahl_max
    fprintf('%5d %1d\n', zahl, pruefe(zahl))
end

    function janein = pruefe(zahl)
    qsumme = 0;
        iteration = 0 ;
        while (zahl ~= 1) & (zahl ~= 4)   % stop, wenn z=1 oder z=4
            qsumme = 0;
            while zahl > 0
                ziffer = rem(zahl, 10);
                qsumme = qsumme + ziffer*ziffer;
                zahl = fix( zahl / 10 );
            end
            iteration = iteration+1;
            zahl = qsumme;
        end
        if qsumme == 1;
            janein = 1;
        else 
            janein = 0;
        end
    end
    
end

