- [Kre, p.94, (27.15)].

Suppose A, B, and C are the first, second and third determinants of [Kre, p.94, l.7] respectively. Let D and E be the first and second determinants of [Kre, p.94, (27.15)] respectively. We want to prove AB - C^{2}= D - E. Our strategy is to use the following methods to reduce calculations:- (Using symmetry) E is a 3´3 determinant. Each entry has three
terms. Therefore, E contains 3
^{3}´6 terms. Using symmetry we may write

E = (x_{1}^{2}+y_{1}^{2}+z_{1}^{2})(x_{2}^{2}+y_{2}^{2}+z_{2}^{2})(x_{12}^{2}+y_{12}^{2}+z_{12}^{2})

+2(x_{12}x_{1}+y_{12}y_{1}+z_{12}z_{1})(x_{1}x_{2}+y_{1}y_{2}+z_{1}z_{2})(x_{2}x_{12}+y_{2}y_{12}+z_{2}z_{12})

-(x_{12}x_{1}+y_{12}y_{1}+z_{12}z_{1})^{2}(x_{2}^{2}+y_{2}^{2}+z_{2}^{2})

-(x_{1}x_{2}+y_{1}y_{2}+z_{1}z_{2})^{2}(x_{12}^{2}+y_{12}^{2}+z_{12}^{2})

-(x_{1}^{2}+y_{1}^{2}+z_{1}^{2})(x_{12}x_{2}+y_{12}y_{2}+z_{12}z_{2})^{2}.

- (Using cancellations) Cancellations reduce E to 21 terms. Similarly, cancellations reduce D to 33
terms.

The expansion of AB - C^{2}contains 57 terms. - (Using signs) We divide the 54 (=21+33) terms on the right-hand side into two groups according to the sign of each term. Choose any term from the left-hand side and check to see if the term is also on the right-hand side. If the sign of the term is +, then we need not check the group with minus signs.
- (Using factorizations)

Obviously, the term x_{2}z_{1}y_{1}z_{2}y_{11}y_{22}does not appear in (x_{1}^{2}+y_{1}^{2}+z_{1}^{2})(x_{2}^{2}+y_{2}^{2}+z_{2}^{2})(x_{11}x_{22}+y_{11}y_{22}+z_{11}z_{22}) (a part of the expansion of D).

- (Using symmetry) E is a 3´3 determinant. Each entry has three
terms. Therefore, E contains 3
- [Ber, p. 49, (4.1.2)].

Proof. We preserve intermediate steps in order to reveal cancellations. Let

x^{1}= r sinq cosf,

x^{2}= r sinq sinf,

and x^{3}= r cosq.

dx^{1}= dr sinq cosf + r d (sinq cosf),

dx^{2}= dr sinq sinf + r d (sinq sinf),

and dx^{3}= dr cosq + r d (cosq).

(dx^{1})^{2}= (dr sinq cosf)^{2}+ r^{2}[d (sinq cosf)]^{2}+ 2rdr sinq cosf (cosq dq cosf - sinq sinf df),

(dx^{2})^{2}= (dr sinq sinf)^{2}+ r^{2}[d (sinq sinf)]^{2}+ 2rdr sinq sinf (cosq dq sinf + sinq cosf df),

and (dx^{3})^{2}= (dr cosq)^{2}+ r^{2}[d (cosq)]^{2}- 2rdr sinq cosq dq.

- (d cos a/du)
^{2}+ (d cos b/du)^{2}+ (d cos g/du)^{2}= (1+p^{2}+q^{2})^{ -2 }T^{ -4}[For, p.17, l.10].

Hint. T'^{2}= (p'p"+q'q"+p^{2}q'q"-pqq'p"-pqp'q"+q^{2}p'p")^{2}/[p'^{2}+q'^{2}+(pq'-qp')^{2}].

TT' = pp"+qq"+p^{2}q'q"-pqq'p"-pqp'q"+q^{2}p'p" (get rid of the denominator of T' or consider TT' the derivative of T^{2}/2).

T^{4}[(d cos a/du)^{2}+ (d cos b/du)^{2}+ (d cos g/du)^{2}]

= (Tq"-q'T')^{2}+ (-Tp"+p'T')^{2}+ (Tpq"-Tqp"-pq'T'+qp'T')^{2}

= (p'p"+q'q"+p^{2}q'q"-pqq'p"-pqp'q"+q^{2}p'p")^{2}(the sum of the coefficients of T'^{2}cancels the denominator of T'^{2})

+ T^{2}(p"^{2}+q"^{2}+p^{2}q"^{2}+q^{2}p"^{2}-2pqp"q")

- 2(TT')^{2}(recognize the pattern: the sum of the coefficients of TT' is -2TT').

- (Systematic reduction due to the intrinsic structure)

Prove ds ds sin c = V(dpdq - dpdq) [For, p.35, l.-9].

Proof. cos^{2}c +sin^{2}c

= E[E(dp/ds)^{2}+2F(dp/ds)(dq/ds)+G(dq/ds)^{2}](dp/ds)^{2}

+ 2F[E(dp/ds)^{2}+2F(dp/ds)(dq/ds)+G(dq/ds)^{2}](dp/ds)(dq/ds)

+ G[E(dp/ds)^{2}+2F(dp/ds)(dq/ds)+G(dq/ds)^{2}](dq/ds)^{2}

= (ds^{2}/ds^{2})[E(dp/ds)^{2}+2F(dp/ds)(dq/ds)+G(dq/ds)^{2}] = 1.

- High-level cancellations versus low-level cancellations.

Low-level cancellations are term by term cancellations. High-level cancellations are cancellations of groups of terms. It would be too much work to perform cancellations at the lowest level even though success would be guaranteed. If possible, we would like to use higher-level cancellations because these cancellations are more efficient.

Example 1. Prove V(¶^{2}V/¶ x_{1}^{2}) = GV^{2}[For, p.37, l.1].

Proof. V(¶^{2}V/¶ x_{1}^{2}) = y_{2}^{2}+z_{2}^{2}-y_{2}^{2}Z^{2}+2y_{2}z_{2}YZ-z_{2}^{2}Y^{2}.

GX^{2}= x_{2}^{2}+y_{2}^{2}+z_{2}^{2}-x_{2}^{2}Y^{2}-y_{2}^{2}Y^{2}-z_{2}^{2}Y^{2}-x_{2}^{2}Z^{2}-y_{2}^{2}Z^{2}-z_{2}^{2}Z^{2}.

Up to this stage, we retain Y and Z as they are. We need not use [For, p.36, l.15] to expand them. Due to these high-level cancellations, it is enough to prove x_{2}^{2}X^{2}= (y_{2}Y+z_{2}Z)^{2}. Only for this step is it necessary to expand X, Y and Z using [For, p.36, l.15].

Example 2. Prove V(¶^{2}V/¶ x_{1}¶ y_{1}) = GXY [For, p.37, l.4].

In this case, unlike Example 1, if we want to compare the terms on the two sides for possible cancellations, we must expand X, Y and Z using [For, p.36, l.15]. We have no way to retain high-level quantities X, Y and Z and enjoy cancellations at the same time.

- Using the hierarchy of mathematical structures as the guide to find the right direction for cancellations.

Example. Prove the coefficient of a_{1}is 0 [For, p.53, l.9].

A good start is half of the battle. If we try to perform cancellations using the high level relations in [For, p.48, l.-8-p.49. l.6], it will lead nowhere. Consequently, we follow the hierarchy to go to a lower level.

(¶/¶q)(s/V^{2}) = (1/V^{2})(¶s/¶q)-2(s/V^{3})(¶V/¶q)

= (F_{2}M+FM_{2}-G_{2}L-GL_{2})/V^{2}-2s(G'+D")/V^{2}.

Then we use [For, p.44, l.-1; p.45, l.8] to expand G_{2}and use [For, p.44, l.-4; p.45, l.7-l.8] to expand F_{2}.

- Do not open Pandora's box. Otherwise, you will unnecessarily divide a proof
into numerous cases.

Example. Prove x_{1}x_{2}+y_{1}y_{2}+z_{1}z_{2}= 0 [Wea1, vol. 1, p.125, (5); p.131, (16)].

Proof. The trick is to use the identity x_{1}x_{2}= [(xx_{1})(xx_{2})]/(x^{2}),

where xx_{1}and xx_{2}can be easily derived from [Wea1, vol. 1, p.130, (15)].

Remark. Do not try to find x_{1}from x_{1}^{2}because you will have a difficult time determining the sign of x_{1}. Even if you are able to figure out the sign of x_{1}case by case, you still unnecessarily complicate the proof by dividing it into numerous cases.

- Once you recognize a pattern, you will see the light at the end of the tunnel.

Example. Prove [Wea1, vol. 1, p.185, (7)].

Keep (b+b') as one unit. If it is a factor, do not try to expand the product by breaking the unit apart. After you substitute [Wea1, vol. 1, p.185, (5)] into the left-hand side of [Wea1, vol. 1, p.185, (7)], you try to find the numerator by calculating its coefficients of (du)^{4}, (dv)^{4}, (du)^{3}(dv), (du)(dv)^{3}, and (du)^{2}(dv)^{2}. You hope every coefficient is zero, but none of them are. The calculations for the first pair are easiest. The calculations for the second pair are more difficult. The calculations for the term (du)^{2}(dv)^{2}are most difficult. After you calculate the coefficients for the first pair, you will find their coefficients are the corresponding coefficients of

(-1/4)[Wea1, vol. 1, p.185, (6)]^{2}.

Then the rest calculations become routine and any error that you commit can be easily debugged.

- (Creating checkpoints) The proof of [Gon1, p.429, (5.17-20)] requires long
calculations. It is frustrating if one discovers that one's final answer does
not match the answer given in the book. One might miss a term; one might copy a
term incorrectly. One does not know where one has made a mistake. If one has
made a mistake, then all one's work beyond that point is simply a waste.
Consequently, it is important to create some checkpoints along the way to help
one detect one's mistakes at an early stage and to establish the reliability
of one's work by passing some tests at the checkpoint. If one
can pass these tests, then it is less likely that
one's work before the checkpoint contains any mistakes. For example, one
may want to prove the simpler formula [Gon1, p.429, (5.17-21)] first (the first
checkpoint) because its highest power is 6 instead of 7. If one can prove [Gon1,
p.429, (5.17-21], the proof at least ensures that one's Taylor series expansion of Tan
^{2}(z/2) is correct. Then one may use this correct expansion to prove the more complicated formula [Gon1, p.429, (5.17-20)]. After one expresses the coefficient f(l) of z^{7}in terms of l, one may want to check if f(1) = - (7!)^{-1}(the second checkpoint) before one actually substitutes l = 1 - 2k^{2}into f(l). When one substitutes l = 1 - 2k^{2}into f(l), one should check the coefficient of k^{4}first (the third checkpoint). This will give one another chance to ensure that the coefficients of the powers of l in f(l) are correct before one calculates the more complicated coefficient of k^{2}.

- (Reducing gamma functions to sine functions)

[Guo, p.160, (3)] can be proved as follows:

By [Guo, p.99, (2)], it suffices to prove

[sin (pa)]^{-1}[sin (pb)]^{-1}([sin (p[g-a])]^{-1}[sin (p[g-b])]^{-1}- [sin (pg)]^{-1}[sin (p[g-a-b])]^{-1})

= ([sin (pg)][sin (p[g-a])][sin (p[g-b])][sin (p[a+b-g])])^{-1}. - (Choice will become clear after cancellation)

Prove the equality given in [Wat1, p.291, l.-14-l.-13].

Proof. Let A satisfy

(d^{2}I/dz^{2})+{(1+a-a')(z-*a*)^{-1}+(1+b-b')(z-b)^{-1}+(1+g-g')(z-c)^{-1}}(dI/dz)+{A[(z-*a*)(z-b)(z-c)]^{-1}}I = 0.

It suffices to prove A = (a+b+g){(a+b+g+1)z + S*a*(a+b'+g'-1)}.

(1). u = (z-*a*)^{a}(z-b)^{b}(z-c)^{g}.

du/dz = [a(z-*a*)^{a-1}(z-b)^{b}(z-c)^{g}+b(z-*a*)^{a}(z-b)^{b-1}(z-c)^{g}+g(z-*a*)^{a}(z-b)^{b}(z-c)^{g-1}]I+(z-*a*)^{a}(z-b)^{b}(z-c)^{g}(dI/dz).

d^{2}I/dz^{2}= {a(a-1)(z-*a*)^{a-2}(z-b)^{b}(z-c)^{g}+b(b-1)(z-*a*)^{a}(z-b)^{b-2}(z-c)^{g}+g(g-1)(z-*a*)^{a}(z-b)^{b}(z-c)^{g-1}

+2[ab(z-*a*)^{a-1}(z-b)^{b-1}(z-c)^{g}+ag(z-*a*)^{a-1}(z-b)^{b}(z-c)^{g-1}+bg(z-*a*)^{a}(z-b)^{b-1}(z-c)^{g-1}]}I

+2[a(z-*a*)^{a-1}(z-b)^{b}(z-c)^{g}+b(z-*a*)^{a}(z-b)^{b-1}(z-c)^{g}+g(z-*a*)^{a}(z-b)^{b}(z-c)^{g-1}](dI/dz)

+(z-*a*)^{a}(z-b)^{b}(z-c)^{g}(d^{2}I/dz^{2}).

(2). The part of A that the I term of {(1-a-a')(z-a)^{-1}+(1-b-b')(z-b)^{-1}+(1-g-g')(z-c)^{-1}}du/dz contributes is

{(1-a-a')(z-*a*)^{-1}+(1-b-b')(z-b)^{-1}+(1-g-g')(z-c)^{-1}}{a(z-*a*)^{-1}+b(z-b)^{-1}+g(z-c)^{-1}}(z-*a*)(z-b)(z-c)

= (1-a-a')(z-*a*)^{-1}a(z-b)(z-c)+(1-b-b')(z-b)^{-1}b(z-*a*)(z-c)+(1-g-g')(z-c)^{-1}g(z-*a*)(z-b)

+(1-a-a')(bz-bc+gz-gb)+(1-b-b')(az-ac+gz-g*a*)+(1-g-g')(az-ab+bz-ba).

Let B = (1-a-a')(z-*a*)^{-1}a(z-b)(z-c)+(1-b-b')(z-b)^{-1}b(z-*a*)(z-c)+(1-g-g')(z-c)^{-1}g(z-*a*)(z-b)

and C = (1-a-a')(bz-bc+gz-gb)+(1-b-b')(az-ac+gz-g*a*)+(1-g-g')(az-ab+bz-ba).

(3). The part of A that the I term of d^{2}I/dz^{2}contributes is

[a(a-1)(z-*a*)^{-2}+b(b-1)(z-b)^{-2}+g(g-1)(z-c)^{-2}

+2{ab[(z-a)(z-b)]^{-1}+ag[(z-*a*)(z-c)]^{-1}+ag[(z-*a*)(z-c)]^{-1}}](z-*a*)(z-b)(z-c)

= a(a-1)(z-*a*)^{-1}(z-b)(z-c)+b(b-1)(z-b)^{-1}(z-*a*)(z-c)+g(g-1)(z-*c*)^{-1}(z-a)(z-b)

+[2abz-2abc+2agz-2agb+2bgz-2bg*a*].

Let D = a(a-1)(z-*a*)^{-1}(z-b)(z-c)+b(b-1)(z-b)^{-1}(z-*a*)(z-c)+g(g-1)(z-*c*)^{-1}(z-a)(z-b),

E = 2abz-2abc+2agz-2agb+2bgz-2bg*a*,

and F = aa'(*a*-b)(*a*-c)(z-*a*)^{-1}+bb'(b-c)(b-*a*)(z-b)^{-1}+gg'(c-*a*)(c-b)(z-c)^{-1}.

Then G=B+D+F=aa'(-z-*a*+b+c)+bb'(-z-b+*a*+c)+gg'(-z-c+*a*+b).

(4). J =*a*(a+b'+g'-1)(a+b+g)+ b(b+g'+a'-1)(a+b+g)+c(g+a'+b'-1)(a+b+g)

=*a*a(a+b'+g'-1+b+g)+bb(a+b+g'+a'-1+g)+cg(a+b+g+a'+b'-1)+H

= (-aa'*a*-bb'b-gg'c)+H.

A=E+G+C

=[-aa'-bb'-gg'-a'g-b'a-b'g-ag'-g'b+2a+2b+2g]z+(-aa'a-bb'b-gg'c)+H (after cancellation)

= (a+b+g)(2-a'-b'-g')z+(-aa'*a*-bb'b-gg'c)+H

= (a+b+g)(a+b+g+1)z+J.

Remark. Before cancellation, we do not know whether we should use 2-a'-b'-g' or a+b+g+1. After cancellation, the answer is clear.

- Suppose we want to use [Wat1, p.308, (II) and (III)] to prove [Wat1, p.308, (I)]. See [Wat1, p.309, l.20]. Compare terms in [Wat1, p.308, (III)] with those in [Wat1, p.308, (I)]. All the terms in [Wat1, p.308, (III)] appear in
[Wat1, p.308, (I)] except P'
_{n-1}(z). After differentiating [Wat1, p.308, (II)], we try to use the resulting equality and [Wat1, p.308, (III)] to cancel the term P'_{n-1}(z). - Prove the formula given in [Wat1, p.317, l.13].

Proof. I. Express the integrals given in [Wat1, p.317, l.10] in terms of the values of the gamma function [Wat1, p.253, l.-6; p.254, l.-9].

II. Find the power series expansion in z^{-2}for F(2^{-1}(n+1), (n/2)+1; n+(3/2); z^{-2}). Use [Guo, p.94, (6)] to express the coefficients of the series in terms of the values of the gamma function.

III. Use the duplication formula [Wat1, p.240, l.-3].

Remark 1. The gamma function is a useful tool to replace any quantity that contains many factors such as a factorial or a notation that can be expressed in terms of the gamma function [Guo, p.94, (6)].

Remark 2. In order to expand the domain of Q_{n}with respect to n, it is unnecessary to repeat the same procedure as in [Wat1, p.317, l.-12-l.-4] .
All we have to do is find its analytic continuation with respect to n using
[Wat1, p.288, l.3, l.6, & l.10]. - Prove the equality given in [Wat1, p.319, l.-14].

Proof. 1 - t^{2 }= 4e^{q}(z+1)^{1/2}(z-1)^{1/2}[e^{q}(z-1)^{1/2}+ (z+1)^{1/2}]^{-2}.

z - t = [ze^{q}(z+1)^{1/2}+z(z-1)^{1/2}- e^{q}(z+1)^{1/2}+(z-1)^{1/2}]/[e^{q}(z+1)^{1/2}(z-1)^{1/2}].

dt = 2e^{q}(z+1)^{1/2}(z-1)^{1/2}[e^{q}(z-1)^{1/2}+ (z+1)^{1/2}]^{-2}.

[e^{q}(z-1)^{1/2}+ (z+1)^{1/2}] [ze^{q}(z+1)^{1/2}+z(z-1)^{1/2}- e^{q}(z+1)^{1/2}+(z-1)^{1/2}]

= e^{2q}(z^{2}-1)+2e^{q}(z+1)^{1/2}(z-1)^{1/2}+ z^{2 }- 1

= (2e^{q}(z+1)^{1/2}(z-1)^{1/2})[(e^{q}+e^{-q})/2].

Remark. In order to fully utilize cancellation, we should not expand [e^{q}(z-1)^{1/2}+ (z+1)^{1/2}]^{2}except when calculating 1 - t^{2}. When a textbook omits long calculations, it can blur the train of thought or the key to solution. - In order to fully take advantage of cancellation, we must group factors properly.

Example. Prove the formula given in [Wat1, p.329, l.6].

Proof. Let A = e^{i(f-w)}cosh (x/2) + sinh (x/2),

B = e^{i(f-w)}sinh (x/2) + cosh (x/2),

C = e^{if}cosh (x'/2) + sinh (x'/2),

D = cosh (x'/2)+ e^{if}sinh (x'/2),

E = sinh (x/2) sinh (x/2) - e^{iw}cosh (x/2) sinh (x'/2),

F = cosh (x/2) sinh (x/2) - e^{iw}sinh (x/2) sinh (x'/2), and

G = e^{iw}sinh x sinh^{2}(x'/2) + e^{-iw}sinh x cosh^{2}(x'/2) - cosh x sinh x'

Then dt/df = ie^{if}GD^{-2}; EF = 2^{-1}e^{iw}G;

(t-1)(t+1)(t-z)^{-1}= 2e^{iw}AB(CD)^{-1}

= 2[x+(x^{2}-1)^{1/2}cos (w-f)][x'+(x'^{2}-1)^{1/2}cos (w-f)]^{-1};

(t-z)^{-1}dt = [Di e^{if}G df][CGD^{2}]^{-1}

=[i e^{if}df][e^{if}(x'+(x'^{2}-1)^{1/2}cos f)]^{-1}.

Remark. Analytic continuation is a process of generalization. The process is ineffective: there are many ways to select the contour of the integral given in [Wat1, p.328, l.-8]; some of them are poor choices for calculation. Appropriate specification is the key to obtaining the most*effective*analytic continuation. We should parameterize t with a parameter f which can reduce calculation by a tremendous amount of cancellation. The choice given in [Wat1, p.328, l.-6] meets this requirement.¬ - If one can ensure that the essential formulas are correct, then the remaining task will be easy.

In proving the equality given in [Wat1, p.340, l.2-l.3], if one can obtain the formula for d^{2}v/dz^{2}correctly, the remaining task will be easy. It is difficult to make serious mistakes except for proving the following two formulas:

dv/dt =(k-1/2+m) ò -
If one tries to catch a snake, everything will become easy once one holds its neck. Similarly, if one tries to solve a problem, everything will become easy once one knows how to start.

Example. Prove [Guo, p.351, (5) & (6)].

Proof. Even if [Guo, p.351, l.3] provides a hint, it still requires some manipulation of [Guo, p.350, (3); p.351, (4)] to get start.

[d/(zdz)]^{n }{(sin z)/z} = (-1)^{n }(nz/2)^{1/2 }z^{-n-1 }J_{n+1/2}(z). ¬

Remark. We may prove [Guo, p.351, (6)] in a similar manner.

- Prove the equality given in [Wat1, p.410, l.-6-l.-3].

Proof. In this proof, all we need to do is find a general formula to express

cos^{2m}z as S_{n=0}^{m}a_{mn}cos 2nz.

We can write down the recursion formula for a_{mn}using mathematical induction: cos^{2m}z = [(cos 2z +1)/2] cos^{2m-2}z.

Remark. This proof requires the application of the above formulas many times. One should use the formulas rather than repeat the calculations to produce the same formulas many times. Avoiding unnecessary calculations greatly reduces the
trouble for checking errors in a mammoth calculation. One may prove the statements given in [Wat1, p.410, l.-2-l.-1] using §32.2 Fourier Expansions; Functions with q<0 of the following webpage: -
Prove the formula given in [Wat1, p.411, Example 1, (i)].

Proof. All one needs to do is generalize the method of proving the equality given in [Wat1, p.410, l.-6-l.-3]. There are many aspects that need to be taken care of.

I. cos^{2r}z = 2^{-2r }(e^{iz }+ e^{-iz})^{2r}= 2^{-2r}C(2r, r) + 2^{1-2r}S_{k=0}^{r-1}C(2r, k) cos [(2r - 2k) z].

II. The q^{r}-term of the coefficient of cos 2rz in the expansion of ce_{0}(z, q) is the coefficient of cos 2rz in b_{r}(z).

(2p)^{-1}ò_{[-p,p]}S_{n=0}^{r}32^{n}[(2n)!]^{-1}cos^{2n}z cos^{2n}q b_{r-n}(q) dq = S_{n=0}^{r}a_{n}b_{r-n}(z) (*).

Hence the coefficient of cos 2rz in b_{r}(z) is 2^{r+1}[r!r!]^{-1}.

III. Now consider the coefficient of the q^{r+2}-term of the coefficient of cos 2rz in the expansion of ce_{0}(z, q).

Since b_{r+1}cannot contain the term cos 2rz, it suffices to consider the coefficient of cos 2rz in b_{r+2}(z).

Replace r with r+2 in (*) and consider only the coefficient of cos 2rz in the expansion of each of the following terms:

(2p)^{-1}ò_{[-p,p]}[32^{r+2}[[2(r+2)]!]^{-1}cos^{2(r+2)}z cos^{2(r+2)}+ 32^{r+1}[[2(r+1)]!]^{-1}cos^{2(r+1)}z cos^{2(r+1)}q b_{1}(q)

+ 32^{r}[(2r)!]^{-1}cos^{2r}z cos^{2r}q b_{2}(q) +…] dq = …+ a_{2}b_{r}(z) + … + b_{r+2}(z).¬ - (Double indices)

Suppose we count summands row by row first, and then count them column by column, and we want to prove the two resulting sums are equal. It would be difficult if one tries to manage the double indices. It is better to prove the equality for a simple case; the solution will become clear as one proceeds.

Example 1. (The resulting sums are the same no matter whether we count summands row by row or column by column) Prove the statement given in [Cod, p.146, l.-8-l.-7].

Proof. I. Case r = 3.

(P_{1}^{~}Q_{0}-A_{0}P_{1}^{~})D_{2}= (A_{1}P_{0}^{~}-P_{0}^{~}Q_{1})D_{2}(1)

(P_{2}^{~}Q_{0}-A_{0}P_{2}^{~})D_{1}= [(A_{1}P_{1}^{~}-P_{1}^{~}Q_{1})+(A_{2}P_{0}^{~}-P_{0}^{~}Q_{2})]D_{1}(2)

(P_{3}^{~}Q_{0}-A_{0}P_{3}^{~})D_{0}= [(A_{1}P_{2}^{~}-P_{2}^{~}Q_{1})+(A_{2}P_{1}^{~}-P_{1}^{~}Q_{2})+(A_{3}P_{0}^{~}-P_{0}^{~}Q_{3})]D_{0}(3)

The left-hand side of [Cod, p.144, (2.7(ii)] counts summands row by row, while the right-hand of [Cod, p.144, (2.7(ii)] counts summands column by column.

II. P_{k}Q_{0}-A_{0}P_{k}

= S_{i=0}^{k}(P_{i}^{~}A_{0}-A_{0}P_{i}^{~})D_{k-i}

= S_{i=0}^{k}S_{l=0}^{i}(A_{k-i}P_{i-l}^{~}-P_{i-l}^{~}Q_{k-i})D_{l}(The sum by adding columns = The sum by adding rows)

= S_{i=0}^{k}(A_{k-i}P_{i}-P_{i}Q_{k-i}).

Example 2. (The number of double indices are the same no matter whether we count the indices row by row or column by column) Prove [Cod, p.146, (2.17)].

Proof. By counting the lattice points row by row or column by column, we have the same index set. Consequently,

S_{l=1}^{r}S_{k=1}^{r+1-l}= S_{k=1}^{r}S_{l=1}^{r+1-k}. - Prove that y = (cf
_{1}+f_{2})(cf_{3}+f_{4})^{-1}satisfies a Riccati equation [Inc1, p.24, l.-15-l.-13].

Proof. Eliminate C between the expression for y and the derived expression for y'. It suffices to note that

C = (yf_{4}-f_{2})(f_{1}-yf_{3})^{-1}and Cf_{3}+f_{4 }= (f_{4}f_{1}-f_{2}f_{3})(f_{1}-yf_{3})^{-1}. The numerator of Cf_{3}+f_{4 }does not contain y. - The general solution of z
^{2}(d^{2}u/d z^{2})+(2a-2bn-1)z(du/dz)+{b^{2}g^{2}z^{2b}+a(a-2bn)}u = 0 is u = z^{bn-a}C_{n}(gz^{b}) [Wat, p.97, l.11-l.14].

Remark. Let z = gx^{b}and u = yx^{-a}. Then dz = gbx^{b-1}dx.

du/dx=(du/dz)(dz/dx).

When calculating du/dz, we may treat x^{-a}as a function of z.

du/dz = x^{-a}(dy/dz)-ax^{-a-1}(dx/dz)y

= x^{-a}(dy/dz)-ax^{-a-1}[(gb)^{-1}x^{1-b}]y.

We may also express x^{-a}in terms of z first.

du/dz=(dy/dz)(z/r)^{-a/b}+y(-a/b)(z/g)^{-a/b-1}g^{-1};

then express z in terms of x. The calculation in the former approach is simpler. - u = z
^{2bn-a}is a solution of z^{2}(d^{2}u/d z^{2})+(2a-2bn-1)z(du/dz)+a(a-2bn)u = 0 [Wat, p.97, l.15-l.16].

Proof. du/dz = (2bn-a)z^{2bn-a-1}and d^{2}u/d z^{2}= (2bn-a)(2bn-a-1)z^{2bn-a-2}.

z^{2}(d^{2}u/d z^{2})+(2a-2bn-1)z(du/dz)+a(a-2bn)u

= [(2bn-a)(2bn-a-1)+(2a-2bn+1)(2bn-a)-a(2bn-a)]z^{2bn-a}

= [(2bn-a)^{2}-(2bn-a)-(2bn-a)+a+1-a(2bn-a)]z^{2bn-a}

= (2bn-a)[(2bn-a)-1-(2bn-a)+a+1-a]=0.

Remark. If one were to expand the brackets before cancelling like terms, it would take more calculations to obtain the same result. - d
^{2}{y/c(z)}/d{y(z)}^{2 }- (2n-1)[y(z)]^{-1}(d{y/c(z)}/dy(z))+y/c(z) = 0 can be reduced to

d^{2}y/dz^{2}-{y"(z)/y'(z)+(2n-1)y'(z)/y(z)+2c'(z)/c(z)}dy/dz

+[{y"(z)/y'(z)+(2n-1)y'(z)/y(z)+2c'(z)/c(z)}c'(z)/c(z) - c"(z)/c(z)+{y'(z)}^{2}]y = 0 [Wat, p.97, l.7-l.9].

Proof. dy(z)/dz = y'(z).

d{y/c(z)}/dz = [c(z)dy/dz-yy'(z)]/c^{2}(z).

d{y/c(z)}/dy(z) = [c(z)dy/dz-yy'(z)]/[y'(z)c^{2}(z)].

It suffices to note that d[d{y/c(z)}/dy(z)]/dz = y'(z)[d^{2}{y/c(z)}/d{y(z)}^{2}]. - Let s = 2im[zw]
^{1/2}/t and f(w) is an arbitrary function of w. Then

u = t^{-1}z^{-n/2}ò_{0}^{¥}exp [-m(z+w)/t]C_{n}(s)f(w)dw is a solution of

z(¶^{2}u/¶z^{2})+(1+n)(¶u/¶z)-m(¶u/¶t) = 0 [Wat, p.99, l.8-l.14].

Proof. z(¶^{2}u/¶z^{2})+(1+n)(¶u/¶z)-m(¶u/¶t)

= (4zt)^{-1}z^{-n/2}ò_{0}^{¥}exp [-m(z+w)/t]{s^{2}(¶^{2}C_{n}(s)/¶s^{2})+s(¶C_{n}(s)/¶s)+(s^{2}-n^{2})C_{n}(s)}f(w)dw = 0. - S(bn-gm)(a
^{2}-l)^{-1}(b^{2}-l)^{-1 }= l^{2}(a^{2}-l)^{-1}+m^{2}(b^{2}-l)^{-1}+n^{2}(c^{2}-l)^{-1}[Bel63, p.180, l.16].

Proof. By the equation given in [Bel63, p.180, l.14-l.15],

(l^{2}(a^{2}-l)^{-1}+m^{2}(b^{2}-l)^{-1}+n^{2}(c^{2}-l)^{-1})(a^{2}(a^{2}-l)^{-1}+b^{2}(b^{2}-l)^{-1}+g^{2}(c^{2}-l)^{-1})-(l^{2}(a^{2}-l)^{-1}+m^{2}(b^{2}-l)^{-1}+n^{2}(c^{2}-l)^{-1})

=((al(a^{2}-l)^{-1}+bm(b^{2}-l)^{-1}+)^{2}+gn(c^{2}-l)^{-1})^{2}.

If we subtract the term on the right-hand side of the above equation from the first term on the left-hand side, we obtain the term on the left-hand side of the desired equation. - Au
^{2}+Bv^{2}+Cw^{2}+2Fvw+2Gwu+2Huv-dD = -S [Bel63, p.217, l.-2].

Proof. It is unnecessary to lower the order of the determinant by one at a time. We expand S along the 4th column. In order to find the cofactor of u, we expand the cofactor along the third row. By inspection, we obtain -(uA+vH+wG). - Let A=(x
^{2}+y^{2}+z^{2})(*a*^{2}x^{2}+b^{2}y^{2}+c^{2}z^{2})-*a*^{2}(b^{2}+c^{2})x^{2}-c^{2}(*a*^{2}+b^{2})z^{2}+*a*^{2}b^{2}c^{2},

B=b^{2}(r^{2}-*a*^{2})(r^{2}-c^{2})+(*a*^{2}-b^{2})(r^{2}-c^{2})x^{2}-(b^{2}-c^{2})(r^{2}-*a*^{2})z^{2}, and

C={b(r^{2}-*a*^{2})+xxc^{-1}(*a*^{2}-c^{2})}{b(r^{2}-c^{2})+zz*a*^{-1}(*a*^{2}-c^{2})}+(*a*^{2}-c^{2})(*a*c)^{-1}{(r^{2}-c^{2})xxc^{-1}-(r^{2}-*a*^{2})zz*a*^{-1}}{*a*xx+czz-*a*bc}. Then

A=B=C [Bel63, p.269, l.-16-l.-12].

Proof. I. A=B.

Proof. The constant term for A and B:*a*^{2}b^{2}c^{2}.

The terms of the second degree for B: -*a*^{2}b^{2}r^{2}-b^{2}c^{2}r^{2}+(*a*^{2}-b^{2})(-c^{2})x^{2}-(b^{2}-c^{2})(-*a*^{2})z^{2}.

The terms of degrees higher than 2 for B: r^{4}b^{2}+(*a*^{2}-b^{2})r^{2}x^{2}-(b^{2}-c^{2})r^{2}z^{2}=r^{2}[r^{2}b^{2}+(*a*^{2}-b^{2})x^{2}-(b^{2}-c^{2})z^{2}].

II. B=C.

Proof. Let D=b^{2}(r^{2}-*a*^{2})(r^{2}-c^{2}). It suffices to prove B-D=C-D.

C-D=b(r^{2}-c^{2})(*a*^{2}-c^{2})xxc^{-1}-b(r^{2}-*a*^{2})(*a*^{2}-c^{2})zz*a*^{-1}-xzxz(*a*c)^{-1}(*a*^{2}-c^{2})^{2}

+(*a*^{2}-c^{2})(*a*c)^{-1}{(r^{2}-c^{2})xxc^{-1}-(r^{2}-*a*^{2})zz*a*^{-1}}{-*a*bc}+(*a*^{2}-c^{2})(*a*c)^{-1}{(r^{2}-c^{2})xxc^{-1}-(r^{2}-*a*^{2})zz*a*^{-1}}{*a*xx+czz}

=-xzxz(*a*c)^{-1}(*a*^{2}-c^{2})^{2}+(*a*^{2}-c^{2})(*a*c)^{-1}{(r^{2}-c^{2})xxc^{-1}-(r^{2}-*a*^{2})zz*a*^{-1}}{*a*xx+czz} (by cancellation)

=(*a*^{2}-c^{2})c^{-2}x^{2}(r^{2}-c^{2})x^{2}-(*a*^{2}-c^{2})*a*^{-2}z^{2}(r^{2}-*a*^{2})z^{2}.

The result follows from [Bel63, p.268, l.-S3]. - Given j
_{k}(k=0,…,n-2), find j_{n-1}using S_{k=0}^{n-1}[k!(n-k)!]^{-1}j_{k}[Guo, p.2, (6)].

Solution. We may reduce the calculations by the following method:

j_{n-1}=-S_{k=0}^{n-2}[C_{k}^{n}/n]j_{k}=-S_{k=0}^{n-2}[(n-1)(n-2)…(n-(k-1))][k(k-1)…2]^{-1}]j_{k}. - Links {1, 2}.

d

Remark. Watson leaves out a factor, (-1)

http://www.elsevierdirect.com/companions/9780123846549/Chap_Mathieu.pdf

In the beginning of this webpage discusses the origins and applications of Mathieu's equation, but they are not the main advantages of this webpage. Its main advantage consists in providing various ideas for solving Mathieu's equation. Without doing exercises, one has no motive to study these new things and cannot sharpen one's problem solving abilities. Mathieu had a deep understanding in theoretic physics. Between 1873 and 1890, he wrote six volumes of mathematical physics (differential equations, capillarity, electrostatics, magnetostatics, electrodynamics and elasticity). He solved differential equations from the viewpoint of physics. Before quantum mechanics was born from laboratory experiments, his mathematical work had already laid the important part of its theoretic foundation. [Coh, vol. 2, chap. XI, Stationary Perturbation Theory] is basically an abstraction and a simplification of the solution to Mathieu's equation. Mathieu is less well-known in the field of quantum mechanics probably because his mathematical work is too difficult to understand for most physicists.