Added a gamma function implementation based directly on specfun Fortran code.