b = 0; for j = 1:100 while b + 1/j < 5 b = b+1/j; end end