p=0;
for j=1:82;
p=p+(1/j);
end