Goodwin–Staton integral
In mathematics, the Goodwin–Staton integral is the special function [1]
The restriction on the argument excludes the non-positive real axis, where the denominator can vanish along the path of integration. The integral is named after E. T. Goodwin and J. Staton, who published a numerical table of the corresponding real integral in 1948.[2]
Relation to other special functions
For positive real , the Goodwin–Staton integral is related to Dawson's integral and the exponential integral by[3]
Differential equations
Differentiating the defining integral and integrating by parts gives the first-order inhomogeneous differential equation
Differentiating this relation twice and eliminating the inhomogeneous term gives the third-order linear homogeneous equation
Properties on the positive real axis
For , differentiation under the integral sign gives
Hence
so is a completely monotonic function on . In particular, is positive, strictly decreasing, and convex there.
Series near the origin
D. S. Jones studied the generalized Goodwin–Staton integral and obtained convergent series representations.[4] Specializing Jones's equation (2.11) gives, for ,
where is the principal logarithm and is the Euler–Mascheroni constant.
Consequently, as within the cut plane,
In particular,
Asymptotic expansion
For large positive , the integral has the asymptotic expansion[4]
Thus its first terms are
The NIST Digital Library of Mathematical Functions refers to Olver for corresponding expansions with bounds on the remainder for both real and complex values of the argument.[5]
Generalized Goodwin–Staton integral
Jones considered the generalized integral[4]
where is real with and lies off the positive real axis. The ordinary Goodwin–Staton integral is the special case
The generalized integral is a Stieltjes transform and satisfies the recurrence
as well as
Jones derived convergent series, asymptotic expansions, and explicit error bounds, and also treated oscillatory and logarithmic variants.[4] Mamedov later gave another evaluation of generalized Goodwin–Staton integrals using the binomial expansion theorem.[6]
Numerical evaluation
For positive real arguments, numerical evaluation can be reduced to Dawson's integral and the exponential integral through the relation above.[3] For small and moderate arguments, Jones's convergent series can be used, while for large arguments the asymptotic expansion is more suitable.[4] The DLMF section on computation of the Goodwin–Staton integral cites the original table of Goodwin and Staton.[7]