Just wrote an R function to calculate the generalized hypergeometric function \(pFq\) for real numbers. The function is evaluated based on the series expansion of Pochhammer symbol. The function also allows you to have vectorized input.
You can find the code at
The function depends on the Pochhammer function which is available at