FORTRAN routines for computation of Special Functions

Occasionally I have to evaluate some obscure (at least to me) mathematical functions in Fortran.  Attempting to code these functions is no trivial task to get it right.  In many cases, these functions are series solutions requiring an unknown number of terms for computational accuracy.  A while back I found this site which has an extensive list of subroutines for evaluating mathematical functions in Fortran.  The subroutines are also optimized for to ensure the efficient calculation of an accurate answer, making their functions much better than anything I could have created.  So thanks to Jianming Jin and the Department of Electrical and Computer Engineering at University of Illinois at Urbana-Champaign.