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