I have a task - to plot graphics of the function:

$$ I(E) = \frac{16i \pi k \mu}{(\beta - ik)^{4}} \frac{1}{1 + i\eta} {}_2 F_1(1+i\eta, 2; 2 + i \eta; x) $$


$$ x = \left( \frac{\beta + ik}{\beta - ik} \right)^2 $$ $$ k = \sqrt{\frac{\mu E}{20.9006}} $$ $$ \eta = 0.15748 \sqrt{\frac{\mu}{E}} $$ $$ \mu = \frac{m_d m_t}{m_d + m_t} \approx 1.2 $$

For task simplicity, $\beta = 1$. But in future I have to find it from the next equation:

$$ \left( \frac{\beta - ik}{\beta + ik} \right)^{2i\eta} = e^{4 \eta arctg{\frac{k}{\beta}}}, \beta \in \mathbb{R} $$

Energy parameter $E \in (0, 1]$ MeV, $i$ - is the imaginary unit.

Now my problem is the numerical calculation of the hypergeometric function ${}_2 F_1 (a, b; c; x)$. It is known that the hypergeometric function has a lot of representations. For example, there is a formula of Euler. But it requires that $Re (c) > Re (b) > 0$ and $|x| < 1$.

Question: What method or representation should I use to compute such type of hypergeometric function?

Service wolfram-alpha computes this function with no error...

  • $\begingroup$ Numerical recipes advice to solve the hypergeometric differential equation along a path (which does not intersect the branch cut at $x\in[1,\infty)$) starting from the circle $|x|=1/2$ where the series converges. $\endgroup$ – Fabian May 14 '12 at 12:17
  • $\begingroup$ Your numerator and denominator parameters are complex? Interesting. Where did this come from? $\endgroup$ – J. M. ain't a mathematician May 14 '12 at 12:19
  • $\begingroup$ @Fabian: to quote their words, "it's like swatting a fly with a golden brick". It works, but I'm sure there are more efficient ways of going about this. I'll have to do some digging... $\endgroup$ – J. M. ain't a mathematician May 14 '12 at 12:20
  • $\begingroup$ @J.M. for sure. But then one has to use different methods depending on $x$ (which I guess the OP wanted to avoid). $\endgroup$ – Fabian May 14 '12 at 12:22
  • $\begingroup$ @J.M. no, they are real(if you speak about $\eta, \beta, k, mu$). this is from physics - it is some kind of differential cross section for Helium. or similar. $\endgroup$ – stukselbax May 14 '12 at 14:01

Since $E > 0$, $x \not= 1$. Thus, the hypergeometric function ${}_2 F_1(1+i \eta, 2 ; 2+ i \eta ; x)$ is defined. Using Euler's integral representation with $a=2$, $b=1+ i \eta$ and $c = 2 + i \eta$, we get: $$ {}_2 F_1(1+i \eta, 2 ; 2+ i \eta ; x) = (1+ i \eta) \int_0^1 t^{i \eta} (1- t x)^{-2} \mathrm{d} t $$ The integral is convergent for real $\eta$ and complex $x$ not on $[1, \infty)$.

Apart from the integral representation, you may find the work by Joshua Willis on "Acceleration of generalized hypergeometric functions through precise remainder asymptotics" right on for your problem.

  • $\begingroup$ Thanks a lot. I use exactly this integral to compute hypergeometric function. I compute it using Simpson method, and all works for me. My problem was in computation of Gamma functions which have to stand before the integral - but you have reduce the fraction and simplified my task! $\endgroup$ – stukselbax May 15 '12 at 8:38

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.