How would we calculate the hypergeometric functions. One way would be to try and copy the code available in the mpmath package (or this FORTRAN code), although this would probably be rather complicated. Another option would be to use mpmath to calculate values over a large enough range then include the values in the code to be used as a look-up table. This should be rather quick.