I'm having tremendous difficulty evaluating 2f1(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 hypergeo 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? Thanks.
Edit: it seems there is far more code for this in MATLAB than R. Any thoughts as to why?