s=[0,0]; for i=1:1000; if s(1)<5 s=[s(1)+1/i,i]; else s=s; end end k=0; for j=1:82 k=k+1/j; end