The Mathieu functions MathieuC(a, q, x) and MathieuS(a, q, x) are solutions of the Mathieu differential equation:

MathieuC and MathieuS are even and odd functions of x, respectively.

For countably many values of a (as a function of q), MathieuC and MathieuS are 2*Pi-periodic. For $n\=0,1,2,...$, MathieuA(n, q) is the nth such characteristic value for MathieuC, and for $n\=1,...$, MathieuB(n, q) is the nth characteristic value for MathieuS. The resulting Mathieu functions are:

where $\mathrm{c1}$ and $\mathrm{c2}$ are normalization constants depending on n and q.

If the index n is even, then both MathieuCE and MathieuSE are Pi-periodic; they are 2*Pi-periodic otherwise. MathieuCE and MathieuSE are even and odd functions of x, respectively.

MathieuFloquet(a, q, x) is a Floquet solution of Mathieu's equation. It has the form:

where $\mathrm{\nu }=\mathrm{MathieuExponent}\left(a\,q\right)$ is the characteristic exponent and $P\left(x\right)$ is a Pi periodic function.

MathieuCPrime, MathieuSPrime, MathieuCEPrime, MathieuSEPrime, and MathieuFloquetPrime are the first derivatives with respect to x of the corresponding Mathieu functions. Note that all higher order derivatives can be written in terms of the 0th and 1st derivatives.

The odd and even Mathieu functions are related to the Floquet solution via:

The various Mathieu functions are normalized as follows.

where P is as given in the definition of the Floquet solution.

The normalizations of MathieuCE, MathieuSE, and their derivatives coincide with the ones in [1] (see references below), except for $n=0$, where the function implemented in Maple is equal to the function in reference [1] multiplied by sqrt(2). In this way, we have, in Maple, for instance, MathieuCE(0, 1, 0) = 1 while in [1] the right-hand side is equal to $\sqrt{2}$. For details on this normalization, see reference [2].

MathieuExponent is an inverse to both MathieuA and MathieuB in the following sense.

where $k$ is an integer.

For q = 0, the Mathieu functions assume special values:

$\mathrm{MathieuC}\left(1.23\,4.56\,7.89\right)$

$147.6289914$

$\mathrm{MathieuFloquet}\left(0.1\,0.2\,0.3\right)$

$0.8965972662+0.1161779284\mathrm{I}$

$\mathrm{MathieuS}\left(a\,q\,-x\right)$

$-\mathrm{MathieuS}\left(a\,q\,x\right)$

$\mathrm{MathieuCE}\left(n\,0\,x\right)$

$\cos \left(nx\right)$

$\mathrm{MathieuExponent}\left(a\,0\right)$

$\sqrt{a}$

$\mathrm{MathieuSE}\left(n\,q\,x+\frac{25}{11}\mathrm{\pi }\right)$

$\mathrm{MathieuSE}\left(n\,q\,x+\frac{3\mathrm{\pi }}{11}\right)$

$\mathrm{diff}\left(\mathrm{MathieuC}\left(a\,q\,x\right)\,x\right)$

$\mathrm{MathieuCPrime}\left(a\,q\,x\right)$

$\mathrm{diff}\left(\,x\right)$

$\left(2q\cos \left(2x\right)-a\right)\mathrm{MathieuC}\left(a\,q\,x\right)$

$\mathrm{diff}\left(\mathrm{MathieuCE}\left(n\,q\,x\right)\,x\,x\,x\right)$

$-4q\sin \left(2x\right)\mathrm{MathieuCE}\left(n\,q\,x\right)+\left(2q\cos \left(2x\right)-\mathrm{MathieuA}\left(n\,q\right)\right)\mathrm{MathieuCEPrime}\left(n\,q\,x\right)$

$\mathrm{series}\left(\mathrm{MathieuA}\left(1\,q\right)\,q\right)$

$1+q-\frac{1}{8}{q}^2-\frac{1}{64}{q}^3-\frac{1}{1536}{q}^4+\frac{11}{36864}{q}^5+O\left(q^6\right)$

$\mathrm{series}\left(\mathrm{MathieuSE}\left(2\,q\,x\right)\,q\,3\right)$

$\sin \left(2x\right)-\frac{1}{12}\sin \left(4x\right)q+\left(-\frac{\sin \left(2x\right)}{288}+\frac{\sin \left(6x\right)}{384}\right)q^2+O\left(q^3\right)$

$\mathrm{series}\left(\mathrm{MathieuCPrime}\left(2\,q\,x\right)\,x\,5\right)$

$\left(2q-2\right)x+\left(\frac{2}{3}{q}^2-\frac{8}{3}q+\frac{2}{3}\right)x^3+O\left(x^5\right)$

[1] Abramowitz, M., and Stegun, I., eds. Handbook of Mathematical Functions. New York: Dover Publications.

[2] Frenkel, D., and Portugal, R. "Algebraic methods to compute Mathieu functions." Journal of Physics A: Mathematical and General, Vol. 34. (2001): 3541-3551.

[3] McLachlan, N. W. Theory and Applications of Mathieu Functions. Oxford University Press.