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