using sn6l
n=20000000000
@printf("gamma(%d) = %.14f\n",n,Hpar(1,n)-log(n))