Fortran Wiki
log_gamma