I'm having tremendous difficulty evaluating $_2F_1(a,b;c;z)$ with the
hypergeo package in R. In my case, values of $a$, $b$, $c$ are always positive real numbers. Even so, the hypergeometric function is incredibly sensitive to their values. I am not looking for extreme precision; I can use Excel to get a rough estimation of the Guass hypergeometric that is fine for my purposes.
Any suggestions for an implementation in R that will give a fast, error-free, if not super accurate Gaussian hypergeometric computation of positive real numbers with a wide range of values?
Edit: it seems there is far more code for this in MATLAB than R. Any thoughts as to why?