MATHEMATICS OF COMPUTATION Volume 79, Number 270, April 2010, Pages 871–915 S 0025-5718(09)02280-7 Article electronically published on September 24, 2009
ON THE NUMERICAL EVALUATION OF FREDHOLM DETERMINANTS
FOLKMAR BORNEMANN
Abstract. Some significant quantities in mathematics and physics are most naturally expressed as the Fredholm determinant of an integral operator, most notably many of the distribution functions in random matrix theory. Though their numerical values are of interest, there is no systematic numerical treat- ment of Fredholm determinants to be found in the literature. Instead, the few numerical evaluations that are available rely on eigenfunction expansions of the operator, if expressible in terms of special functions, or on alternative, nu- merically more straightforwardly accessible analytic expressions, e.g., in terms of Painlev´e transcendents, that have masterfully been derived in some cases. In this paper we close the gap in the literature by studying projection methods and, above all, a simple, easily implementable, general method for the numer- ical evaluation of Fredholm determinants that is derived from the classical Nystr¨om method for the solution of Fredholm equations of the second kind. Using Gauss–Legendre or Clenshaw–Curtis as the underlying quadrature rule, we prove that the approximation error essentially behaves like the quadrature error for the sections of the kernel. In particular, we get exponential con- vergence for analytic kernels, which are typical in random matrix theory. The application of the method to the distribution functions of the Gaussian unitary ensemble (GUE), in the bulk scaling limit and the edge scaling limit, is dis- cussed in detail. After extending the method to systems of integral operators, we evaluate the two-point correlation functions of the more recently studied Airy and Airy1 processes.
1. Introduction Fredholm’s (1903) landmark paper1 on linear integral equations is generally con- sidered to be the forefather of those mathematical concepts that finally led to mod- ern functional analysis and operator theory; see the historical accounts in Kline (1972, pp. 1058–1075) and Dieudonn´e (1981, Chap. V). Fredholm was interested in the solvability of what is now called a Fredholm equation of the second kind, b (1.1) u(x)+z K(x, y)u(y) dy = f(x)(x ∈ (a, b)), a and explicit formulas for the solution thereof, for a right hand side f and a ker- nel K, both assumed to be continuous functions. He introduced his now famous
Received by the editor June 24, 2008 and, in revised form, March 16, 2009. 2000 Mathematics Subject Classification. Primary 65R20, 65F40; Secondary 47G10, 15A52. 1Simon (2005, p. VII) writes: “This deep paper is extremely readable and I recommend it to those wishing a pleasurable afternoon.” An English translation of the paper can be found in Birkhoff (1973, pp. 449–465).