(phixonline)--> atom n, rn = 0, lim = power(2,iff(machine_bits()=32?53:64)) for i=1 to 10 do n = iff(i=1?2:n*n-n+1) printf(1,iff(n<=lim?"%d: %d\n":"%d: %g\n"),{i,n}) rn += 1/n end for printf(1,"sum of reciprocals: %g\n",{rn})