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