kernelopts(printbytes=false):
Digits:=15:
n:=10000000:
hn:=0:
for k from n by -1 to 1 do
    hn:=hn+1.0/k
end:
printf("gamma(%d) = %.14f\n",n,hn-log(n));