In the previous post, I alluded to using Hermite polynomials in conjunction with higher-order Laplace approximation. In this post I’ll expand on what that means.

Hermite polynomials are orthogonal polynomials over the real line with respect to the weight given by the standard normal distribution. (There are two conventions for defining Hermite polynomials, what Wikipedia calls the physicist convention and the probabilist convention. We’re using the latter here.)

The first few Hermite polynomials are 1, *x*, *x*^{2} – 1, *x*^{3} – 3*x*, … . You can find the rest of the Hermite polynomials via the recurrence relation

*H*_{n+1}(*x*) = *x* *H*_{n} – *n* *H*_{n-1}(*x*).

You could think of the Hermite polynomials as **the right basis to use** when working with normal probability distributions. Writing a polynomial as a linear combination of Hermite polyn0mials is a change of basis that makes integration easy:

Here [*k* even] is the function that returns 1 if *k* is even and 0 otherwise. This notation was introduced by Kenneth Iverson in the APL programming language and has become moderately common in mathematics.