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