Skip to article frontmatterSkip to article content

Hermite integral formula

1Complex Analysis

We first recall some results from complex variables.

The necessary and sufficient condition is given by

2Laurent expansion and residue

If f(z)f(z) is not analytic at a single point z0z_0, then we can use a Laurent series. Let f(z)f(z) be analytic in the region

D={z:0<zz0<ρ}D = \{ z : 0 < |z-z_0| < \rho \}

and let z0z_0 be an isolated singular point of f(z)f(z). The Laurent expansion of f(z)f(z) in DD is given by

f(z)=n=Cn(zz0)nf(z) = \sum_{n=-\infty}^\infty C_n (z-z_0)^n

with

Cn=12πiCf(z)dz(zz0)n+1C_n = \frac{1}{2\pi \ii} \oint_C \frac{f(z) \ud z}{(z-z_0)^{n+1}}

where CC is a simple closed contour lying in DD. The coefficient C1C_{-1} is called the residue of f(z)f(z) at z0z_0. When n=1n=-1, we have

Cf(z)dz=2πiC1\oint_C f(z) \ud z = 2\pi \ii C_{-1}

3Lagrange interpolation

Let x0,x1,,xnx_0,x_1,\ldots,x_n be a set of n+1n+1 distinct interpolation points. Define

(x)=(xx0)(xxn)\ell(x) = (x-x_0)\ldots(x-x_n)

The Lagrange polynomials are given by

j(x)=(x)(xj)(xxj)\ell_j(x) = \frac{\ell(x)}{\ell'(x_j)(x-x_j)}

and the interpolant of degree nn is

p(x)=j=0nf(xj)j(x)p(x) = \sum_{j=0}^n f(x_j) \ell_j(x)

Let xx be any point which is distinct from the interpolation points. Let Γj\Gamma_j be a contour in the complex plane that encloses xjx_j but none of the other points, nor the point xx. Using the Cauchy integral formula, we see that

12πiΓjdt(t)(xt)=12πiΓjϕ(t)txjdtwhereϕ(t)=1(tx0)(txj1)(txj+1)(txn)(xt)=ϕ(xj)=1(xjx0)(xjxj1)(xjxj+1)(xjxn)(xxj)=1(xj)(xxj)\begin{aligned} & \frac{1}{2\pi\ii} \oint_{\Gamma_j} \frac{\ud t}{\ell(t) (x-t)} \\ =& \frac{1}{2\pi\ii} \oint_{\Gamma_j} \frac{\phi(t)}{t-x_j} \ud t \\ & \textrm{where} \qquad \phi(t) = \frac{1}{(t-x_0)\ldots(t-x_{j-1})(t-x_{j+1})\ldots(t-x_n)(x-t)} \\ =& \phi(x_j) \\ =& \frac{1}{(x_j-x_0)\ldots(x_j-x_{j-1})(x_j-x_{j+1})\ldots(x_j-x_n)(x-x_j)} \\ =& \frac{1}{\ell'(x_j)(x-x_j)} \end{aligned}

Hence

j(x)=(x)(xj)(xxj)=12πiΓj(x)dt(t)(xt)\ell_j(x) = \frac{\ell(x)}{\ell'(x_j)(x-x_j)} = \frac{1}{2\pi\ii} \oint_{\Gamma_j} \frac{\ell(x) \ud t}{\ell(t) (x-t)}

Now let Γ\Gamma' be a curve which encloses all of the points {xj}\{x_j\} but not the point xx and let ff be analytic on and interior to Γ\Gamma'. Then we can combine the Γj\Gamma_j integrals to get an expression for the interpolant

p(x)=12πijΓjf(xj)(x)dt(t)(xt)=12πijΓjf(t)(x)dt(t)(xt)=12πiΓ(x)f(t)dt(t)(xt)\begin{aligned} p(x) =& \frac{1}{2\pi\ii} \sum_j \oint_{\Gamma_j} f(x_j) \frac{\ell(x) \ud t}{\ell(t) (x- t)} \\ =& \frac{1}{2\pi\ii} \sum_j \oint_{\Gamma_j} f(t) \frac{\ell(x) \ud t}{\ell(t) (x-t)} \\ =& \frac{1}{2\pi\ii} \oint_{\Gamma'} \frac{\ell(x) f(t) \ud t}{\ell(t) (x-t)} \end{aligned}

by using the same type of argument we used above.

Now suppose we enlarge the contour of integration to a new contour Γ that encloses xx as well as {xj}\{x_j\}, and we assume ff is analytic on and inside Γ. Then

12πiΓ(x)f(t)dt(t)(xt)=12πiΓ(x)f(t)dt(t)(xt)+12πiΓ(x)f(t)dt(t)(xt)=p(x)f(x)\begin{aligned} \frac{1}{2\pi\ii} \oint_{\Gamma} \frac{\ell(x) f(t) \ud t}{\ell(t) (x-t)} =& \frac{1}{2\pi\ii} \oint_{\Gamma'} \frac{\ell(x) f(t) \ud t}{\ell(t) (x-t)} + \frac{1} {2\pi\ii} \oint_{\Gamma''} \frac{\ell(x) f(t) \ud t}{\ell(t) (x-t)} \\ =& p(x) - f(x) \end{aligned}

where we used the Cauchy integral formula to get f(x)-f(x).

4Effect of point distribution {xj}\{x_j\}

On a fixed contour Γ, the quantities f(t)f(t) and txt-x in the error formula are independent of {xj}\{x_j\}. The error depends on the ratio

(x)(t)=j=0n(xxj)j=0n(txj)\frac{\ell(x)}{\ell(t)} = \frac{ \prod\limits_{j=0}^n (x-x_j)}{\prod\limits_{j=0}^n(t- x_j)}

If Γ is far away from {xj}\{x_j\}, then for each tΓt \in \Gamma,

(x)(t)(x)tn+1,t\left| \frac{\ell(x)}{\ell(t)} \right| \approx \frac{|\ell(x)|}{|t|^{n+1}}, \qquad |t| \to \infty

this ratio will shrink exponentially as nn \to \infty, and if this happens, we may conclude that p(x)p(x) converges exponentially to f(x)f(x) as nn \to \infty. The crucial condition is that it must be possible to continue ff analytically far out to Γ.

Footnotes
  1. In practice, rounding errors may spoil the convergence.