Computation of Hypergeometric Functions

Total Page:16

File Type:pdf, Size:1020Kb

Computation of Hypergeometric Functions Computation of Hypergeometric Functions by John Pearson Worcester College Dissertation submitted in partial fulfilment of the requirements for the degree of Master of Science in Mathematical Modelling and Scientific Computing University of Oxford 4 September 2009 Abstract We seek accurate, fast and reliable computations of the confluent and Gauss hyper- geometric functions 1F1(a; b; z) and 2F1(a; b; c; z) for different parameter regimes within the complex plane for the parameters a and b for 1F1 and a, b and c for 2F1, as well as different regimes for the complex variable z in both cases. In order to achieve this, we implement a number of methods and algorithms using ideas such as Taylor and asymptotic series com- putations, quadrature, numerical solution of differential equations, recurrence relations, and others. These methods are used to evaluate 1F1 for all z 2 C and 2F1 for jzj < 1. For 2F1, we also apply transformation formulae to generate approximations for all z 2 C. We discuss the results of numerical experiments carried out on the most effective methods and use these results to determine the best methods for each parameter and variable regime investigated. We find that, for both the confluent and Gauss hypergeometric functions, there is no simple answer to the problem of their computation, and different methods are optimal for different parameter regimes. Our conclusions regarding the best methods for computation of these functions involve techniques from a wide range of areas in scientific computing, which are discussed in this dissertation. We have also prepared MATLAB code that takes account of these conclusions. Acknowledgements I would like to offer special thanks to my project supervisors at the University of Oxford, Mason Porter and Sheehan Olver, for their invaluable help and support during the last few months. I would also like to thank the Numerical Algorithms Group, in particular David Sayers and Mick Pont, for providing the topic for this project, for supplying me with a copy of the NAG Toolbox, and for their assistance throughout. Thanks go to the National Institute of Standards and Technology, and in particular to Dan Lozier and Frank Olver, for kindly providing me with an advance copy of the new book, `NIST Digital Library of Mathematical Functions', a preliminary version of which is now hosted at http://dlmf.nist.gov/, which has given me many ideas about special functions and their computation. I am also grateful to Andy Wathen for giving me his thoughts on this dissertation. Further, I would like to thank the Engineering and Physical Sciences Research Council and the Numerical Algorithms Group for the provision of funding which has enabled me to take this course. Finally, I would like to express my gratitude to those who have taught and supervised me this year and throughout my time at Oxford, to my family for their continued support, and to my friends for making this year a very enjoyable one. Contents 1 Introduction 1 2 Background on hypergeometric functions 2 2.1 Relevant properties of hypergeometric functions . 2 2.2 Motivation for the computation of hypergeometric functions . 4 3 Computation of the confluent hypergeometric function 1F1(a; b; z) 6 3.1 Properties of 1F1(a; b; z) ............................. 7 3.2 Taylor series . 9 3.3 Writing the confluent hypergeometric function as a single fraction . 12 3.4 Buchholz polynomials . 14 3.5 Asymptotic series . 17 3.6 Quadrature methods . 19 3.7 Solving the confluent hypergeometric differential equation . 21 3.8 Recurrence relations . 24 3.9 Summary and analysis of results . 27 4 Computation of the Gauss hypergeometric function 2F1(a; b; c; z) 29 4.1 Properties of 2F1 ................................. 29 4.2 Taylor series . 31 4.3 Writing the Gauss hypergeometric function as a single fraction . 34 4.4 Quadrature methods . 35 4.5 Solving the hypergeometric differential equation . 37 4.6 Transformation formulae . 41 4.7 Analytic continuation formulae for z near e±iπ=3 . 43 4.8 Recurrence relations . 45 4.9 Summary and analysis of results . 49 5 Conclusions, Discussion and Future Considerations 50 A Transformation formulae for 2F1(a; b; c; z) when b − a 2 Z or c − a − b 2 Z 53 B List of test cases used for 1F1(a; b; z) 58 C List of test cases used for 2F1(a; b; c; z) 60 D Methods of testing the robustness of code selected 62 E Numerical results for 1F1(a; b; z) 66 F Numerical results for 2F1(a; b; c; z) 72 G Other methods considered for evaluating 1F1(a; b; z) 76 G.1 Series in terms of beta random variables . 76 G.2 Expansion in terms of incomplete gamma functions . 77 G.3 Asymptotic expansion for large jbj and jzj ................... 78 G.4 Hyperasymptotic expansions . 80 G.5 Other quadrature methods . 82 G.6 Other differential equation methods . 87 G.7 Pad´eapproximants and rational approximation . 89 G.8 Other expansions for 1F1(a; b; z)......................... 92 G.9 Multiplication formula . 93 H Other methods considered for evaluating 2F1(a; b; c; z) 94 H.1 Other quadrature methods . 94 H.2 Other differential equation methods . 96 H.3 Pad´eapproximants and rational approximation . 97 H.4 Chebyshev expansion for 2F1(a; b; c; z)...................... 98 I Evaluation of 0F1(; a; z) and other special functions required for this project 99 J Some examples of hypergeometric functions from practical applications 104 K List of code written for this project 109 Bibliography 113 1 Introduction The computation of the hypergeometric function pFq, a special function encountered in a variety of applications, is frequently sought. However, aside from the most basic hyper- geometric functions, this is an extremely difficult task in practice. The reason for this is that the non-trivial structure of the series that defines the function creates many numerical issues such as cancellation and round-off error, which become especially significant for cer- tain ranges of the parameters and the variable. This results in many methods of numerical computation being ineffective for all but the simplest parameter and variable ranges. We focus in this dissertation on computing the two most commonly used hypergeomet- ric functions, the confluent hypergeometric function 1F1(a; b; z) and the Gauss hypergeo- metric function 2F1(a; b; c; z), which both suffer from these problems. The program which will be used to compute these functions will be MATLAB, which employs double precision arithmetic, so it is especially important for effective methods to be found to overcome the numerical issues involved. The goal of this project is to carry out a comprehensive survey of methods for computing 1F1 and 2F1, and to determine how to choose appropriate methods for different parameter and variable ranges, resulting in reliable and fast computation for as large a range of the parameters (a, b for 1F1 and a, b, c for 2F1) and variable z as possible. The algorithms used are required to be accurate, fast and robust for specific parameter and variable regions for which they have been selected. In this dissertation, we will test a large variety of approaches, such as a range of se- ries methods, ways of numerically solving the differential equations that the hypergeometric functions satisfy and quadrature methods, as well as employing recurrence relations to at- tempt to use computations of relatively simple cases to obtain computations for extreme parameter cases. We will also test asymptotic series for 1F1. For 2F1, we also explore the use of transformation formulae and expansions for special parameter cases. We aim to combine techniques that work for specific parameter regimes in order to compile a package of methods that is effective for as large a part of the complex plane for each parameter and variable as possible. 1 This dissertation will explain the more successful methods tested for each function, with details of other methods implemented or analysed discussed in Appendices G and H. Details will be given of the background to these methods, how they were implemented, and the results obtained from test cases taken from a wide range of examples in the complex plane. An explanation will be provided as to why the methods might be successful, and a conclusion reached as to which methods should be used in what part of the complex plane. The conclusions will be based partly on the review of the author(s) who proposed the method and partly on results and observations from our implementation of them. 2 Background on hypergeometric functions In this section, we will introduce properties of the generalized hypergeometric function that will be exploited in this project. The motivation for computing hypergeometric functions will be discussed, with details given of some of the practical applications of these functions (with further examples given in Appendix J). We will also discuss the numerical issues that make the problem of their computation a difficult one. 2.1 Relevant properties of hypergeometric functions The hypergeometric function pFq is defined in [37] as follows for a1; :::; ap; b1; :::; bq; z 2 C: 1 j X (a1)j:::(ap)j z F (a ; :::; a ; b ; :::; b ; z) = ; (2.1) p q 1 p 1 q (b ) :::(b ) j! j=0 1 j q j where, for some parameter µ, the Pochhammer symbol (µ)j is defined as (µ)0 = 1; (µ)j = µ(µ + 1):::(µ + j − 1); j = 1; 2; ::: : For the remainder of this dissertation, the values of aj; bj; z will be complex unless other- wise specified. However, the hypergeometric function is not defined if any bj; j = 1; :::; q are real and equal to a non-positive integer, and there are numerical issues in its computation if one or more values of bj are close to a non-positive integer.
Recommended publications
  • Tsinghua Lectures on Hypergeometric Functions (Unfinished and Comments Are Welcome)
    Tsinghua Lectures on Hypergeometric Functions (Unfinished and comments are welcome) Gert Heckman Radboud University Nijmegen [email protected] December 8, 2015 1 Contents Preface 2 1 Linear differential equations 6 1.1 Thelocalexistenceproblem . 6 1.2 Thefundamentalgroup . 11 1.3 The monodromy representation . 15 1.4 Regular singular points . 17 1.5 ThetheoremofFuchs .. .. .. .. .. .. 23 1.6 TheRiemann–Hilbertproblem . 26 1.7 Exercises ............................. 33 2 The Euler–Gauss hypergeometric function 39 2.1 The hypergeometric function of Euler–Gauss . 39 2.2 The monodromy according to Schwarz–Klein . 43 2.3 The Euler integral revisited . 49 2.4 Exercises ............................. 52 3 The Clausen–Thomae hypergeometric function 55 3.1 The hypergeometric function of Clausen–Thomae . 55 3.2 The monodromy according to Levelt . 62 3.3 The criterion of Beukers–Heckman . 66 3.4 Intermezzo on Coxeter groups . 71 3.5 Lorentzian Hypergeometric Groups . 75 3.6 Prime Number Theorem after Tchebycheff . 87 3.7 Exercises ............................. 93 2 Preface The Euler–Gauss hypergeometric function ∞ α(α + 1) (α + k 1)β(β + 1) (β + k 1) F (α, β, γ; z) = · · · − · · · − zk γ(γ + 1) (γ + k 1)k! Xk=0 · · · − was introduced by Euler in the 18th century, and was well studied in the 19th century among others by Gauss, Riemann, Schwarz and Klein. The numbers α, β, γ are called the parameters, and z is called the variable. On the one hand, for particular values of the parameters this function appears in various problems. For example α (1 z)− = F (α, 1, 1; z) − arcsin z = 2zF (1/2, 1, 3/2; z2 ) π K(z) = F (1/2, 1/2, 1; z2 ) 2 α(α + 1) (α + n) 1 z P (α,β)(z) = · · · F ( n, α + β + n + 1; α + 1 − ) n n! − | 2 with K(z) the Jacobi elliptic integral of the first kind given by 1 dx K(z)= , Z0 (1 x2)(1 z2x2) − − p (α,β) and Pn (z) the Jacobi polynomial of degree n, normalized by α + n P (α,β)(1) = .
    [Show full text]
  • Transformation Formulas for the Generalized Hypergeometric Function with Integral Parameter Differences
    CORE Metadata, citation and similar papers at core.ac.uk Provided by Abertay Research Portal Transformation formulas for the generalized hypergeometric function with integral parameter differences A. R. Miller Formerly Professor of Mathematics at George Washington University, 1616 18th Street NW, No. 210, Washington, DC 20009-2525, USA [email protected] and R. B. Paris Division of Complex Systems, University of Abertay Dundee, Dundee DD1 1HG, UK [email protected] Abstract Transformation formulas of Euler and Kummer-type are derived respectively for the generalized hypergeometric functions r+2Fr+1(x) and r+1Fr+1(x), where r pairs of numeratorial and denominatorial parameters differ by positive integers. Certain quadratic transformations for the former function, as well as a summation theorem when x = 1, are also considered. Mathematics Subject Classification: 33C15, 33C20 Keywords: Generalized hypergeometric function, Euler transformation, Kummer trans- formation, Quadratic transformations, Summation theorem, Zeros of entire functions 1. Introduction The generalized hypergeometric function pFq(x) may be defined for complex parameters and argument by the series 1 k a1; a2; : : : ; ap X (a1)k(a2)k ::: (ap)k x pFq x = : (1.1) b1; b2; : : : ; bq (b1)k(b2)k ::: (bq)k k! k=0 When q = p this series converges for jxj < 1, but when q = p − 1 convergence occurs when jxj < 1. However, when only one of the numeratorial parameters aj is a negative integer or zero, then the series always converges since it is simply a polynomial in x of degree −aj. In (1.1) the Pochhammer symbol or ascending factorial (a)k is defined by (a)0 = 1 and for k ≥ 1 by (a)k = a(a + 1) ::: (a + k − 1).
    [Show full text]
  • A Note on the 2F1 Hypergeometric Function
    A Note on the 2F1 Hypergeometric Function Armen Bagdasaryan Institution of the Russian Academy of Sciences, V.A. Trapeznikov Institute for Control Sciences 65 Profsoyuznaya, 117997, Moscow, Russia E-mail: [email protected] Abstract ∞ α The special case of the hypergeometric function F represents the binomial series (1 + x)α = xn 2 1 n=0 n that always converges when |x| < 1. Convergence of the series at the endpoints, x = ±1, dependsP on the values of α and needs to be checked in every concrete case. In this note, using new approach, we reprove the convergence of the hypergeometric series for |x| < 1 and obtain new result on its convergence at point x = −1 for every integer α = 0, that is we prove it for the function 2F1(α, β; β; x). The proof is within a new theoretical setting based on a new method for reorganizing the integers and on the original regular method for summation of divergent series. Keywords: Hypergeometric function, Binomial series, Convergence radius, Convergence at the endpoints 1. Introduction Almost all of the elementary functions of mathematics are either hypergeometric or ratios of hypergeometric functions. Moreover, many of the non-elementary functions that arise in mathematics and physics also have representations as hypergeometric series, that is as special cases of a series, generalized hypergeometric function with properly chosen values of parameters ∞ (a ) ...(a ) F (a , ..., a ; b , ..., b ; x)= 1 n p n xn p q 1 p 1 q (b ) ...(b ) n! n=0 1 n q n X Γ(w+n) where (w)n ≡ Γ(w) = w(w + 1)(w + 2)...(w + n − 1), (w)0 = 1 is a Pochhammer symbol and n!=1 · 2 · ..
    [Show full text]
  • Lauricella's Hypergeometric Function
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 7, 452-470 (1963) Lauricella’s Hypergeometric Function FD* B. C. CARLSON Institute for Atomic Research and D@artment of Physics, Iowa State University, Ames, Iowa Submitted by Richard Bellman I. INTRODUCTION In 1880 Appell defined four hypergeometric series in two variables, which were generalized to 71variables in a straightforward way by Lauricella in 1893 [I, 21. One of Lauricella’s series, which includes Appell’s function F, as the case 71= 2 (and, of course, Gauss’s hypergeometric function as the case n = l), is where (a, m) = r(u + m)/r(a). The FD f unction has special importance for applied mathematics and mathematical physics because elliptic integrals are hypergeometric functions of type FD [3]. In Section II of this paper we shall define a hypergeometric function R(u; b,, . ..) b,; zi, “., z,) which is the same as FD except for small but important modifications. Since R is homogeneous in the variables zi, ..., z,, it depends in a nontrivial way on only 71- 1 ratios of these variables; indeed, every R function with n variables is expressible in terms of an FD function with 71- 1 variables, and conversely. Although R and FD are therefore equivalent, R turns out to be more convenient for both theory and applica- tion. The choice of R was initially suggested by the observation that the Euler transformations of FD would be greatly simplified by introducing homo- geneous variables.
    [Show full text]
  • Numerical Evaluation of the Gauss Hypergeometric Function with the Hypergeo Package
    Numerical evaluation of the Gauss hypergeometric function with the hypergeo package Robin K. S. Hankin Auckland University of Technology Abstract This paper introduces the hypergeo package of R routines, for numerical calculation of hypergeometric functions. The package is focussed on efficient and accurate evaluation of the hypergeometric function over the whole of the complex plane within the constraints of fixed-precision arithmetic. The hypergeometric series is convergent only within the unit circle, so analytic continuation must be used to define the function outside the unit circle. This short document outlines the numerical and conceptual methods used in the package; and justifies the package philosophy, which is to maintain transparent and verifiable links between the software and AMS-55. The package is demonstrated in the context of game theory. Keywords: Hypergeometric functions, numerical evaluation, complex plane, R, residue theo- rem. 1. Introduction k The geometric series k∞=0 tk with tk = z may be characterized by its first term and the con- k 1 stant ratio of successiveP terms tk+1/tk = z, giving the familiar identity k∞=0 z = (1 − z)− . Observe that while the series has unit radius of convergence, the rightP hand side is defined over the whole complex plane except for z = 1 where it has a pole. Series of this type may be generalized to a hypergeometric series in which the ratio of successive terms is a rational function of k: t +1 P (k) k = tk Q(k) where P (k) and Q(k) are polynomials. If both numerator and denominator have been com- pletely factored we would write t +1 (k + a1)(k + a2) ··· (k + ap) k = z tk (k + b1)(k + b2) ··· (k + bq)(k + 1) (the final term in the denominator is due to historical reasons), and if we require t0 = 1 then we write ∞ k a1,a2,...,ap tkz = aFb ; z (1) b1,b2,...,bq Xk=0 2 Numerical evaluation of the Gauss hypergeometric function with the hypergeo package − z when defined.
    [Show full text]
  • Representations of Hypergeometric Functions for Arbitrary Parameter
    Representations of hypergeometric functions for arbitrary parameter values and their use D.B. Karpa∗ and J.L. L´opezb aFar Eastern Federal University, Vladivostok, Russia bDpto. de Ingenier´ıa Matem´atica e Inform´atica, Universidad P´ublica de Navarra and INAMAT, Navarra, Spain Abstract. Integral representations of hypergeometric functions proved to be a very useful tool for studying their properties. The purpose of this paper is twofold. First, we extend the known representations to arbitrary values of the parameters and show that the extended representations can be interpreted as examples of regularizations of integrals containing Meijer’s G function. Second, we give new applications of both, known and extended representations. These include: inverse factorial series expansion for the Gauss type function, new information about zeros of the Bessel and Kummer type functions, connection with radial positive definite functions and generalizations of Luke’s inequalities for the Kummer and Gauss type functions. Keywords: generalized hypergeometric function, Meijer’s G function, integral representation, radial positive definite function, inverse factorial series, Hadamard finite part MSC2010: 33C20, 33C60, 33F05, 42A82, 65D20 1 Introduction Throughout the paper we will use the standard definition of the generalized hypergeometric function pFq as the sum of the series ∞ a (a ) (a ) (a ) arXiv:1609.02340v2 [math.CA] 5 Oct 2016 F z = F (a; b; z)= 1 n 2 n · · · p n zn (1) p q b p q (b ) (b ) (b ) n! n=0 1 n 2 n q n X · · · if p q, z C. If p = q + 1 the above series only converges in the open unit disk and F (z) is ≤ ∈ p q defined as its analytic continuation for z C [1, ).
    [Show full text]
  • Hypergeometric Integral Transforms
    Hypergeometric integral transforms François Rouvière April 15, 2018 The following integral formula plays an essential role in Cormack’s 1981 study [1] of the line Radon transform in the plane: s sT (p=s) rT (p=r) dp 2 n n = , n ; 0 < r < s; (1) 2 2 2 2 Z r s p p r p 2 Z where Tn denotes thep Chebyshevp polynomial de…ned by Tn(cos ) = cos n.A few years later, Cormack [2] extended his method to the Radon transform on certain hypersurfaces in Rn, using a similar formula with Chebyshev polynomials replaced by Gegenbauer’s. All these polynomials may be viewed as hypergeometric functions. Here we shall prove in Section 2 a more general integral formula (Proposition 3) and an inversion formula for a general hypergeometric integral transform (Theorem 4). These results are applied to Chebyshev polynomials in Section 3. For the reader’s convenience we give self-contained proofs, only sketched in the appendices to Cormack’s papers [1] and [2]. No previous knowledge of hypergeometric functions is assumed; the basics are given in Section 1. 1 Hypergeometric functions Let a; b; c; z denote complex numbers. For k N, let 2 (a)k := a(a + 1) (a + k 1) if k 1 , (a)0 := 1: For instance the classical binomial series may be written as k a z (1 z) = (a)k , z < 1; (2) k! j j k N X2 a where (1 z) denotes the principal value. It is useful to note that (a)k = 1 x s 1 (a + k)=(a) if a = N, where (s) := 0 e x dx is Euler’s Gamma function.
    [Show full text]
  • Multivariable Hypergeometric Functions
    Multivariable Hypergeometric Functions Eric M. Opdam Abstract. The goal of this lecture is to present an overview of the modern de- velopments around the theme of multivariable hypergeometric functions. The classical Gauss hypergeometric function shows up in the context of differen- tial geometry, algebraic geometry, representation theory and mathematical physics. In all cases it is clear that the restriction to the one variable case is unnatural. Thus from each of these contexts it is desirable to generalize the classical Gauss function to a class of multivariable hypergeometric func- tions. The theories that have emerged in the past decades are based on such considerations. 1. The Classical Gauss Hypergeometric Function The various interpretations of Gauss’ hypergeometric function have challenged mathematicians to generalize this function. Multivariable versions of this func- tion have been proposed already in the 19th century by Appell, Lauricella, and Horn. Reflecting developments in geometry, representation theory and mathemat- ical physics, a renewed interest in multivariable hypergeometric functions took place from the 1980’s. Such generalizations have been initiated by Aomoto [1], Gelfand and Gelfand [14], and Heckman and Opdam [19], and these theories have been further developed by numerous authors in recent years. The best introduction to this story is a recollection of the role of the Gauss function itself. So let us start by reviewing some of the basic properties of this classical function. General references for this introductory section are [24, 12], and [39]. The Gauss hypergeometric series with parameters a, b, c ∈ C and c ∈ Z≤0 is the following power series in z: ∞ (a) (b) F (a, b, c; z):= n n zn .
    [Show full text]
  • Note on the Hurwitz–Lerch Zeta Function of Two Variables
    S S symmetry Article Note on the Hurwitz–Lerch Zeta Function of Two Variables Junesang Choi 1, Recep ¸Sahin 2, O˘guzYa˘gcı 2 and Dojin Kim 3,* 1 Department of Mathematics, Dongguk University, Gyeongju 38066, Korea; [email protected] 2 Department of Mathematics, Kırıkkale University, Kırıkkale 71450, Turkey; [email protected] (R.S.); [email protected] (O.Y.) 3 Department of Mathematics, Dongguk University, Seoul 04620, Korea * Correspondence: [email protected] Received: 3 August 2020; Accepted: 27 August 2020; Published: 28 August 2020 Abstract: A number of generalized Hurwitz–Lerch zeta functions have been presented and investigated. In this study, by choosing a known extended Hurwitz–Lerch zeta function of two variables, which has been very recently presented, in a systematic way, we propose to establish certain formulas and representations for this extended Hurwitz–Lerch zeta function such as integral representations, generating functions, derivative formulas and recurrence relations. We also point out that the results presented here can be reduced to yield corresponding results for several less generalized Hurwitz–Lerch zeta functions than the extended Hurwitz–Lerch zeta function considered here. For further investigation, among possibly various more generalized Hurwitz–Lerch zeta functions than the one considered here, two more generalized settings are provided. Keywords: beta function; gamma function; Pochhammer symbol; Hurwitz–Lerch zeta function; Hurwitz–Lerch zeta function of two variables; hypergeometric functions; confluent hypergeometric functions; Appell hypergeometric functions; Humbert hypergeometric functions of two variables; integral representations; generating functions; derivative formulas; recurrence relation 1. Introduction and Preliminaries The generalized (or Hurwitz) zeta function z(s, n) is defined by (see, e.g., [1], pp.
    [Show full text]
  • Ramanujan's Master Theorem for the Hypergeometric Fourier Transform
    Ramanujan’s Master theorem for the hypergeometric Fourier transform on root systems Gestur Olafsson, Angela Pasquale To cite this version: Gestur Olafsson, Angela Pasquale. Ramanujan’s Master theorem for the hypergeometric Fourier transform on root systems. Journal of Fourier Analysis and Applications, Springer Verlag, 2013, 19 (6), pp.1150-1183. 10.1007/s00041-013-9290-5. hal-01232021 HAL Id: hal-01232021 https://hal.archives-ouvertes.fr/hal-01232021 Submitted on 21 Nov 2015 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. RAMANUJAN’S MASTER THEOREM FOR THE HYPERGEOMETRIC FOURIER TRANSFORM ON ROOT SYSTEMS G. OLAFSSON´ AND A. PASQUALE Abstract. Ramanujan’s Master theorem states that, under suitable conditions, the Mellin transform of an alternating power series provides an interpolation formula for the coefficients of this series. Ramanujan applied this theorem to compute several definite integrals and power series, which explains why it is referred to as the “Master Theorem”. In this paper we prove an analogue of Ramanujan’s Master theorem for the hypergeometric Fourier trans- form on root systems. This theorem generalizes to arbitrary positive multiplicity functions the results previously proven by the same authors for the spherical Fourier transform on semisimple Riemannian symmetric spaces.
    [Show full text]
  • Some Incomplete Hypergeometric Functions and Incomplete Riemann-Liouville Fractional Integral Operators
    mathematics Article Some Incomplete Hypergeometric Functions and Incomplete Riemann-Liouville Fractional Integral Operators Mehmet Ali Özarslan 1 and Ceren Ustao˘glu 2,* 1 Department of Mathematics, Faculty of Arts and Sciences, Eastern Mediterranean University, Famagusta, TRNC, Mersin 10, Turkey; [email protected] 2 Department of Computer Engineering, Final International University,Toroslar Caddesi, No. 6, Çatalköy, Girne, TRNC, Mersin 10, Turkey * Correspondence: ceren.ustaoglu@final.edu.tr Received: 19 March 2019; Accepted: 22 May 2019; Published: 27 May 2019 Abstract: Very recently, the incomplete Pochhammer ratios were defined in terms of the incomplete beta function By(x, z). With the help of these incomplete Pochhammer ratios, we introduce new incomplete Gauss, confluent hypergeometric, and Appell’s functions and investigate several properties of them such as integral representations, derivative formulas, transformation formulas, and recurrence relations. Furthermore, incomplete Riemann-Liouville fractional integral operators are introduced. This definition helps us to obtain linear and bilinear generating relations for the new incomplete Gauss hypergeometric functions. Keywords: Gauss hypergeometric function; confluent hypergeometric function; Appell’s functions; incomplete fractional calculus; Riemann-Liouville fractional integral; generating functions 1. Introduction and Preliminaries In recent years, some extensions of the well-known special functions have been considered by several authors (see, for example, [2–10] ). The familiar incomplete gamma functions g(s, x) and G(s, x) are defined by: Z x s−1 −t g(s, x) := t e dt (Re(s) > 0; x = 0) 0 and Z ¥ s−1 −t G(s, x) := t e dt (x = 0; Re(s) > 0 when x = 0), x respectively.
    [Show full text]
  • Notes on Differential Equations and Hypergeometric Functions
    Notes on differential equations and hypergeometric functions (NOT for publication) F.Beukers 21st August 2009 Contents 1 Ordinary linear differential equations 3 1.1 Differential equations and systems of equations . 3 1.2 Local theory . 7 1.3 Fuchsian equations . 12 1.4 Riemann-Hilbert correspondence . 14 1.5 Fuchsian equations of order two . 17 2 Gauss hypergeometric functions 21 2.1 Definition, first properties . 21 2.2 Monodromy of the hypergeometric function . 24 2.3 Triangle groups . 31 2.4 Some loose ends . 34 3 Hypergeometric functions nFn−1 36 3.1 Definition, first properties . 36 3.2 Monodromy . 37 3.3 Hypergeometric groups . 39 3.4 Rigidity . 41 4 Explicit monodromy 45 4.1 Introduction . 45 4.2 Mellin-Barnes integrals . 46 2 Chapter 1 Ordinary linear differential equations 1.1 Differential equations and systems of equa- tions A differential field K is a field equipped with a derivation, that is, a map ∂ : K → K which has the following properties, For all a, b ∈ K we have ∂(a + b) = ∂a + ∂b. For all a, b ∈ K we have ∂(ab) = a∂b + b∂a. The subset C := {a ∈ K|∂a = 0} is a subfield of K and is called the field of constants. We shall assume that C is algebraically closed and has characteristic zero. We shall also assume that ∂ is non-trivial that is, there exist a ∈ K such that ∂a 6= 0. Standard examples which will be used in later chapters are C(z), C((z)), C((z))an. They are the field of rational functions, formal Laurent series at z = 0 and Lau- rent series which converge in a punctured disk 0 < |z| < ρ for some ρ > 0.
    [Show full text]