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