Ultra fast Gauss hypergeometric function

This is the implementation of the Gauss hypergeometric function using series representation and suitable transformation formulas.

This implementation uses the series representation of the Gauss hypergeometric function and uses transformations such that the argument of the series, z, is always in the range [0,0.5] and therefore it converges very fast. It is basically the implementation suggested by "https://www.math.ucla.edu/~mason/research/pearson_final.pdf" written by John Pearson, section 4.6 Transformation formulae, equations 4.16-4.20.

It is compatible with all versions of MATLAB (as long as it contains the cumprod() function which is basically version 2006a to this day), however, the test file uses the native hypergeom() function implemented by Mathworks for comparison, beware!