Гипергеометрическая функция

Методы, которые мы применяем к биномиальным коэффициентам, оказываются весьма эффективными в тех случаях, когда они работают. Однако важно помнить, что их успешность часто носит специфический характер — эти методы, как правило, не универсальны и требуют тщательного подбора под конкретную задачу. Занимаясь решением разнообразных задач, мы иногда сталкиваемся с тем, что подходы, которые казались перспективными, неожиданно приводят к тупиковым путям. Биномиальные коэффициенты ведут себя подобно хамелеонам: их внешний вид легко меняется в зависимости от контекста, что делает их одновременно гибкими и сложными для систематизации. Поэтому естественным становится вопрос: существует ли некий унифицирующий объект, который обобщил бы всё многообразие методов суммирования? К счастью, такой объект действительно существует — им являются так называемые гипергеометрические ряды.

Начало изучению гипергеометрических рядов было положено в конце 18 – начале 19 века. Первые систематические исследования появились в работах Леонарда Эйлера в 1760-х годах, где он рассмотрел ряды с рекуррентным соотношением для коэффициентов. Окончательную форму гипергеометрической функции и глубокий анализ её свойств дал Карл Фридрих Гаусс в знаменитой работе 1812 года, посвящённой сходимости и функциональным соотношениям. Позже, в середине XIX века, Бернард Риман развил теорию с помощью комплексного анализа, введя геометрический взгляд на гипергеометрическое дифференциальное уравнение.

И до сих пор гипергеометрические ряды остаются предметом активных исследований — они находят применение в комбинаторике, теории чисел, математической физике и даже в современных вычислительных алгоритмах.

Гипергеометрическая функция и ряд

Гипергеометрическая функция — функция от комплексной переменной zz с p+qp + q параметрами

F(a1, a2, , apb1, b2, , bq | z)   ⁣=def   ⁣n=0a1na2napnb1nb2nbqnznn!\pFq{a_1 ,~ a_2 ,~ \dotsc ,~ a_p}{b_1 ,~ b_2 ,~ \dotsc ,~ b_q}{z} \defeq \sum\limits_{n=0}^\oo \frac{a_1^{\overline{n}} \cdot a_2^{\overline{n}} \dotsm a_p^{\overline{n}}} {b_1^{\overline{n}} \cdot b_2^{\overline{n}} \dotsm b_q^{\overline{n}}} \cdot \frac{z^n}{n!}

Параметры a1,a2,,apa_1, a_2, \dotsc, a_p называются верхними параметрами, а параметры b1,b2,,bqb_1, b_2, \dotsc, b_q называются нижними параметрами. Для того, чтобы знаменатель не превращался в 00, каждый нижний параметр не может быть нулём или целым отрицательным числом.

Строчная запись выглядит так: F(a1,a2,,ap; b1,b2,,bqz)\F(a_1, a_2, \dotsc, a_p; ~ b_1, b_2, \dotsc, b_q \mid z). Обращаю ваше внимание на то, что верхние и нижние параметры разделяются точкой с запятой, а аргумент, как и в полной нотации, отделяется вертикальной линией.

Конкретное значение гипергеометрической функции в точке zz называется гипергеометрическим рядом.

Огромное количество функций являются частными случаями гипергеометрических.

Например, если параметры вообще отсутствуют, мы получаем экспоненту:

F( | z)=n=0znn!=ez\pFq{-}{-}{z} = \sum\limits_{n=0}^\oo \frac{z^n}{n!} = e^z

Можно получить обычный геометрический ряд, сделав в числителе множитель n!=1nn! = 1^{\overline{n}}:

F(1 | z)=n=01nznn!=n=0n!znn!=n=0zn=11z\pFq{1}{-}{z} = \sum\limits_{n=0}^\oo 1^{\overline{n}} \cdot \frac{z^n}{n!} = \sum\limits_{n=0}^\oo n! \cdot \frac{z^n}{n!} = \sum\limits_{n=0}^\oo z^n = \frac{1}{1-z}

А если поставить в числитель произвольный параметр aa, можно получить

F(a | z)=n=0anznn!=n=0(a+n1n)zn=1(1z)a\pFq{a}{-}{z} = \sum\limits_{n=0}^\oo a^{\overline{n}} \cdot \frac{z^n}{n!} = \sum\limits_{n=0}^\oo \binom{a+n-1}{n} \, z^n = \frac{1}{(1-z)^a}

Заменив aa на a-a и zz на z-z можно получить биномиальный ряд

F(a | z)=(1+z)a\pFq{-a}{-}{-z} = (1+z)^a

Целое отрицательное число, стоящее на месте верхнего параметра, делает из бесконечного ряда конечный, поскольку (a)k=0(-a)^{\overline{k}} = 0 при k>ak > a и aN0a \in \NN_0.

Далее я просто приведу частные случаи гипергеометрических функций, чтобы вы убедились в универсальности инструмента:

ln(1+x)=xF(1, 12 | x)li(x)=γ+lnlnx+lnxF(1, 12, 2 | lnx)arcsin(x)=xF(1/2, 1/23/2 | x2)Si(x)=xF(1/2, 13/2, 3/2, 2 | x24)sin(x)=xF(3/2 | x24)cos(x)=F(1/2 | x24)Jν(x)=(x/2)νΓ(ν+1)F(ν+1 | x24)Iν(x)=(x/2)νΓ(ν+1)F(ν+1 | x24)γ(a, x)=xaaF(aa+1 | x)erf(x)=2xπF(1/23/2 | x2)K(k)=π2F(1/2, 1/21 | k2)E(k)=π2F(1/2, 1/21 | k2)\align{ \ln(1+x) &= x \cdot \pFq{1 ,~ 1}{2}{-x} & \li(x) &= \gamma + \ln |\ln x| + \ln x \cdot \pFq{1 ,~ 1}{2 ,~ 2}{\ln x}\\[0.8em]\arcsin(x) &= x \cdot \pFq{1/2 ,~ 1/2}{3/2}{x^2} & \Si(x) &= x \cdot \pFq{1/2 ,~ 1}{3/2 ,~ 3/2 ,~ 2}{-\frac{x^2}{4}}\\[0.8em]\sin(x) &= x \cdot \pFq{-}{3/2}{-\frac{x^2}{4}} & \cos(x) &= \pFq{-}{1/2}{-\frac{x^2}{4}}\\[0.8em]J_\nu(x) &= \frac{(x/2)^\nu}{\Gamma(\nu+1)} \cdot \pFq{-}{\nu+1}{-\frac{x^2}{4}} & I_\nu(x) &= \frac{(x/2)^\nu}{\Gamma(\nu+1)} \cdot \pFq{-}{\nu+1}{\frac{x^2}{4}}\\[0.8em]\gamma(a ,~ x) &= \frac{x^a}{a} \cdot \pFq{a}{a+1}{-x} & \erf(x) &= \frac{2x}{\sqrt{\pi}} \cdot \pFq{1/2}{3/2}{-x^2}\\[0.8em]\K(k) &= \frac{\pi}{2} \cdot \pFq{1/2 ,~ 1/2}{1}{k^2} & \E(k) &= \frac{\pi}{2} \cdot \pFq{-1/2 ,~ 1/2}{1}{k^2} & }

здесь Jν(x)J_\nu(x) — функция Бесселя первого рода, Iν(x)I_\nu(x) — модифицированная функция Бесселя, γ(a,x)\gamma(a, x) — неполная гамма-функция, Si(x)\Si(x) — интегральный синус, K(k)\K(k) и E(k)\E(k) — эллиптические интегралы первого и второго рода.

Пока гипергеометрические ряды нам не дали ничего, кроме громких имён и красивых определений. Однако мы увидели, что множество весьма разных функций можно представить именно в виде гипергеометрических рядов. Это и станет основным объектом нашего внимания.

Какие ряды являются гипергеометрическими? На этот вопрос легко ответить, если рассмотреть отношение последовательных членов гипергеометрического ряда

(a1n+1a2n+1apn+1b1n+1b2n+1bqn+1zn+1(n+1)!)/(a1na2napnb1nb2nbqnznn!)=(n+a1)(n+a2)(n+ap)(n+b1)(n+b2)(n+bq)zn+1\left. \left( \frac{a_1^{\overline{n+1}} \cdot a_2^{\overline{n+1}} \dotsm a_p^{\overline{n+1}}} {b_1^{\overline{n+1}} \cdot b_2^{\overline{n+1}} \dotsm b_q^{\overline{n+1}}} \cdot \frac{z^{n+1}}{(n+1)!} \right) \middle/ \left( \frac{a_1^{\overline{n}} \cdot a_2^{\overline{n}} \dotsm a_p^{\overline{n}}} {b_1^{\overline{n}} \cdot b_2^{\overline{n}} \dotsm b_q^{\overline{n}}} \cdot \frac{z^n}{n!} \right) \right. = \frac{(n+a_1) \, (n+a_2) \dotsm (n+a_p)}{(n+b_1) \, (n+b_2) \dotsm (n+b_q)} \cdot \frac{z}{n+1}

У нас получилась рациональная функция по nn, умноженная на zz. То есть если отношение соседних членов функционального ряда по zz представляет собой рациональную функцию по номеру членов, умноженную на zz, то этот функциональный ряд является гипергеометрической функцией.

Аналогично, если отношение соседних членов обычного числового ряда представляет собой рациональную функцию по номеру членов, то ряд является гипергеометрическим рядом.

Сразу скажу о способе упрощения выражений. Поскольку верхние индексы идут в числитель, а нижние в знаменатель, одинаковые числа можно сокращать. Например,

F(1, 2, 33, 4 | z)=F(1, 2, 33, 4 | z)=F(1, 24 | z)\pFq{1 ,~ 2 ,~ 3}{3 ,~ 4}{z} = \pFq{1 ,~ 2 ,~ \cancel{3}}{\cancel{3} ,~ 4}{z} = \pFq{1 ,~ 2}{4}{z}

Давайте, чтобы не путаться с обозначениями и определениями, разберём несколько примеров.

Первый пример. Посмотрим на числовой ряд n=0tn\sum\limits_{n=0}^\oo t_n, где

tn+1tn=n2+7n+104n2+1\frac{t_{n+1}}{t_n} = \frac{n^2 + 7n + 10}{4n^2 + 1}

Это рациональная функция по номеру nn. Значит, разложив числитель и знаменатель на множители, мы сможем найти параметры гипергеометрического ряда, который выражает исходный числовой ряд.

tn+1tn=n2+7n+104n2+1=(n+2)(n+5)4(ni/2)(n+i/2)=(n+2)(n+5)(n+1)(ni/2)(n+i/2)1/4n+1\frac{t_{n+1}}{t_n} = \frac{n^2 + 7n + 10}{4n^2 + 1} = \frac{(n+2) \, (n+5)}{4 \, (n-i/2) \, (n+i/2)} = \frac{(n+2) \, (n+5) \, (n+1)}{(n-i/2) \, (n+i/2)} \cdot \frac{1/4}{n+1}

Тогда наш исходный ряд легко выражается через гипергеометрический ряд

n=0tn=F(1, 2, 5i/2, i/2 | 14)\sum\limits_{n=0}^\oo t_n = \pFq{1 ,~ 2 ,~ 5}{-i/2 ,~ i/2}{\frac{1}{4}}

Второй пример. Посмотрим на ряд n=0tn\sum\limits_{n=0}^\oo t_n, где

tn+1tn=n22nz2n+6z3n2+3n+2\frac{t_{n+1}}{t_n} = \frac{n^2 - 2nz - 2n + 6z - 3}{n^2 + 3n + 2}

Это рациональная функция по nn, но подходит ли она под наше определение? Если рассматривать zz как аргумент гипергеометрической функции, то нет, она под свойство не подходит. Здесь zz входит линейно в числитель, и его нельзя выделить. Значит, n=0tn\sum\limits_{n=0}^\oo t_n не является гипергеометрической функцией по zz. Тем не менее, этот ряд является гипергеометрическим рядом:

tn+1tn=n22nz2n+6z3n2+3n+2=(n2z+1)(n3)(n+2)(n+1)=(n2z+1)(n3)(n+2)1n+1\frac{t_{n+1}}{t_n} = \frac{n^2 - 2nz - 2n + 6z - 3}{n^2 + 3n + 2} = \frac{(n-2z+1) \, (n-3)}{(n+2) \, (n+1)} = \frac{(n-2z+1) \, (n-3)}{(n+2)} \cdot \frac{1}{n+1}

Теперь, выделив отдельно верхние и нижние параметры, получим, что

n=0tn=F(12z, 32 | 1)\sum\limits_{n=0}^\oo t_n = \pFq{1-2z ,~ -3}{2}{1}

Ещё раз подчеркну, что это не гипергеометрическая функция от переменной zz, а гипергеометрический ряд, то есть конкретное значение гипергеометрической функции в точке 11. Здесь переменная zz присутствует в записи одного из параметров.

Третий пример. Рассмотрим разложение тангенса в ряд Тейлора вокруг нуля:

tanz=n=1B2n(4)n(14n)(2n)!z2n1\tan z = \sum\limits_{n=1}^\oo \frac{B_{2n} \, (-4)^n \, (1 - 4^n)}{(2n)!} \cdot z^{2n-1}

Попробуем вычислить отношение соседних ненулевых коэффициентов. Обозначим unu_n — коэффициент при z2n1z^{2n-1}, то есть Тогда отношение соседних коэффициентов равно

un+1un=B2n+2B2n(4)n+1(14n+1)(4)n(14n)(2n)!(2n+2)!=B2n+2B2n4(4n+11)4n11(2n+1)(2n+2)\frac{u_{n+1}}{u_n} = \frac{B_{2n+2}}{B_{2n}} \cdot \frac{(-4)^{n+1}(1 - 4^{n+1})}{(-4)^n(1 - 4^n)} \cdot \frac{(2n)!}{(2n+2)!} = \frac{B_{2n+2}}{B_{2n}} \cdot \frac{4(4^{n+1} - 1)}{4^n - 1} \cdot \frac{1}{(2n+1)(2n+2)}

Асимптотически числа Бернулли растут как B2n4πn(n/πe)2n|B_{2n}| \sim 4 \sqrt{\pi n} \, (n / \pi e)^{2n}, и их отношение B2n+2/B2nB_{2n+2} / B_{2n} ведёт себя как n2n^2 при больших nn. Но тут даже без асимптотики ясно, что отношение un+1/unu_{n+1}/u_n не является рациональной функцией от nn, поскольку числа Бернулли не задаются рациональной рекурсией.

Следовательно, ряд для tanz\tan z не является гипергеометрическим рядом, поскольку отношение его коэффициентов не является рациональной функцией от номера. Это особенно примечательно, ведь tanz=sinz/cosz\tan z = \sin z / \cos z, а обе функции в правой части гипергеометрические. Однако, как мы убедились только что, отношение гипергеометрических функций может и не является гипергеометрической функцией.

Четвёртый пример. А что там с синусом? Там всё гораздо лучше тангенса.

sinz=n=1(1)n+1(2n1)!z2n1=n=0(1)n(2n+1)!z2n+1=zn=0(1)n(2n+1)!z2n\sin z = \sum\limits_{n=1}^\oo \frac{(-1)^{n+1}}{(2n-1)!} \cdot z^{2n-1} = \sum\limits_{n=0}^\oo \frac{(-1)^n}{(2n+1)!} \cdot z^{2n+1} = z \cdot \sum\limits_{n=0}^\oo \frac{(-1)^n}{(2n+1)!} \cdot z^{2n}

Обозначим tnt_nnn-ый член ряда. Тогда отношение соседних коэффициентов равно

tn+1tn=(1)n+1z2n+2(2n+3)!/(1)nz2n(2n+1)!=z2(2n+3)(2n+2)=1n+3/2z2/4n+1\frac{t_{n+1}}{t_n} = \left. \frac{(-1)^{n+1} \cdot z^{2n+2}}{(2n+3)!} \middle/ \frac{(-1)^n \cdot z^{2n}}{(2n+1)!} \right. = -\frac{z^2}{(2n+3) \, (2n+2)} = -\frac{1}{n+3/2} \cdot \frac{z^2/4}{n+1}

Тогда сам синус выражается как

sinz=zF(3/2 | z24)\sin z = z \cdot \pFq{-}{3/2}{-\frac{z^2}{4}}

Обратите внимание на то, что я вынес множитель zz. При z=0z=0 значение синуса равно 00, а значение любого гипергеометрического ряда в точке 00 равно 11. После вынесения zz в разложении синуса я как раз получил ряд, значение которого в точке 00 равно 11, а значит приведение этого ряда к гипергеометрическому будет корректно.

Давайте посмотрим, как выглядит правило суммирования по диагонали для биномиальных коэффициентов в гипергеометрической форме

k=0n(r+kk)=(r+n+1n)\sum\limits_{k=0}^n \binom{r+k}{k} = \binom{r+n+1}{n}

Для начала сделаем из конечной суммы бесконечную, заменив kk на nkn-k

k=0n(r+kk)=k=0(r+nknk)=k=0(r+nk)!r!(nk)!\sum\limits_{k=0}^n \binom{r+k}{k} = \sum\limits_{k=0}^\oo \binom{r+n-k}{n-k} = \sum\limits_{k=0}^\oo \frac{(r+n-k)!}{r! \, (n-k)!}

Несмотря на то, что сумма формально бесконечная, она на самом деле конечная: множитель 1/(nk)!1/(n-k)! равен 00 при всех k>nk > n. Подробно этот неестественный факт разобран в статье «Гамма-функция».

Посмотрим теперь на отношение соседних членов. Будьте внимательны: здесь индексной переменной является kk, а не nn.

(r+nk1)!r!(nk1)!/(r+nk)!r!(nk)!=nkr+nk=kn(knr)(k+1)1k+1\left. \frac{(r+n-k-1)!}{r! \, (n-k-1)!} \middle/ \frac{(r+n-k)!}{r! \, (n-k)!} \right. = \frac{n-k}{r+n-k} = \frac{k-n}{(k-n-r) \, (k+1)} \cdot \frac{1}{k+1}

При этом первый (свободный) член ряда равен (r+nn)\binom{r+n}{n}. Значит, исходная сумма превращается в гипергеометрический ряд

k=0n(r+kk)=(r+nn)F(1, nnr | 1)\sum\limits_{k=0}^n \binom{r+k}{k} = \binom{r+n}{n} \cdot \pFq{1 ,~ -n}{-n-r}{1}

Тогда, подставляя значение суммы из правила суммирования по диагонали, получаем

(r+nn)F(1, nnr | 1)=(r+n+1n)\binom{r+n}{n} \cdot \pFq{1 ,~ -n}{-n-r}{1} = \binom{r+n+1}{n}

И можно на поделить на биномиальный коэффициент, и получить чуть более приятную формулу

F(1, nnr | 1)=r+n+1r+1если (r+nn)0\pFq{1 ,~ -n}{-n-r}{1} = \frac{r+n+1}{r+1} \quad\text{если}~ \binom{r+n}{n} \neq 0

Давайте теперь получим явную формулу для одного частного случая гипергеометрической функции — функции с двумя верхними параметрами aa, bb и одним нижним параметром cc. Эта функция, на самом деле, исторически была первой ...

Посмотрим на свёртку Вандермонда

k=0n(rk)(snk)=(r+sn)для r,s   ⁣   ⁣C и nZ\sum\limits_{k=0}^n \binom{r}{k} \binom{s}{n-k} = \binom{r+s}{n} \quad\text{для}~ r, s \;\! \in \;\! \CC ~\text{и}~ n \in \ZZ

Сразу запишем это тождество в гипергеометрической форме. Отношение соседних членов равно

(rk+1)(snk1)/(rk)(snk)=(rk)(nk)(k+1)(sn+k+1)=(kr)(kn)k+sn+11k+1\left. \binom{r}{k+1} \binom{s}{n-k-1} \middle/ \binom{r}{k} \binom{s}{n-k} \right. = \frac{(r-k) \, (n-k)}{(k+1) \, (s-n+k+1)} = \frac{(k-r) \, (k-n)}{k+s-n+1} \cdot \frac{1}{k+1}

При этом первый (свободный) член равен (sn)\binom{s}{n}. Получаем, что

k=0n(rk)(snk)=(sn)F(r, nsn+1 | 1)\sum\limits_{k=0}^n \binom{r}{k} \binom{s}{n-k} = \binom{s}{n} \cdot \pFq{-r ,~ -n}{s-n+1}{1}

По тождеству свёртки Вандермонда левая часть равна (r+sn)\binom{r+s}{n}. Поэтому

(r+sn)=(sn)F(r, nsn+1 | 1)\binom{r+s}{n} = \binom{s}{n} \cdot \pFq{-r ,~ -n}{s-n+1}{1}

Теперь, разделив обе части на (sn)\binom{s}{n}, получим выражение

F(r, nsn+1 | 1)=(r+sn)/(sn)=(r+s)nsn\pFq{-r ,~ -n}{s-n+1}{1} = \left. \binom{r+s}{n} \middle/ \binom{s}{n} \right. = \frac{(r+s)^{\underline{n}}}{s^{\underline{n}}}

Если заменить параметры a=ra = -r, b=nb = -n и c=sn+1c = s-n+1, то можно получить более естественный вид гипергеометрической функции

F(a, bc | 1)=(cab1b)/(cb1b)=Γ(c)Γ(cab)Γ(ca)Γ(cb)при bZ и b0 и a,cC\pFq{a ,~ b}{c}{1} = \left. \binom{c-a-b-1}{-b} \middle/ \binom{c-b-1}{-b} \right. = \frac{\Gamma(c) \cdot \Gamma(c-a-b)}{\Gamma(c-a) \cdot \Gamma(c-b)} \quad\text{при}~ b \in \ZZ ~\text{и}~ b \le 0 ~\text{и}~ a, c \in \CC

Здесь предлагаю немного отвлечься и поговорить о сходимости и об обобщениях. Не знаю, как вам, но мне доставляет моральный дискомфорт полученная формула, ведь у нас там параметры aa и cc являются любыми комплексными числами, а параметр bb только отрицательным целым числом. Хочется какой-то полной симметрии. Давайте её получим.

Давайте узнаем, в каких случаях наш частный гипергеометрический ряд сходится. Итак, смотрим на гипергеометрический ряд

F(a, bc | z)=n=0anbnccznn!при z=1\pFq{a ,~ b}{c}{z} = \sum\limits_{n=0}^\oo \frac{a^{\overline{n}} \cdot b^{\overline{n}}}{c^{\overline{c}}} \cdot \frac{z^n}{n!} \quad\text{при}~ |z| = 1

Вычислим на асимптотику коэффициента

anbnccn!na+bc1Γ(a)Γ(b)при n\frac{a^{\overline{n}} \cdot b^{\overline{n}}}{c^{\overline{c}} \cdot n!} \sim \frac{n^{a+b-c-1}}{\Gamma(a) \cdot \Gamma(b)} \quad\text{при}~ n \to \oo

Ряд n=0nα\sum\limits_{n=0}^\oo n^\alpha сходится при Reα<1\Re \alpha < -1. А это значит, что наш гипергеометрический ряд сходится при

Re(a+bc1)<1    Re(cab)>0\Re (a+b-c-1) < -1 \implies \Re (c-a-b) > 0

Далее нам нужно будет использовать частный случай интегрального представления Эйлера, которое мы в общем виде изучим позднее. А сейчас я просто покажу, как нужно проводить рассуждения. Заведём переменную t[0,1]t \in [0, 1]. Для этой переменной можно разложить выражение (1t)a(1-t)^{-a} в биномиальный ряд

(1t)a=n=0ann!tn(1-t)^{-a} = \sum\limits_{n=0}^\oo \frac{a^{\overline{n}}}{n!} \cdot t^n

Этот ряд сходится равномерно при всех значениях tt, то есть то, что я дальше буду делать, законно. Теперь рассмотрим интеграл

01tb1(1t)cb1(1t)adt\int\limits_0^1 t^{b-1} \, (1-t)^{c-b-1} \, (1-t)^{-a} \, dt

Подставим вместо (1t)a(1-t)^{-a} наш ряд и поменяем знаки интегрирования и суммирования местами

01tb1(1t)cb1(1t)adt=01tb1(1t)cb1n=0ann!tndt=n=0ann!01tb1(1t)cb1tndt==n=0ann!01tn+b1(1t)cb1dt=n=0ann!B(n+b,cb)=n=0ann!Γ(n+b)Γ(cb)Γ(n+c)==Γ(cb)Γ(b)Γ(c)n=0anbncn1n!=Γ(cb)Γ(b)Γ(c)n=0anbncn1n!=Γ(cb)Γ(b)Γ(c)F(a, bc | 1)\align{ &\int\limits_0^1 t^{b-1} \, (1-t)^{c-b-1} \, (1-t)^{-a} \, dt = \int\limits_0^1 t^{b-1} \, (1-t)^{c-b-1} \, \sum\limits_{n=0}^\oo \frac{a^{\overline{n}}}{n!} \cdot t^n \cdot dt = \sum\limits_{n=0}^\oo \frac{a^{\overline{n}}}{n!} \int\limits_0^1 t^{b-1} \, (1-t)^{c-b-1} \, t^n \, dt = \\ &\qquad = \sum\limits_{n=0}^\oo \frac{a^{\overline{n}}}{n!} \cdot \int\limits_0^1 t^{n+b-1} \, (1-t)^{c-b-1} \, dt = \sum\limits_{n=0}^\oo \frac{a^{\overline{n}}}{n!} \cdot \Beta(n+b, c-b) = \sum\limits_{n=0}^\oo \frac{a^{\overline{n}}}{n!} \cdot \frac{\Gamma(n+b) \cdot \Gamma(c-b)}{\Gamma(n+c)} = \\ &\qquad = \Gamma(c-b) \cdot \frac{\Gamma(b)}{\Gamma(c)} \cdot \sum\limits_{n=0}^\oo \frac{a^{\overline{n}} \cdot b^{\overline{n}}}{c^{\overline{n}}} \cdot \frac{1}{n!} = \Gamma(c-b) \cdot \frac{\Gamma(b)}{\Gamma(c)} \cdot \sum\limits_{n=0}^\oo \frac{a^{\overline{n}} \cdot b^{\overline{n}}}{c^{\overline{n}}} \cdot \frac{1}{n!} = \Gamma(c-b) \cdot \frac{\Gamma(b)}{\Gamma(c)} \cdot \pFq{a ,~ b}{c}{1} }

С другой стороны, исходный интеграл можно сразу схлопнуть в бета-функцию

01tb1(1t)cb1(1t)adt=01tb1(1t)cab1dt=B(cab,b)=Γ(cab)Γ(b)Γ(ca)\int\limits_0^1 t^{b-1} \, (1-t)^{c-b-1} \, (1-t)^{-a} \, dt = \int\limits_0^1 t^{b-1} \, (1-t)^{c-a-b-1} \, dt = \Beta(c-a-b, b) = \frac{\Gamma(c-a-b) \cdot \Gamma(b)}{\Gamma(c-a)}

Приравняв два вида исходного интеграла, получаем выражение для гипергеометрической функции

Γ(cb)Γ(b)Γ(c)F(a, bc | 1)=Γ(cab)Γ(b)Γ(ca)    F(a, bc | 1)=Γ(c)Γ(cab)Γ(ca)Γ(cb)\Gamma(c-b) \cdot \frac{\Gamma(b)}{\Gamma(c)} \cdot \pFq{a ,~ b}{c}{1} = \frac{\Gamma(c-a-b) \cdot \Gamma(b)}{\Gamma(c-a)} \implies \pFq{a ,~ b}{c}{1} = \frac{\Gamma(c) \cdot \Gamma(c-a-b)}{\Gamma(c-a) \cdot \Gamma(c-b)}

У нас получилась то же самое, что получалось при рассмотрении свёртки Вандермонда, однако здесь больше нет требования целочисленности параметра bb. Единственное требование на параметры, которое тут есть — гипергеометрический ряд должен сходиться, то есть Rec>Rea+Reb\Re c > \Re a + \Re b.

Дифференциальные соотношения

При исследовании функций нам обычно интересны не только их конкретные значения, но и их рост, поведение при разных аргументах и так далее. На большинство интересующих нас вопросов могут ответить производные и дифференциальные уравнения, которые получаются из этих функций.

Давайте продифференцируем гипергеометрическую функцию

ddzF(a1, a2, , apb1, b2, , bq | z)=ddzn=0a1na2napnb1nb2nbqnznn!=n=1a1na2napnb1nb2nbqnzn1(n1)!==n+1=1a1n+1a2n+1apn+1b1n+1b2n+1bqn+1znn!=n=0a1(a1+1)na2(a2+1)nap(ap+1)nb1(b1+1)nb2(b2+1)nbq(bq+1)nznn!==a1a2apb1b2bqn=0(a1+1)n(a2+1)n(ap+1)n(b1+1)n(b2+1)n(bq+1)nznn!==a1a2apb1b2bqF(a1+1, a2+1, , ap+1b1+1, b2+1, , bq+1 | z)\align{ &\frac{d}{dz} \, \pFq{a_1 ,~ a_2 ,~ \dotsc ,~ a_p}{b_1 ,~ b_2 ,~ \dotsc ,~ b_q}{z} = \frac{d}{dz} \, \sum\limits_{n=0}^\oo \frac{a_1^{\overline{n}} \cdot a_2^{\overline{n}} \dotsm a_p^{\overline{n}}} {b_1^{\overline{n}} \cdot b_2^{\overline{n}} \dotsm b_q^{\overline{n}}} \cdot \frac{z^n}{n!} = \sum\limits_{n=1}^\oo \frac{a_1^{\overline{n}} \cdot a_2^{\overline{n}} \dotsm a_p^{\overline{n}}} {b_1^{\overline{n}} \cdot b_2^{\overline{n}} \dotsm b_q^{\overline{n}}} \cdot \frac{z^{n-1}}{(n-1)!} = \\ &\qquad = \sum\limits_{n+1=1}^\oo \frac{a_1^{\overline{n+1}} \cdot a_2^{\overline{n+1}} \dotsm a_p^{\overline{n+1}}} {b_1^{\overline{n+1}} \cdot b_2^{\overline{n+1}} \dotsm b_q^{\overline{n+1}}} \cdot \frac{z^n}{n!} = \sum\limits_{n=0}^\oo \frac{a_1 \, (a_1+1)^{\overline{n}} \cdot a_2 \, (a_2+1)^{\overline{n}} \dotsm a_p \, (a_p+1)^{\overline{n}}} {b_1 \, (b_1+1)^{\overline{n}} \cdot b_2 \, (b_2+1)^{\overline{n}} \dotsm b_q \, (b_q+1)^{\overline{n}}} \cdot \frac{z^n}{n!} = \\ &\qquad = \frac{a_1 \cdot a_2 \dotsm a_p}{b_1 \cdot b_2 \dotsm b_q} \cdot \sum\limits_{n=0}^\oo \frac{(a_1+1)^{\overline{n}} \cdot (a_2+1)^{\overline{n}} \dotsm (a_p+1)^{\overline{n}}} {(b_1+1)^{\overline{n}} \cdot (b_2+1)^{\overline{n}} \dotsm (b_q+1)^{\overline{n}}} \cdot \frac{z^n}{n!} = \\ &\qquad = \frac{a_1 \cdot a_2 \dotsm a_p}{b_1 \cdot b_2 \dotsm b_q} \cdot \pFq{a_1+1 ,~ a_2+1 ,~ \dotsc ,~ a_p+1}{b_1+1 ,~ b_2+1 ,~ \dotsc ,~ b_q+1}{z} }

Получается, что при дифференцировании параметры выносятся и увеличиваются на 11.