As a general reference to the functions described in this section (and much more), we refer the reader to Whittaker and Watson [WW15].
Return the Gamma function Γ(f) of the series f. The series f must be defined over the free real or complex field, the valuation of f must be 0 and the constant term of f must be 1.
Given a real or complex number s (not equal to 0, - 1, - 2, ... ), calculate the value Γ(s) of the gamma function at s. For s with positive real part this is the value of the integral from 0 to infinity of ux - 1e - u with respect to u. For other s (not a non-positive integer) the function is defined by analytic continuation, and it satisfies the product formula (1/sΓ(s))=eγ s∏n=1^∞(1 + (s/n))e - s/n. The function Γ also satisfies Γ(s)Γ(1 - s)=(π/(sin(π s))), and Γ(s + 1)=sΓ(s).
Complementary: BoolElt Default: false
Gamma: FldReElt Default:
For real numbers s, t this returns the value of the incomplete gamma function γ(s, t)=int0t us - 1e - udu. The optional argument Complementary can be used to find the complement intt^∞us - 1e - udu instead. There is a second optional argument that may be used in the computation of the incomplete gamma value; the free real value of Gamma should be the value of Γ(s), in which case γ(s, t) may be computed as the difference between the given value for Γ(s) and that of the complementary γ at s, t. Pari is used here.
For free real s (such that s + (1/2) is not a non-positive integer) this returns the value of Γ(s + (1/2)). For integer values of s this is faster than Gamma(s+(1/2)), because Legendre's doubling formula Γ(s + (1/2))=21 - 2sSqrt(π)(Γ(2s)/Γ(s)) is used. Pari is used here.
Return the Log-Gamma function (Log)(Γ(f)) of the series f. f must be defined over a real or complex field, the valuation of f must be 0 and the constant term of f must be 1.
For real or complex s (not a non-positive integer) return the value of the principal branch of the logarithm of the gamma function of s.
For real or complex s (not a non-positive integer) return the principal value of the logarithmic derivative Ψ(s)=(d logΓ(s)/ds)=(Γ'(s)/Γ(s)), of the gamma function, which allows the expansion Ψ(s)= - γ - (1/s) + s∑n=1^∞(1/n(s + n)); here γ is Euler's gamma. Pari is used here.
Given a small integer n and a real number r, calculate the value of the Bessel function y = Jn(r), of the first kind of order n. Results for negative arguments are defined by: J - n(r) = Jn( - r) = ( - 1)n Jn(r). The Bessel function of the first kind of order n is defined by Jn(x)=(1/2πi)((z/2))nint - ∞0^ + u - n - 1e^(u - (z2/4t))du, and satisfies Jn(x)=∑k=0^∞(( - 1)kzn + 2k/2n + 2kk!Γ(n + k + 1)).
Given a small integer n and a real number r, calculate the value of the Bessel function y = Yn(r), of the second kind of order n. Results for negative arguments are defined by: Y - n(r) = - ( - 1)n Yn(r), Yn( - r) is not a real number. The Bessel function of the second kind of order n satisfies the Bessel differential equation.
Given a small integer n and a real number s, calculate the value of the Bessel function of the first kind of half integral index n + (1/2), J_(n + (1/2)), defined as above. Pari is used here.
Given a complex n and a positive real s, compute the value of the modified Bessel function of the second kind Kn(s), which may be defined by Kn(s)=(π/2)(inJ - n(i s) - i - nJn(s))cot(nπ). The function KBessel2 is an alternative (often faster) implementation of this function. Pari is used here.