sum = 0;
Nmax = 0;
for j = 1:100,
  if ( sum < 5 )
    sum = sum + 1/j;
    Nmax = j;
  end
end
sum = sum - 1/Nmax 
Nmax = Nmax - 1