my $log_nfact =
  $n * log ($n) - $n + 0.5 * (log(2*PI*$n + 1/3));