28 Mathieu Functions and Hill’s EquationMathieu Functions of Integer Order

§28.6 Expansions for Small q

Contents
  1. §28.6(i) Eigenvalues
  2. §28.6(ii) Functions cen and sen

§28.6(i) Eigenvalues

Leading terms of the power series for am⁡(q) and bm⁡(q) for m≤6 are:

28.6.1 a0⁡(q) =−12⁢q2+7128⁢q4−292304⁢q6+68687188 74368⁢q8+⋯,
28.6.2 a1⁡(q) =1+q−18⁢q2−164⁢q3−11536⁢q4+1136864⁢q5+495 89824⁢q6+5594 37184⁢q7−83353 89440⁢q8+⋯,
28.6.3 b1⁡(q) =1−q−18⁢q2+164⁢q3−11536⁢q4−1136864⁢q5+495 89824⁢q6−5594 37184⁢q7−83353 89440⁢q8+⋯,
28.6.4 a2⁡(q) =4+512⁢q2−76313824⁢q4+10 02401796 26240⁢q6−16690 6840145 86471 42400⁢q8+⋯,
28.6.5 b2⁡(q) =4−112⁢q2+513824⁢q4−289796 26240⁢q6+2139145 86471 42400⁢q8+⋯,
28.6.6 a3⁡(q) =9+116⁢q2+164⁢q3+1320480⁢q4−516384⁢q5−1961235 92960⁢q6−6091048 57600⁢q7+⋯,
28.6.7 b3⁡(q) =9+116⁢q2−164⁢q3+1320480⁢q4+516384⁢q5−1961235 92960⁢q6+6091048 57600⁢q7+⋯,
28.6.8 a4⁡(q) =16+130⁢q2+4338 64000⁢q4−570127216 00000⁢q6+⋯,
28.6.9 b4⁡(q) =16+130⁢q2−3178 64000⁢q4+1004927216 00000⁢q6+⋯,
28.6.10 a5⁡(q) =25+148⁢q2+117 74144⁢q4+11 47456⁢q5+378918 13888⁢q6+⋯,
28.6.11 b5⁡(q) =25+148⁢q2+117 74144⁢q4−11 47456⁢q5+378918 13888⁢q6+⋯,
28.6.12 a6⁡(q) =36+170⁢q2+187439 04000⁢q4+67 436179293 59872 00000⁢q6+⋯,
28.6.13 b6⁡(q) =36+170⁢q2+187439 04000⁢q4−58 616339293 59872 00000⁢q6+⋯.

Leading terms of the of the power series for m=7,8,9,… are:

28.6.14 am⁡(q)bm⁡(q)}=m2+12⁢(m2−1)⁢q2+5⁢m2+732⁢(m2−1)3⁢(m2−4)⁢q4+9⁢m4+58⁢m2+2964⁢(m2−1)5⁢(m2−4)⁢(m2−9)⁢q6+⋯.

For more details on these expansions and recurrence relations for the coefficients see Frenkel and Portugal (2001, §2).

The coefficients of the power series of a2⁢n⁡(q), b2⁢n⁡(q) and also a2⁢n+1⁡(q), b2⁢n+1⁡(q) are the same until the terms in q2⁢n−2 and q2⁢n, respectively. Then

28.6.15 am⁡(q)−bm⁡(q)=2⁢qm(2m−1⁢(m−1)!)2⁢(1+O⁡(q2)).

Higher coefficients in the foregoing series can be found by equating coefficients in the following continued-fraction equations:

28.6.16 a−(2⁢n)2−q2a−(2⁢n−2)2−q2a−(2⁢n−4)2−⁢⋯⁢q2a−22−2⁢q2a=−q2(2⁢n+2)2−a−q2(2⁢n+4)2−a−⋯,
a=a2⁢n⁡(q),
28.6.17 a−(2⁢n+1)2−q2a−(2⁢n−1)2−⁢⋯⁢q2a−32−q2a−12−q=−q2(2⁢n+3)2−a−q2(2⁢n+5)2−a−⋯,
a=a2⁢n+1⁡(q),
28.6.18 a−(2⁢n+1)2−q2a−(2⁢n−1)2−⁢⋯⁢q2a−32−q2a−12+q=−q2(2⁢n+3)2−a−q2(2⁢n+5)2−a−⋯,
a=b2⁢n+1⁡(q),
28.6.19 a−(2⁢n+2)2−q2a−(2⁢n)2−q2a−(2⁢n−2)2−⁢⋯⁢q2a−22=−q2(2⁢n+4)2−a−q2(2⁢n+6)2−a−⋯,
a=b2⁢n+2⁡(q).

Numerical values of the radii of convergence ρn(j) of the power series (28.6.1)–(28.6.14) for n=0,1,…,9 are given in Table 28.6.1. Here j=1 for a2⁢n⁡(q), j=2 for b2⁢n+2⁡(q), and j=3 for a2⁢n+1⁡(q) and b2⁢n+1⁡(q). (Table 28.6.1 is reproduced from Meixner et al. (1980, §2.4).)

Table 28.6.1: Radii of convergence for power-series expansions of eigenvalues of Mathieu’s equation.
n ρn(1) ρn(2) ρn(3)
0 or 1 1.46876 86138 6.92895 47588 3.76995 74940
2 7.26814 68935 16.80308 98254 11.27098 52655
3 16.47116 58923 30.09677 28376 22.85524 71216
4 30.42738 20960 48.13638 18593 38.52292 50099
5 47.80596 57026 69.59879 32769 58.27413 84472
6 69.92930 51764 95.80595 67052 82.10894 36067
7 95.47527 27072 125.43541 1314 110.02736 9210
8 125.76627 89677 159.81025 4642 142.02943 1279
9 159.47921 26694 197.60667 8692 178.11513 940

It is conjectured that for large n, the radii increase in proportion to the square of the eigenvalue number n; see Meixner et al. (1980, §2.4). It is known that

28.6.20 lim infn→∞ρn(j)n2≥k⁢k′⁢(K⁡(k))2=2.04183 4⁢…,

where k is the unique root of the equation 2⁢E⁡(k)=K⁡(k) in the interval (0,1), and k′=1−k2. For E⁡(k) and K⁡(k) see §19.2(ii).

§28.6(ii) Functions cen and sen

Leading terms of the power series for the normalized functions are:

28.6.21 21/2⁢ce0⁡(z,q) =1−12⁢q⁢cos⁡2⁢z+132⁢q2⁢(cos⁡4⁢z−2)−1128⁢q3⁢(19⁢cos⁡6⁢z−11⁢cos⁡2⁢z)+⋯,
28.6.22 ce1⁡(z,q) =cos⁡z−18⁢q⁢cos⁡3⁢z+1128⁢q2⁢(23⁢cos⁡5⁢z−2⁢cos⁡3⁢z−cos⁡z)−11024⁢q3⁢(19⁢cos⁡7⁢z−89⁢cos⁡5⁢z−13⁢cos⁡3⁢z+2⁢cos⁡z)+⋯,
28.6.23 se1⁡(z,q) =sin⁡z−18⁢q⁢sin⁡3⁢z+1128⁢q2⁢(23⁢sin⁡5⁢z+2⁢sin⁡3⁢z−sin⁡z)−11024⁢q3⁢(19⁢sin⁡7⁢z+89⁢sin⁡5⁢z−13⁢sin⁡3⁢z−2⁢sin⁡z)+⋯,
28.6.24 ce2⁡(z,q) =cos⁡2⁢z−14⁢q⁢(13⁢cos⁡4⁢z−1)+1128⁢q2⁢(13⁢cos⁡6⁢z−769⁢cos⁡2⁢z)+⋯,
28.6.25 se2⁡(z,q) =sin⁡2⁢z−112⁢q⁢sin⁡4⁢z+1128⁢q2⁢(13⁢sin⁡6⁢z−49⁢sin⁡2⁢z)+⋯.

For m=3,4,5,…,

28.6.26 cem⁡(z,q)=cos⁡m⁢z−q4⁢(1m+1⁢cos⁡(m+2)⁢z−1m−1⁢cos⁡(m−2)⁢z)+q232⁢(1(m+1)⁢(m+2)⁢cos⁡(m+4)⁢z+1(m−1)⁢(m−2)⁢cos⁡(m−4)⁢z−2⁢(m2+1)(m2−1)2⁢cos⁡m⁢z)+⋯.

For the corresponding expansions of sem⁡(z,q) for m=3,4,5,… change cos to sin everywhere in (28.6.26). For more details on these expansions and recurrence relations for the coefficients see Frenkel and Portugal (2001, §2).’

The radii of convergence of the series (28.6.21)–(28.6.26) are the same as the radii of the corresponding series for an⁡(q) and bn⁡(q); compare Table 28.6.1 and (28.6.20).