p=0; for i=1:100 p=p+(1/i); if p>=5 break end end p=0; for i=1:82 p=p+(1/i); end