Appell series

From testwiki
Jump to navigation Jump to search

Template:Short description Template:Dablink In mathematics, Appell series are a set of four hypergeometric series F1, F2, F3, F4 of two variables that were introduced by Template:Harvs and that generalize Gauss's hypergeometric series 2F1 of one variable. Appell established the set of partial differential equations of which these functions are solutions, and found various reduction formulas and expressions of these series in terms of hypergeometric series of one variable.

Definitions

The Appell series F1 is defined for |x| < 1, |y| < 1 by the double series

F1(a,b1,b2;c;x,y)=∑m,n=0∞(a)m+n(b1)m(b2)n(c)m+nm!n!xmyn,

where (q)n is the rising factorial Pochhammer symbol. For other values of x and y the function F1 can be defined by analytic continuation. It can be shown[1] that

F1(a,b1,b2;c;x,y)=∑r=0∞(a)r(b1)r(b2)r(c−a)r(c+r−1)r(c)2rr!xryr2F1(a+r,b1+r;c+2r;x)2F1(a+r,b2+r;c+2r;y).

Similarly, the function F2 is defined for |x| + |y| < 1 by the series

F2(a,b1,b2;c1,c2;x,y)=∑m,n=0∞(a)m+n(b1)m(b2)n(c1)m(c2)nm!n!xmyn

and it can be shown[2] that

F2(a,b1,b2;c1,c2;x,y)=∑r=0∞(a)r(b1)r(b2)r(c1)r(c2)rr!xryr2F1(a+r,b1+r;c1+r;x)2F1(a+r,b2+r;c2+r;y).

Also the function F3 for |x| < 1, |y| < 1 can be defined by the series

F3(a1,a2,b1,b2;c;x,y)=∑m,n=0∞(a1)m(a2)n(b1)m(b2)n(c)m+nm!n!xmyn,

and the function F4 for |x|Template:Frac + |y|Template:Frac < 1 by the series

F4(a,b;c1,c2;x,y)=∑m,n=0∞(a)m+n(b)m+n(c1)m(c2)nm!n!xmyn.

The four series can each be represented as a simple series where the elements are Gaussian hypergeometric functions in terms of y:[3]

F1(a,b1,b2;c;x,y)=∑r=0∞(a)r(b1)r(c)rr!xr2F1(a+r,b2,c+r,y),
F2(a,b1,b2;c1,c2;x,y)=∑r=0∞(a)r(b1)r(c1)rr!xr2F1(a+r,b2,c2,y),
F3(a1,a2,b1,b2;c;x,y)=∑r=0∞(a1)r(b1)r(c)rr!xr2F1(a2,b2,c+r,y),
F4(a,b;c1,c2;x,y)=∑r=0∞(a)r(b)r(c1)rr!xr2F1(a+r,b+r,c2,y).

Similar expressions can be obtained upon exchange of the variables x and y and their respective parameters, e.g., c1 and c2 for F4.

Recurrence relations

Like the Gauss hypergeometric series 2F1, the Appell double series entail recurrence relations among contiguous functions. For example, a basic set of such relations for Appell's F1 is given by:

(a−b1−b2)F1(a,b1,b2,c;x,y)−aF1(a+1,b1,b2,c;x,y)+b1F1(a,b1+1,b2,c;x,y)+b2F1(a,b1,b2+1,c;x,y)=0,
cF1(a,b1,b2,c;x,y)−(c−a)F1(a,b1,b2,c+1;x,y)−aF1(a+1,b1,b2,c+1;x,y)=0,
cF1(a,b1,b2,c;x,y)+c(x−1)F1(a,b1+1,b2,c;x,y)−(c−a)xF1(a,b1+1,b2,c+1;x,y)=0,
cF1(a,b1,b2,c;x,y)+c(y−1)F1(a,b1,b2+1,c;x,y)−(c−a)yF1(a,b1,b2+1,c+1;x,y)=0.

Any other relation[4] valid for F1 can be derived from these four.

Similarly, all recurrence relations for Appell's F3 follow from this set of five:

cF3(a1,a2,b1,b2,c;x,y)+(a1+a2−c)F3(a1,a2,b1,b2,c+1;x,y)−a1F3(a1+1,a2,b1,b2,c+1;x,y)−a2F3(a1,a2+1,b1,b2,c+1;x,y)=0,
cF3(a1,a2,b1,b2,c;x,y)−cF3(a1+1,a2,b1,b2,c;x,y)+b1xF3(a1+1,a2,b1+1,b2,c+1;x,y)=0,
cF3(a1,a2,b1,b2,c;x,y)−cF3(a1,a2+1,b1,b2,c;x,y)+b2yF3(a1,a2+1,b1,b2+1,c+1;x,y)=0,
cF3(a1,a2,b1,b2,c;x,y)−cF3(a1,a2,b1+1,b2,c;x,y)+a1xF3(a1+1,a2,b1+1,b2,c+1;x,y)=0,
cF3(a1,a2,b1,b2,c;x,y)−cF3(a1,a2,b1,b2+1,c;x,y)+a2yF3(a1,a2+1,b1,b2+1,c+1;x,y)=0.

Derivatives and differential equations

For Appell's F1, the following derivatives result from the definition by a double series:

∂n∂xnF1(a,b1,b2,c;x,y)=(a)n(b1)n(c)nF1(a+n,b1+n,b2,c+n;x,y)
∂n∂ynF1(a,b1,b2,c;x,y)=(a)n(b2)n(c)nF1(a+n,b1,b2+n,c+n;x,y)

From its definition, Appell's F1 is further found to satisfy the following system of second-order differential equations:

x(1−x)∂2F1(x,y)∂x2+y(1−x)∂2F1(x,y)∂x∂y+[c−(a+b1+1)x]∂F1(x,y)∂x−b1y∂F1(x,y)∂y−ab1F1(x,y)=0
y(1−y)∂2F1(x,y)∂y2+x(1−y)∂2F1(x,y)∂x∂y+[c−(a+b2+1)y]∂F1(x,y)∂y−b2x∂F1(x,y)∂x−ab2F1(x,y)=0

A system partial differential equations for F2 is

x(1−x)∂2F2(x,y)∂x2−xy∂2F2(x,y)∂x∂y+[c1−(a+b1+1)x]∂F2(x,y)∂x−b1y∂F2(x,y)∂y−ab1F2(x,y)=0
y(1−y)∂2F2(x,y)∂y2−xy∂2F2(x,y)∂x∂y+[c2−(a+b2+1)y]∂F2(x,y)∂y−b2x∂F2(x,y)∂x−ab2F2(x,y)=0

The system have solution

F2(x,y)=C1F2(a,b1,b2,c1,c2;x,y)+C2x1−c1F2(a−c1+1,b1−c1+1,b2,2−c1,c2;x,y)+C3y1−c2F2(a−c2+1,b1,b2−c2+1,c1,2−c2;x,y)+C4x1−c1y1−c2F2(a−c1−c2+2,b1−c1+1,b2−c2+1,2−c1,2−c2;x,y)

Similarly, for F3 the following derivatives result from the definition:

∂∂xF3(a1,a2,b1,b2,c;x,y)=a1b1cF3(a1+1,a2,b1+1,b2,c+1;x,y)
∂∂yF3(a1,a2,b1,b2,c;x,y)=a2b2cF3(a1,a2+1,b1,b2+1,c+1;x,y)

And for F3 the following system of differential equations is obtained:

x(1−x)∂2F3(x,y)∂x2+y∂2F3(x,y)∂x∂y+[c−(a1+b1+1)x]∂F3(x,y)∂x−a1b1F3(x,y)=0
y(1−y)∂2F3(x,y)∂y2+x∂2F3(x,y)∂x∂y+[c−(a2+b2+1)y]∂F3(x,y)∂y−a2b2F3(x,y)=0

A system partial differential equations for F4 is

x(1−x)∂2F4(x,y)∂x2−y2∂2F4(x,y)∂y2−2xy∂2F4(x,y)∂x∂y+[c1−(a+b+1)x]∂F4(x,y)∂x−(a+b+1)y∂F4(x,y)∂y−abF4(x,y)=0
y(1−y)∂2F4(x,y)∂y2−x2∂2F4(x,y)∂x2−2xy∂2F4(x,y)∂x∂y+[c2−(a+b+1)y]∂F4(x,y)∂y−(a+b+1)x∂F4(x,y)∂x−abF4(x,y)=0

The system has solution

F4(x,y)=C1F4(a,b,c1,c2;x,y)+C2x1−c1F4(a−c1+1,b−c1+1,2−c1,c2;x,y)+C3y1−c2F4(a−c2+1,b−c2+1,c1,2−c2;x,y)+C4x1−c1y1−c2F4(2+a−c1−c2,2+b−c1−c2,2−c1,2−c2;x,y)

Integral representations

The four functions defined by Appell's double series can be represented in terms of double integrals involving elementary functions only Template:Harv. However, Template:Harvs discovered that Appell's F1 can also be written as a one-dimensional Euler-type integral:

F1(a,b1,b2,c;x,y)=Γ(c)Γ(a)Γ(c−a)∫01ta−1(1−t)c−a−1(1−xt)−b1(1−yt)−b2dt,ℜc>ℜa>0.

This representation can be verified by means of Taylor expansion of the integrand, followed by termwise integration.

Interestingly, from p. 4 of Y. Brychkov and N. Saad's 2015 paper, "On some formulas for the Appell function F3(a,a',b,b';c;w,z)", Appell's F3 function has the following one-dimensional integral representation in terms of two Gauss hypergeometric functions:

F3(a,a′,b,b′;c;w,z)=Γ(c)Γ(c′)Γ(c−c′)∫01tc′−1(1−t)c−c′−12F1(a,b;c′;wt)2F1(a′,b′;c−c′;z(1−t))dt,

where

(a,a′,b,b′,c,c′,w,z)∈ℂ∧c∉ℤ(−,0)∧|w|<1∧|z|<1∧ℜc>ℜc′>0.

Special cases

Picard's integral representation implies that the incomplete elliptic integrals F and E as well as the complete elliptic integral Π are special cases of Appell's F1:

F(ϕ,k)=∫0ϕdθ1−k2sin2θ=sin⁡(ϕ)F1(12,12,12,32;sin2ϕ,k2sin2ϕ),|ℜϕ|<π2,
E(ϕ,k)=∫0ϕ1−k2sin2θdθ=sin⁡(ϕ)F1(12,12,−12,32;sin2ϕ,k2sin2ϕ),|ℜϕ|<π2,
Π(n,k)=∫0π/2dθ(1−nsin2θ)1−k2sin2θ=π2F1(12,1,12,1;n,k2).

Furthermore, the above stated integral representation of Appell's F3 function gives us:

∫0xt(a2±t2)βarcsin⁡(bt)dt=13bx3(a2±x2)βF3(−β,12,1,12;52;x2x2±a2,b2x2)
∫0xt(a2±t2)β1−b2t2arcsin⁡(bt)dt=13bx3(a2±x2)βF3(−β,1,1,1;52;a2x2±a2,b2x2)
∫0aK(bt)K(ca−t)dt=π2a4F3(12,12,12,12;2;ac2,ab2)
∫0aE(ct)K(ba−t)dt=π2a4F3(12,−12,12,12;2;ab2,ac2)
∫0aE(ct)K(ba−t)1−c2tdt=π2a4F3(12,12,12,32;2;ab2,ac2)
∫0aE(bt)E(ca−t)dt=π2a4F3(−12,−12,12,12;2;ac2,ab2)
∫0aE(bt)E(ca−t)1−b2tdt=π2a4F3(−12,12,12,32;2;ac2,ab2)
∫0aE(bt)E(ca−t)(1−b2t)(1−c2(a−t))dt=π2a4F3(12,12,32,32;2;ac2,ab2)

where K(k) and E(k) represent the complete elliptic integrals of the first and second kind, respectively.

Template:Main

References

Template:Reflist

  1. ↑ See Burchnall & Chaundy (1940), formula (30).
  2. ↑ See Burchnall & Chaundy (1940), formula (26) or Erdélyi (1953), formula 5.12(9).
  3. ↑ See Appell & Kampé de Fériet (1926), Chapter I, formulae 15–18, page 15
  4. ↑ For example, (y−x)F1(a,b1+1,b2+1,c,x,y)=yF1(a,b1,b2+1,c,x,y)−xF1(a,b1+1,b2,c,x,y)