f = zeros(1,1994); f(1) = 1; for n = 2 : 1994 x = n; if rem(x,2) == 0 x = x / 2; f(n) = f(x); else x = x-1; f(n) = f(x)+1; end end max(f)