SoatDev IT Consulting
SoatDev IT Consulting
  • About us
  • Expertise
  • Services
  • How it works
  • Contact Us
  • News
  • April 16, 2024
  • Rss Fetcher

It’s occasionally necessary to evaluate a hypergeometric function at a large negative argument. I was working on a project today that involved evaluating F(a, b; c; z) where z is a large negative number.

The hypergeometric function F(a, b; c; z) is defined by a power series in z whose coefficients are functions of a, b, and c. However, this power series has radius of convergence 1. This means you can’t use the series to evaluate F(a, b; c; −x) for x > 1.

It’s important to keep in mind the difference between a function and its power series representation. The former may exist where the latter does not. A simple example is the function f(x) = 1/(1 − x). The power series for this function has radius 1, but the function is defined everywhere except at x = 1.

Although the series defining F(a, b; c; z) is confined to the unit disk, the function itself is not. It can be extended analytically beyond the unit disk, usually with a branch cut along the real axis for z ≥ 1.

It’s good to know that our function can be evaluated for large negative x, but how do we evaluate it?

Linear transformation formulas

Hypergeometric functions satisfy a huge number of identities, the simplest of which are known as the linear transformation formulas even though they are not linear transformations of z. They involve bilinear transformations z, a.k.a. fractional linear transformations, a.k.a. Möbius transformations. [1]

One such transformation is the following, found in A&S 15.3.4 [2].

F(a, b; c; z) = (1-z)^{-a} Fleft(a, c-b; c; frac{z}{z-1} right)

If z < 1, then 0 < z/(z − 1) < 1, which is inside the radius of convergence. However, as z goes off to −∞, z/(z − 1) approaches 1, and the convergence of the power series will be slow.

A more complicated, but more efficient, formula is A&S 15.3.7, a linear transformation formula relates F at z to two other hypergeometric functions evaluated at 1/z. Now when z is large, 1/z is small, and these series will converge quickly.

begin{align*}
F(a, b; c; z) &= frac{Gamma(c) Gamma(b-a)}{Gamma(b) Gamma(c-a)} ,(-z)^{-aphantom{b}} Fleft(a, 1-c+a; 1-b+a; frac{1}{z}right) \
&+ frac{Gamma(c) Gamma(a-b)}{Gamma(a) Gamma(c-b)} ,(-z)^{-bphantom{a}} Fleft(,b, 1-c+b; 1-a+b; ,frac{1}{z}right) \
end{align*}

Related posts

  • Foreshadowing hypergoemetric functions
  • How to put a series in hypergeometric form
  • Branch cuts

[1] It turns out these transformations are linear, but not as functions of a complex argument. They’re linear as transformations on a projective space. More on that here.

[2] A&S refers to the venerable Handbook of Mathematical Functions by Abramowitz and Stegun.

The post Hypergeometric function of a large negative argument first appeared on John D. Cook.

Previous Post
Next Post

Recent Posts

  • Old Mutual Warns of Rising Deepfake Scams in South Africa
  • Cut Costs or Improve Customer Support? GenAI Says You Can Both
  • SANTACO Western Cape Signs Tech Partnership with Quickloc8
  • Deel wants Rippling to hand over any agreements involving paying the alleged spy
  • Telecommunications & IT Companies Form the Next Frontier for Growth & Cyber Crime across Africa

Categories

  • Industry News
  • Programming
  • RSS Fetched Articles
  • Uncategorized

Archives

  • May 2025
  • April 2025
  • February 2025
  • January 2025
  • December 2024
  • November 2024
  • October 2024
  • September 2024
  • August 2024
  • July 2024
  • June 2024
  • May 2024
  • April 2024
  • March 2024
  • February 2024
  • January 2024
  • December 2023
  • November 2023
  • October 2023
  • September 2023
  • August 2023
  • July 2023
  • June 2023
  • May 2023
  • April 2023

Tap into the power of Microservices, MVC Architecture, Cloud, Containers, UML, and Scrum methodologies to bolster your project planning, execution, and application development processes.

Solutions

  • IT Consultation
  • Agile Transformation
  • Software Development
  • DevOps & CI/CD

Regions Covered

  • Montreal
  • New York
  • Paris
  • Mauritius
  • Abidjan
  • Dakar

Subscribe to Newsletter

Join our monthly newsletter subscribers to get the latest news and insights.

© Copyright 2023. All Rights Reserved by Soatdev IT Consulting Inc.