# Mathieu function

(Redirected from Mathieu differential equation)

In mathematics, Mathieu functions are solutions of Mathieu's differential equation

${\displaystyle {\frac {d^{2}y}{dx^{2}}}+(a-2q\cos 2x)y=0,}$

where ${\displaystyle a}$ and ${\displaystyle q}$ are parameters. They were first introduced by Émile Léonard Mathieu, who encountered them while studying vibrating elliptical drumheads.[1][2] They have applications in many fields of the physical sciences, such as optics, quantum mechanics, and general relativity. In general, they tend to occur in problems involving some sort of periodic motion or in the analysis of partial differential equation boundary value problems possessing elliptic symmetry.[3]

## Definition

### Mathieu functions

In some usages, Mathieu function refers to solutions of the Mathieu differential equation for arbitrary values of ${\displaystyle a}$  and ${\displaystyle q}$ . When no confusion can arise, other authors use Mathieu function to refer specifically to ${\displaystyle \pi }$ - or ${\displaystyle 2\pi }$ -periodic solutions, which exist only for special values of ${\displaystyle a}$  and ${\displaystyle q}$ .[4] More precisely, for given (real) ${\displaystyle q}$  such periodic solutions exist for an infinite number of values of ${\displaystyle a}$ , called characteristic numbers, conventionally indexed as two separate sequences ${\displaystyle a_{n}(q)}$  and ${\displaystyle b_{n}(q)}$ , for ${\displaystyle n=1,2,3,\ldots }$ . The corresponding functions are denoted ${\displaystyle {\text{ce}}_{n}(x,q)}$  and ${\displaystyle {\text{se}}_{n}(x,q)}$ , respectively. They are sometimes also referred to as cosine-elliptic and sine-elliptic, or Mathieu functions of the first kind.

As a result of assuming that ${\displaystyle q}$  is real, both the characteristic numbers and associated functions are real-valued.[5]

${\displaystyle {\text{ce}}_{n}(x,q)}$  and ${\displaystyle {\text{se}}_{n}(x,q)}$  can be further classified by parity and periodicity (both with respect to ${\displaystyle x}$ ), as follows:[4]

Function Parity Period
${\displaystyle {\text{ce}}_{n},\ n{\text{ even}}}$  even ${\displaystyle 2\pi }$
${\displaystyle {\text{ce}}_{n},\ n{\text{ odd}}}$  even ${\displaystyle \pi }$
${\displaystyle {\text{se}}_{n},\ n{\text{ even}}}$  odd ${\displaystyle 2\pi }$
${\displaystyle {\text{se}}_{n},\ n{\text{ odd}}}$  odd ${\displaystyle \pi }$

The indexing with the integer ${\displaystyle n}$ , besides serving to arrange the characteristic numbers in ascending order, is convenient in that ${\displaystyle {\text{ce}}_{n}(x,q)}$  and ${\displaystyle {\text{se}}_{n}(x,q)}$  become proportional to ${\displaystyle \cos nx}$  and ${\displaystyle \sin nx}$  as ${\displaystyle q\rightarrow 0}$ . With ${\displaystyle n}$  being an integer, this gives rise to the classification of ${\displaystyle {\text{ce}}_{n}}$  and ${\displaystyle {\text{se}}_{n}}$  as Mathieu functions (of the first kind) of integral order. For general ${\displaystyle a}$  and ${\displaystyle q}$ , solutions besides these can be defined, including Mathieu functions of fractional order as well as non-periodic solutions.

### Modified Mathieu functions

Closely related are the modified Mathieu functions, which are solutions of Mathieu's modified differential equation

${\displaystyle {\frac {d^{2}y}{dx^{2}}}-(a-2q\cosh 2x)y=0,}$

which can be related to the original Mathieu equation by taking ${\displaystyle x\to \pm xi}$ . Accordingly, the modified Mathieu functions of the first kind of integral order, denoted by ${\displaystyle {\text{Ce}}_{n}(x,q)}$  and ${\displaystyle {\text{Se}}_{n}(x,q)}$ , are defined from[6]

{\displaystyle {\begin{aligned}{\text{Ce}}_{n}(x,q)&={\text{ce}}_{n}(ix,q).\\{\text{Se}}_{n}(x,q)&=-i{\text{se}}_{n}(ix,q).\end{aligned}}}

These functions are real-valued when ${\displaystyle x}$  is real.

### Normalization

A common normalization,[7] which will be adopted throughout this article, is to demand

${\displaystyle \int _{0}^{2\pi }{\text{ce}}_{n}(x,q)^{2}dx=\int _{0}^{2\pi }{\text{se}}_{n}(x,q)^{2}dx=\pi }$

as well as require ${\displaystyle {\text{ce}}_{n}(x,q)\rightarrow +\cos nx}$  and ${\displaystyle {\text{se}}_{n}(x,q)\rightarrow +\sin nx}$  as ${\displaystyle q\rightarrow 0}$ .

## Floquet theory

Many properties of the Mathieu differential equation can be deduced from the general theory of ordinary differential equations with periodic coefficients, called Floquet theory. The central result is Floquet's theorem:

Floquet's theorem[8] Mathieu's equation always has at least one solution ${\displaystyle y(x)}$  such that ${\displaystyle y(x+\pi )=\sigma y(x)}$ , where ${\displaystyle \sigma }$  is a constant which depends on the parameters of the equation and may be real or complex.

It is natural to associate the characteristic numbers ${\displaystyle a(q)}$  with those values of ${\displaystyle a}$  which result in ${\displaystyle \sigma =\pm 1}$ .[9] Note, however, that the theorem only guarantees the existence of at least one solution satisfying ${\displaystyle y(x+\pi )=\sigma y(x)}$ , when Mathieu's equation in fact has two independent solutions for any given ${\displaystyle a}$ , ${\displaystyle q}$ . Indeed, it turns out that with ${\displaystyle a}$  equal to one of the characteristic numbers, Mathieu's equation has only one periodic solution (that is, with period ${\displaystyle \pi }$  or ${\displaystyle 2\pi }$ ), and this solution is one of the ${\displaystyle {\text{ce}}_{n}(x,q)}$ , ${\displaystyle {\text{se}}_{n}(x,q)}$ . The other solution is nonperiodic, denoted ${\displaystyle {\text{fe}}_{n}(x,q)}$  and ${\displaystyle {\text{ge}}_{n}(x,q)}$ , respectively, and referred to as a Mathieu function of the second kind.[10] This result can be formally stated as Ince's theorem:

Ince's theorem[11] Define a basically periodic function as one satisfying ${\displaystyle y(x+\pi )=\pm y(x)}$ . Then, except in the trivial case ${\displaystyle q=0}$ , Mathieu's equation never possesses two (independent) basically periodic solutions for the same values of ${\displaystyle a}$  and ${\displaystyle q}$ .

An example ${\displaystyle P(a,q,x)}$  from Floquet's theorem, with ${\displaystyle a=1}$ , ${\displaystyle q=1/5}$ , ${\displaystyle \mu \approx 1+0.0995i}$  (real part, red; imaginary part, green)

An equivalent statement of Floquet's theorem is that Mathieu's equation admits a complex-valued solution of form

${\displaystyle F(a,q,x)=\exp(i\mu \,x)\,P(a,q,x),}$

where ${\displaystyle \mu }$  is a complex number, the Floquet exponent (or sometimes Mathieu exponent), and ${\displaystyle P}$  is a complex valued function periodic in ${\displaystyle x}$  with period ${\displaystyle \pi }$ . An example ${\displaystyle P(a,q,x)}$  is plotted to the right.

## Other types of Mathieu functions

### Mathieu functions of the second kind

Since Mathieu's equation is a second order differential equation, it always possesses two independent solutions. Floquet's theory says that if ${\displaystyle a}$  is equal to a characteristic number, one of these solutions is periodic, and the other nonperiodic. The periodic solution is one of the ${\displaystyle {\text{ce}}_{n}(x,q)}$  and ${\displaystyle {\text{se}}_{n}(x,q)}$ , called a Mathieu function of the first kind of integral order. The nonperiodic one is denoted either ${\displaystyle {\text{fe}}_{n}(x,q)}$  and ${\displaystyle {\text{ge}}_{n}(x,q)}$ , respectively, and is called Mathieu function of the second kind of integral order. The nonperiodic solutions are unstable, that is, they diverge as ${\displaystyle z\rightarrow \pm \infty }$ .[12]

The second solutions corresponding to the modified Mathieu functions ${\displaystyle {\text{Ce}}_{n}(x,q)}$  and ${\displaystyle {\text{Se}}_{n}(x,q)}$  are naturally defined as ${\displaystyle {\text{Fe}}_{n}(x,q)=-i{\text{fe}}_{n}(xi,q)}$  and ${\displaystyle {\text{Ge}}_{n}(x,q)={\text{ge}}_{n}(xi,q)}$ .

### Mathieu functions of fractional order

Mathieu functions of fractional order can be defined as those solutions ${\displaystyle {\text{ce}}_{p}(x,q)}$  and ${\displaystyle {\text{se}}_{p}(x,q)}$ , ${\displaystyle p}$  a non-integer, which turn into ${\displaystyle \cos px}$  and ${\displaystyle \sin px}$  as ${\displaystyle q\rightarrow 0}$ .[6] If ${\displaystyle p}$  is irrational, they are non-periodic; however, they remain bounded as ${\displaystyle x\rightarrow \infty }$ .

An important property of the solutions ${\displaystyle {\text{ce}}_{p}(x,q)}$  and ${\displaystyle {\text{se}}_{p}(x,q)}$ , for ${\displaystyle p}$  non-integer, is that they exist for the same value of ${\displaystyle a}$ . In contrast, when ${\displaystyle p}$  is an integer, ${\displaystyle {\text{ce}}_{p}(x,q)}$  and ${\displaystyle {\text{se}}_{p}(x,q)}$  never occur for the same value of ${\displaystyle a}$ . (See Ince's Theorem above.)

These classifications are summarized in the table below. The modified Mathieu function counterparts are defined similarly.

Classification of Mathieu functions[13]
Order First kind Second kind
Integral ${\displaystyle {\text{ce}}_{n}(x,q)}$  ${\displaystyle {\text{fe}}_{n}(x,q)}$
Integral ${\displaystyle {\text{se}}_{n}(x,q)}$  ${\displaystyle {\text{ge}}_{n}(x,q)}$
Fractional

(${\displaystyle p}$  non-integral)

${\displaystyle {\text{ce}}_{p}(x,q)}$  ${\displaystyle {\text{se}}_{p}(x,q)}$

## Explicit representation and computation

### First kind

Mathieu functions of the first kind can be represented as Fourier series:[4]

{\displaystyle {\begin{aligned}{\text{ce}}_{2n}(x,q)&=\sum _{r=0}^{\infty }A_{2r}^{(2n)}\cos 2rx\\{\text{ce}}_{2n+1}(x,q)&=\sum _{r=0}^{\infty }A_{2r+1}^{(2n+1)}\cos(2r+1)x\\{\text{se}}_{2n+1}(x,q)&=\sum _{r=0}^{\infty }B_{2r+1}^{(2n+1)}\sin(2r+1)x\\{\text{se}}_{2n+2}(x,q)&=\sum _{r=0}^{\infty }B_{2r+2}^{(2n+2)}\sin(2r+2)x\\\end{aligned}}}

By substitution into the Mathieu equation, the expansion coefficients ${\displaystyle A}$  and ${\displaystyle B}$  can be shown to obey three-term recurrence relations in the ${\displaystyle r}$  index. For instance, for each ${\displaystyle {\text{ce}}_{2n}}$  one finds[14]

{\displaystyle {\begin{aligned}aA_{0}-qA_{2}&=0\\(a-4)A_{2}-q(A_{4}+2A_{0})&=0\\(a-4r^{2})A_{2r}-q(A_{2r+2}+A_{2r-2})&=0,\quad r\geq 2\end{aligned}}}

where the ${\displaystyle n}$  superscripts have been omitted for clarity. Being a second-order recurrence in the index ${\displaystyle 2r}$ , one can always find two independent solutions ${\displaystyle X_{2r}}$  and ${\displaystyle Y_{2r}}$  such that the general solution can be expressed as a linear combination of the two: ${\displaystyle A_{2r}=c_{1}X_{2r}+c_{2}Y_{2r}}$ . Moreover, in this particular case, an asymptotic analysis[15] shows that one possible choice of fundamental solutions has the property

{\displaystyle {\begin{aligned}X_{2r}&=r^{-2r-1}\left(-{\frac {e^{2}q}{4}}\right)^{r}\left[1+{\mathcal {O}}(r^{-1})\right]\\Y_{2r}&=r^{2r-1}\left(-{\frac {4}{e^{2}q}}\right)^{r}\left[1+{\mathcal {O}}(r^{-1})\right]\end{aligned}}}

In particular, ${\displaystyle X_{2r}}$  is finite whereas ${\displaystyle Y_{2r}}$  diverges. Writing ${\displaystyle A_{2r}=c_{1}X_{2r}+c_{2}Y_{2r}}$ , we therefore see that in order for the Fourier series representation of ${\displaystyle {\text{ce}}_{2n}}$  to converge, ${\displaystyle a}$  must be chosen such that ${\displaystyle c_{2}=0}$ . These choices of ${\displaystyle a}$  correspond to the characteristic numbers.

In general, however, the solution of a three-term recurrence with variable coefficients cannot be represented in a simple manner, and hence there is no simple way to determine ${\displaystyle a}$  from the condition ${\displaystyle c_{2}=0}$ . Moreover, even if the approximate value of a characteristic number is known, it cannot be used to obtain the coefficients ${\displaystyle A_{2r}}$  by numerically iterating the recurrence towards increasing ${\displaystyle r}$ . The reason is that as long as ${\displaystyle a}$  only approximates a characteristic number, ${\displaystyle c_{2}}$  is not identically ${\displaystyle 0}$  and the divergent solution ${\displaystyle Y_{2r}}$  eventually dominates for large enough ${\displaystyle r}$ .

To overcome these issues, more sophisticated semi-analytical/numerical approaches are required, for instance using a continued fraction expansion,[16][4] casting the recurrence as a matrix eigenvalue problem,[17] or implementing a backwards recurrence algorithm.[15] The complexity of the three-term recurrence relation is one of the reasons there are few simple formulas and identities involving Mathieu functions.[18]

Mathieu functions and their associated characteristic numbers ${\displaystyle a}$  (in physics literature called eigenvalues) for small values of ${\displaystyle q}$  (and thus in rising powers of ${\displaystyle q}$ ) for small integral values of the Floquet exponent ${\displaystyle \mu }$  are most easily obtained perturbatively [19].

In practice, Mathieu functions and the corresponding characteristic numbers can be readily calculated using pre-packaged software, such as Mathematica, Maple, MATLAB, and SciPy.

### Second kind

Mathieu functions of the second kind can be expressed in terms of Bessel functions.[20][21] For instance,

${\displaystyle {\text{fe}}_{2n}(x,q)=-\pi n\sum _{r=0}^{\infty }(-1)^{r+n}A_{2r}^{(2n)}{\text{Im}}[J_{r}({\sqrt {q}}e^{ix})]{\text{Im}}[Y_{r}({\sqrt {q}}e^{-ix})]}$

where ${\displaystyle n,q>0}$ , and ${\displaystyle J_{r}(x)}$  and ${\displaystyle Y_{r}(x)}$  are Bessel functions of the first and second kind, respectively.

## Properties

There are relatively few analytic expressions and identities involving Mathieu functions. Moreover, unlike many other special functions, the solutions of Mathieu's equation cannot in general be expressed in terms of hypergeometric functions. This can be seen by transformation of Mathieu's equation to algebraic form, using the change of variable ${\displaystyle t=\cos(x)}$ :

${\displaystyle (1-t^{2}){\frac {d^{2}y}{dt^{2}}}-t\,{\frac {dy}{dt}}+(a+2q(1-2t^{2}))\,y=0.}$

Since this equation has an irregular singular point at infinity, it cannot be transformed into an equation of the hypergeometric type.[18]

### Qualitative behavior

Sample plots of Mathieu functions of the first kind

Plot of ${\displaystyle {\text{ce}}_{1}(x,q)}$  for varying ${\displaystyle q}$

For small ${\displaystyle q}$ , ${\displaystyle {\text{ce}}_{n}}$  and ${\displaystyle {\text{se}}_{n}}$  behave similarly to ${\displaystyle \cos nx}$  and ${\displaystyle \sin nx}$ . For arbitrary ${\displaystyle q}$ , they may deviate significantly from their trigonometric counterparts; however, they remain periodic in general. Moreover, for any real ${\displaystyle q}$ , ${\displaystyle {\text{ce}}_{m}(x,q)}$  and ${\displaystyle {\text{se}}_{m+1}(x,q)}$  have exactly ${\displaystyle m}$  simple zeros in ${\displaystyle 0 , and as ${\displaystyle q\rightarrow \infty }$  the zeros cluster about ${\displaystyle x=\pi /2}$ .[22][23]

For ${\displaystyle q>0}$  and as ${\displaystyle x\rightarrow \infty }$  the modified Mathieu functions tend to behave as damped periodic functions.

In the following, the ${\displaystyle A}$  and ${\displaystyle B}$  factors from the Fourier expansions for ${\displaystyle {\text{ce}}_{n}}$  and ${\displaystyle {\text{se}}_{n}}$  may be referenced (see Explicit representation and computation). They depend on ${\displaystyle q}$  and ${\displaystyle n}$  but are independent of ${\displaystyle x}$ .

### Reflections and translations

Due to their parity and periodicity, ${\displaystyle {\text{ce}}_{n}}$  and ${\displaystyle {\text{se}}_{n}}$  have simple properties under reflections and translations by multiples of ${\displaystyle \pi }$ :[6]

{\displaystyle {\begin{aligned}&{\text{ce}}_{n}(x+\pi )=(-1)^{n}{\text{ce}}_{n}(x)\\&{\text{se}}_{n}(x+\pi )=(-1)^{n}{\text{se}}_{n}(x)\\&{\text{ce}}_{n}(x+\pi /2)=(-1)^{n}{\text{ce}}_{n}(-x+\pi /2)\\&{\text{se}}_{n+1}(x+\pi /2)=(-1)^{n}{\text{se}}_{n+1}(-x+\pi /2)\end{aligned}}}

One can also write functions with negative ${\displaystyle q}$  in terms of those with positive ${\displaystyle q}$ :[4][24]

{\displaystyle {\begin{aligned}&{\text{ce}}_{2n}(x,-q)=(-1)^{n}{\text{ce}}_{2n}(-x+\pi /2,q)\\&{\text{ce}}_{2n+1}(x,-q)=(-1)^{n}{\text{se}}_{2n+1}(-x+\pi /2,q)\\&{\text{se}}_{2n+1}(x,-q)=(-1)^{n}{\text{ce}}_{2n+1}(-x+\pi /2,q)\\&{\text{se}}_{2n+2}(x,-q)=(-1)^{n}{\text{se}}_{2n+2}(-x+\pi /2,q)\end{aligned}}}

Moreover,

{\displaystyle {\begin{aligned}&a_{2n+1}(q)=b_{2n+1}(-q)\\&b_{2n+2}(q)=b_{2n+2}(-q)\end{aligned}}}

### Orthogonality and completeness

Like their trigonometric counterparts ${\displaystyle \cos nx}$  and ${\displaystyle \sin nx}$ , the periodic Mathieu functions ${\displaystyle {\text{ce}}_{n}(x,q)}$  and ${\displaystyle {\text{se}}_{n}(x,q)}$  satisfy orthogonality relations

{\displaystyle {\begin{aligned}&\int _{0}^{2\pi }{\text{ce}}_{n}{\text{ce}}_{m}dx=\int _{0}^{2\pi }{\text{se}}_{n}{\text{se}}_{m}dx=\delta _{nm}\pi \\&\int _{0}^{2\pi }{\text{ce}}_{n}{\text{se}}_{m}dx=0\end{aligned}}}

Moreover, with ${\displaystyle q}$  fixed and ${\displaystyle a}$  treated as the eigenvalue, the Mathieu equation is of Sturm-Liouville form. This implies that the eigenfunctions ${\displaystyle {\text{ce}}_{n}(x,q)}$  and ${\displaystyle {\text{se}}_{n}(x,q)}$  form a complete set, i.e. any ${\displaystyle \pi }$ - or ${\displaystyle 2\pi }$ -periodic function of ${\displaystyle x}$  can be expanded as a series in ${\displaystyle {\text{ce}}_{n}(x,q)}$  and ${\displaystyle {\text{se}}_{n}(x,q)}$ .[3]

### Integral identities

Solutions of Mathieu's equation satisfy a class of integral identities with respect to kernels ${\displaystyle \chi (x,x')}$  that are solutions of

${\displaystyle {\frac {\partial ^{2}\chi }{\partial x^{2}}}-{\frac {\partial ^{2}\chi }{\partial x'^{2}}}=2q\left(\cos 2x-\cos 2x'\right)\chi }$

More precisely, if ${\displaystyle \phi (x)}$  solves Mathieu's equation with given ${\displaystyle a}$  and ${\displaystyle q}$ , then the integral

${\displaystyle \psi (x)\equiv \int _{C}\chi (x,x')\phi (x')dx'}$

where ${\displaystyle C}$  is a path in the complex plane, also solves Mathieu's equation with the same ${\displaystyle a}$  and ${\displaystyle q}$ , provided the following conditions are met:[25]

• ${\displaystyle \chi (x,x')}$  solves ${\displaystyle {\frac {\partial ^{2}\chi }{\partial x^{2}}}-{\frac {\partial ^{2}\chi }{\partial x'^{2}}}=2q\left(\cos 2x-\cos 2x'\right)\chi }$
• In the regions under consideration, ${\displaystyle \psi (x)}$  exists and ${\displaystyle \chi (x,x')}$  is analytic
• ${\displaystyle \left(\phi {\frac {\partial \chi }{\partial x'}}-{\frac {\partial \phi }{\partial x'}}\chi \right)}$  has the same value at the endpoints of ${\displaystyle C}$

Using an appropriate change of variables, the equation for ${\displaystyle \chi }$  can be transformed into the wave equation and solved. For instance, one solution is ${\displaystyle \chi (x,x')=\sinh(2q^{1/2}\sin x\sin x')}$ . Examples of identities obtained in this way are[26]

{\displaystyle {\begin{aligned}{\text{se}}_{2n+1}(x,q)&={\frac {{\text{se}}'_{2n+1}(0,q)}{\pi q^{1/2}B_{1}^{(2n+1)}}}\int _{0}^{\pi }\sinh(2q^{1/2}\sin x\sin x'){\text{se}}_{2n+1}(x',q)dx'\qquad (q>0)\\{\text{Ce}}_{2n}(x,q)&={\frac {{\text{ce}}_{2n}(\pi /2,q)}{\pi A_{0}^{(2n)}}}\int _{0}^{\pi }\cos(2q^{1/2}\cosh x\cos x'){\text{ce}}_{2n}(x',q)dx'\qquad \ \ \ (q>0)\end{aligned}}}

Identities of the latter type are useful for studying asymptotic properties of the modified Mathieu functions.[27]

### Asymptotic expansions

The following asymptotic expansions hold for ${\displaystyle q>0}$ , ${\displaystyle {\text{Im}}(x)=0}$ , ${\displaystyle {\text{Re}}(x)\rightarrow \infty }$ , and ${\displaystyle 2q^{1/2}\cosh x\simeq q^{1/2}e^{x}}$ :[28]

{\displaystyle {\begin{aligned}{\text{Ce}}_{2n}(x,q)&\sim \left({\frac {2}{\pi q^{1/2}}}\right)^{1/2}{\frac {{\text{ce}}_{2n}(0,q){\text{ce}}_{2n}(\pi /2,q)}{A_{0}^{(2n)}}}\cdot e^{-x/2}\sin \left(q^{1/2}e^{x}+{\frac {\pi }{4}}\right)\\{\text{Ce}}_{2n+1}(x,q)&\sim \left({\frac {2}{\pi q^{3/2}}}\right)^{1/2}{\frac {{\text{ce}}_{2n+1}(0,q){\text{ce}}'_{2n+1}(\pi /2,q)}{A_{1}^{(2n+1)}}}\cdot e^{-x/2}\cos \left(q^{1/2}e^{x}+{\frac {\pi }{4}}\right)\\{\text{Se}}_{2n+1}(x,q)&\sim -\left({\frac {2}{\pi q^{3/2}}}\right)^{1/2}{\frac {{\text{se}}'_{2n+1}(0,q){\text{se}}_{2n+1}(\pi /2,q)}{B_{1}^{(2n+1)}}}\cdot e^{-x/2}\cos \left(q^{1/2}e^{x}+{\frac {\pi }{4}}\right)\\{\text{Se}}_{2n+2}(x,q)&\sim \left({\frac {2}{\pi q^{5/2}}}\right)^{1/2}{\frac {{\text{se}}'_{2n+2}(0,q){\text{se}}'_{2n+2}(\pi /2,q)}{B_{2}^{(2n+2)}}}\cdot e^{-x/2}\sin \left(q^{1/2}e^{x}+{\frac {\pi }{4}}\right)\end{aligned}}}

Thus, the modified Mathieu functions decay exponentially for large real argument. Similar asymptotic expansions can be written down for ${\displaystyle {\text{Fe}}_{n}}$  and ${\displaystyle {\text{Ge}}_{n}}$ ; these also decay exponentially for large real argument.

One can also derive asymptotic expansions for large ${\displaystyle q}$ .[29][30] For the characteristic numbers in particular, one has

${\displaystyle a=-2|q|+{\mathcal {O}}|q|^{1/2}}$

## Applications

Mathieu's differential equations appear in a wide range of contexts in engineering, physics, and applied mathematics. Many of these applications fall into one of two general categories: 1) the analysis of partial differential equations in elliptic geometries, and 2) dynamical problems which involve forces that are periodic in either space or time. Examples within both categories are discussed below.

### Partial differential equations

Mathieu functions arise when separation of variables in elliptic coordinates is applied to 1) the Laplace equation in 3 dimensions, and 2) the Helmholtz equation in either 2 or 3 dimensions. Since the Helmholtz equation is a prototypical equation for modeling the spatial variation of classical waves, Mathieu functions can be used to describe a variety of wave phenomena. For instance, in computational electromagnetics they can be used to analyze the scattering of electromagnetic waves off elliptic cylinders, and wave propagation in elliptic waveguides.[31] In general relativity, an exact plane wave solution to the Einstein field equation can be given in terms of Mathieu functions.

More recently, Mathieu functions have been used to solve a special case of the Smoluchowski equation, describing the steady-state statistics of self-propelled particles.[32]

The remainder of this section details the analysis for the two-dimensional Helmholtz equation.[33] In rectangular coordinates, the Helmholtz equation is

${\displaystyle \left({\frac {\partial ^{2}}{\partial x^{2}}}+{\frac {\partial ^{2}}{\partial y^{2}}}\right)\psi +k^{2}\psi =0,}$

Elliptic coordinates in this case are defined by

{\displaystyle {\begin{aligned}x&=c\cosh \mu \cos \nu \\y&=c\sinh \mu \sin \nu \end{aligned}}}

where ${\displaystyle 0\leq \mu <\infty }$ , ${\displaystyle 0\leq \nu <2\pi }$ , and ${\displaystyle c}$  is a positive constant. The Helmholtz equation in these coordinates is

${\displaystyle {\frac {1}{c^{2}(\sinh ^{2}\mu +\sin ^{2}\nu )}}\left({\frac {\partial ^{2}}{\partial \mu ^{2}}}+{\frac {\partial ^{2}}{\partial \nu ^{2}}}\right)\psi +k^{2}\psi =0}$

The constant ${\displaystyle \mu }$  curves are confocal ellipses with focal length ${\displaystyle c}$ ; hence, these coordinates are convenient for solving the Helmholtz equation on domains with elliptic boundaries. Separation of variables via ${\displaystyle \psi (\mu ,\nu )=F(\mu )G(\nu )}$  yields the Mathieu equations

{\displaystyle {\begin{aligned}&{\frac {d^{2}F}{d\mu ^{2}}}-\left(a-{\frac {c^{2}k^{2}}{2}}\cosh 2\mu \right)F=0\\&{\frac {d^{2}G}{d\nu ^{2}}}+\left(a-{\frac {c^{2}k^{2}}{2}}\cos 2\nu \right)G=0\\\end{aligned}}}

where ${\displaystyle a}$  is a separation constant.

As a specific physical example, the Helmholtz equation can be interpreted as describing normal modes of an elastic membrane under uniform tension. In this case, the following physical conditions are imposed:[34]

• Periodicity with respect to ${\displaystyle \nu }$ , i.e. ${\displaystyle \psi (\mu ,\nu )=\psi (\mu ,\nu +2\pi )}$
• Continuity of displacement across the interfocal line: ${\displaystyle \psi (0,\nu )=\psi (0,-\nu )}$
• Continuity of derivative across the interfocal line: ${\displaystyle \psi _{\mu }(0,\nu )=-\psi _{\mu }(0,-\nu )}$

For given ${\displaystyle k}$ , this restricts the solutions to those of the form ${\displaystyle {\text{Ce}}_{n}(\mu ,q){\text{ce}}_{n}(\nu ,q)}$  and ${\displaystyle {\text{Se}}_{n}(\mu ,q){\text{se}}_{n}(\nu ,q)}$ , where ${\displaystyle q=c^{2}k^{2}/2}$ . This is the same as restricting allowable values of ${\displaystyle a}$ , for given ${\displaystyle k}$ . Restrictions on ${\displaystyle k}$  then arise due to imposition of physical conditions on some bounding surface, such as an elliptic boundary defined by ${\displaystyle \mu =\mu _{0}>0}$ . For instance, clamping the membrane at ${\displaystyle \mu =\mu _{0}}$  imposes ${\displaystyle \psi (\mu _{0},\nu )=0}$ , which in turn requires

{\displaystyle {\begin{aligned}{\text{Ce}}_{n}(\mu _{0},q)=0\\{\text{Se}}_{n}(\mu _{0},q)=0\end{aligned}}}

These conditions define the normal modes of the system.

### Dynamical problems

In dynamical problems with periodically varying forces, the equation of motion sometimes takes the form of Mathieu's equation. In such cases, knowledge of the general properties of Mathieu's equation— particularly with regard to stability of the solutions—can be essential for understanding qualitative features of the physical dynamics.[35] A classic example along these lines is the inverted pendulum.[36] Other examples are

### Quantum mechanics

Mathieu functions play a role in certain quantum mechanical systems, particularly those with spatially periodic potentials such as the quantum pendulum and crystalline lattices.

The modified Mathieu equation also arises when describing the quantum mechanics of singular potentials. For the particular singular potential ${\displaystyle V(r)=g^{2}/r^{4}}$  the radial Schrödinger equation

${\displaystyle {\frac {d^{2}y}{dr^{2}}}+\left[k^{2}-{\frac {\ell (\ell +1)}{r^{2}}}-{\frac {g^{2}}{r^{4}}}\right]y=0}$

can be converted into the equation

${\displaystyle {\frac {d^{2}\varphi }{dz^{2}}}+\left[2h^{2}\cosh 2z-\left(\ell +{\frac {1}{2}}\right)^{2}\right]\varphi =0.}$

The transformation is achieved with the following substitutions

${\displaystyle y=r^{1/2}\varphi ,r=\gamma e^{z},\gamma ={\frac {ig}{h}},h^{2}=ikg,h=e^{I\pi /4}(kg)^{1/2}.}$

By solving the Schrödinger equation (for this particular potential) in terms of solutions of the modified Mathieu equation, scattering properties such as the S-matrix and the absorptivity can be obtained.[37]

## Notes

1. ^ Mathieu (1868).
2. ^ Morse and Feshbach (1953).
3. ^ a b Gutiérrez-Vega (2015).
4. Arscott (1964), chapter III
5. ^ Arscott (1964) 43–44
6. ^ a b c McLachlan (1947), chapter II.
7. ^ Arscott (1964); Iyanaga (1980); Gradshteyn(2007); This is also the normalization used by the computer algebra system Maple, see https://www.maplesoft.com/support/help/Maple/view.aspx?path=Mathieu.
8. ^ Arscott (1964), p. 29.
9. ^ It is not true, in general, that a ${\displaystyle 2\pi }$  periodic function has the property ${\displaystyle y(x+\pi )=-y(x)}$ . However, this turns out to be true for functions which are solutions of Mathieu's equation.
10. ^ McLachlan (1951), pp. 141-157, 372
11. ^ Arscott (1964), p. 34
12. ^ McLachlan (1947), p. 144
13. ^ McLachlan (1947), p. 372
14. ^ McLachlan (1947), p. 28
15. ^ a b Wimp (1984), pp. 83-84
16. ^ McLachlan (1947)
17. ^ Chaos-Cador and Ley-Koo (2001)
18. ^ a b Temme (2015), p. 234
19. ^ Müller-Kirsten, Harald J.W., Introduction to Quantum Mechanics: Schrödinger Equation and Path Integral, 2nd ed. (2012), World Scientific, pp. 420-428.
20. ^ McLachlan (1947), p. 254
21. ^ Malits (2010)
22. ^ Meixner and Schäfke (1954), p.134
23. ^ McLachlan (1947), pp. 234–235
24. ^ Gradshteyn (2007), p. 953
25. ^ Arscott (1964), pp. 40-41
26. ^ Gradshteyn (2007), pp. 763–765
27. ^ Arscott (1964), p. 86
28. ^ McLachlan (1947), chapter XI
29. ^ McLachlan (1947), p. 237
30. ^ R.B. Dingle and H. J. W. Müller-Kirsten, Asymptotic Expansions of Mathieu Functions and their Characteristic Numbers, J. reine u. angew. Math. 211 (1962) 11 - 32; H. J. W. Müller, On Asymptotic Expansions of Mathieu Functions, J. reine u. angew. Math. 211 (1962) 179 - 190.
31. ^ Bibby and Peterson (2013); Barakat (1963); Sebak and Shafai (1991); Kretzschmar (1970)
32. ^ Solon et al (2015)
33. ^ see Willatzen and Voon (2011), pp. 61–65
34. ^ McLachlan (1947), pp. 294–297
35. ^ a b Meixner and Schäfke (1954), pp. 324–343
36. ^ Ruby (1996)
37. ^ Müller-Kirsten (2006)

## References

• Arscott, Felix (1964). Periodic differential equations: an introduction to Mathieu, Lamé, and allied functions. Pergamon Press. ISBN 9781483164885.
• Barakat, R. (1963), "Diffraction of Plane Waves by an Elliptic Cylinder", The Journal of the Acoustical Society of America, 35 (12): 1990–1996, doi:10.1121/1.1918878
• Bibby, Malcom M.; Peterson, Andrew F. (2014). Accurate Computation of Mathieu Functions. Morgan & Claypool. ISBN 9781627050852.
• Chaos-Cador, L.; Ley-Koo, E. (2002), "Mathieu functions revisited: matrix evaluation and generating functions", Revista mexicana de física, 48 (1): 67–75
• Gradshteyn, Izrail Solomonovich; et al. (February 2007). Jeffrey, Alan; Zwillinger, Daniel (eds.). Table of Integrals, Series, and Products. Translated by Scripta Technica, Inc. (7 ed.). Academic Press, Inc. ISBN 978-0-12-373637-6. MR 2360010.
• Gutiérrez-Vega, Julio C. (2015), "Mathieu Functions", in Nicholas J. Higham; et al. (eds.), The Princeton Companion to Applied Mathematics, Princeton University Press, pp. 159–160
• Iyanaga, Shōkichi; Kawada, Yukiyosi, eds. (1980) [1977]. Encyclopedic Dictionary of Mathematics, Volume I. Translated from the 2nd Japanese edition, paperback version of the 1977 edition (1st ed.). MIT Press. ISBN 978-0-262-59010-5. MR 0591028.
• Kretzschmar, J.G. (1970), "Wave Propagation in Hollow Conducting Elliptical Waveguides", IEEE Transactions on Microwave Theory and Techniques, 18 (9): 547–554, doi:10.1109/TMTT.1970.1127288
• Malits, Pinchas (2010), "Relations between Mathieu functions of the first and second kind", Integral Transforms and Special Functions, 21 (6): 423–436, doi:10.1080/10652460903360499
• Mathieu, E. (1868), "Mémoire sur Le Mouvement Vibratoire d'une Membrane de forme Elliptique", Journal de Mathématiques Pures et Appliquées: 137–203
• McLachlan, N. W. (1951). Theory and application of Mathieu functions. Oxford University Press. Note: Reprinted lithographically in Great Britain at the University Press, Oxford, 1951 from corrected sheets of the (1947) first edition.
• Meixner, Josef; Schäfke, Friedrich Wilhelm (1954). Mathieusche Funktionen und Sphäroidfunktionen (in German). Berlin: Springer-Verlag.
• Morse, Philip McCord; Feshbach, Herman (1953-01-01). Methods of Theoretical Physics: Pt. 1 (Reprint ed.). Boston, Mass: McGraw-Hill Inc.,US. ISBN 9780070433168.
• Müller-Kirsten, Harald J.W. (2006). Introduction to Quantum Mechanics: Schrödinger Equation and Path Integral. World Scientific. ISBN 981-256-692-9.
• Sebak, A.; Shafai, L. (1991), "Generalized solutions for electromagnetic scattering by elliptical structures", Computer Physics Communications, 68 (1–3): 315–330, doi:10.1016/0010-4655(91)90206-Z
• Solon, A.P.; Cates, M.E.; Tailleur, J. (2015), "Active brownian particles and run-and-tumble particles: A comparative study", The European Physical Journal Special Topics, 224 (7): 1231–1262, arXiv:1504.07391, doi:10.1140/epjst/e2015-02457-0
• Temme, Nico M. (2015), "Special Functions", in Nicholas J. Higham; et al. (eds.), The Princeton Companion to Applied Mathematics, Princeton University Press, p. 234
• Lew Yan Voon LC, Willatzen M (2011). Separable Boundary-Value Problems in Physics. Wiley-VCH. doi:10.1002/9783527634927. ISBN 978-3-527-41020-0. (free online access to the appendix on Mathieu functions)
• Wimp, Jet (1984). Computation with Recurrence Relations. Pitman Publishing. pp. 83–84. ISBN 0-273-08508-5.
• Wolf, G. (2010), "Mathieu Functions and Hill's Equation", in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0521192255, MR 2723248