Fix bug cases reported by Rocco Romeo:
* Add another approximation for x very small.
* Add code to handle cases where tgamma(a) would overflow.
* Add tests for the above.