a=0; sum=0; for j=1:1000; a=a+(1/j); if a<5; N=max(j); end end for i=1:N; sum=sum+(1/i); end