On the asymptotic behavior of the Fourier coefficients of $(1-R\cos\theta)^{-3/2}$ Today in class I showed some ways for dealing with the classical integral $\int_{0}^{2\pi}\frac{d\theta}{(A+B\cos\theta)^2}$ under the constraints $A>B>0$, including


*

*Symmetry and the tangent half-angle substitution;

*Relating the integral to the area enclosed by an ellipse, via the polar form of an ellipse with respect to a focus and the formula $\text{Area}=\pi a b =\frac{1}{2}\int_{0}^{2\pi}\rho(\theta)^2\,d\theta$; 

*Computing the geometric-like series $\sum_{n\geq 0} r^n \sin(n\theta)$ and applying Parseval's identity to it;

*Applying Feynman's trick (differentiation under the integral sign) to $\int_{0}^{2\pi}\frac{d\theta}{1-R\cos\theta}$ which is an elementary integral due to point 1.


I finished the lesson by remarking that the point $4.$ allows to compute $\int_{0}^{2\pi}\frac{d\theta}{\left(1-R\cos\theta\right)^3}$ almost without extra efforts, while the $L^2$ machinery (point 3.) does not seem to grant the same. Apparently I dug my own grave, since someone readily asked (with the assumption $R\in(-1,1)$)

What is the asymptotic behavior of the coefficients $c_n$ in
  $$ \frac{1}{\left(1-R\cos\theta\right)^{3/2}}= c_0+\sum_{n\geq 1}c_n \cos(n\theta) $$
  ? (Q2) Do we get something interesting by following the "unnatural" approach of applying Parseval's identity to such Fourier (cosine) series?

At the moment I replied that the Paley-Wiener theorem ensures an exponential decay of the $c_n$s, and with just a maybe to the second question. Later I figured out a pretty technical proof (through hypergeometric functions) of 
$$ |c_n| \sim K_R\cdot\sqrt{n}\cdot\left(\frac{|R|}{1+\sqrt{1-R^2}}\right)^n \quad \text{as }n\to +\infty$$
and the fact that $c_n$ is given by a linear combination of complete elliptic integrals of the first and second kind. (Q1) I would like to know if there is a more elementary way for deriving the previous asymptotic behaviour. And the outcome of (Q2).
 A: A sketch how to extract the wanted leading order asymptotics without knowledge of special functions:
The $n$-th Fourier coefficent of the function ${f(x,R)=(1-R\cos(x))^{-3/2}}$ is given by 
$$c_n(a)=\frac{ a^{3/2}}{\pi}\text{Re}\int_0^{2 \pi}\frac{e^{ixn}}{(a-\cos(x))^{3/2}}=-\frac{2a^{3/2}}{\pi}\text{Re}\frac{d}{da}\underbrace{\int_0^{2\pi}e^{ixn}F(x,a)dx}_{C_n(a)}
$$
here $a=1/R>1$ and  $(a-\cos(x))^{-1/2}=F(x,a)$
To get the asypmtotics of interest we integrate $e^{izn}F(z,a)$ over a rectangle $\mathcal{D}$ with verticies $[\epsilon,2\pi-\epsilon,2\pi -\epsilon+i\infty,\epsilon+i\infty]$ where $\epsilon\rightarrow0_+$.
Leaving the details to the interested reader, we notice that:


*

*the bottom part of the contour integral yields $C_n(a)$

*the top part vanishs due to the exponential decay at complex infinity

*the two vertical parts cancel out until we reach the branch point of $F(z,a)$ at $z_{0,1}=i\text{arccosh}(a)+2\pi \delta_{n,1}$

*there are no singularities contained inside $\mathcal{D}$


Taking these four points into account we get (the remaining parts of the vertical contours are actually two times the left piece due to a factor of $(-1)(i)^2$ stemming from the paramtrisation $z= i y$ and the phaseshift due to the branchcut
) :
$$
\oint_\mathcal{D}e^{izn}F(z,a)dz=0=C_n(a)-2\int_{\text{arccosh(a)}}^{\infty}\frac{e^{-ny}}{\sqrt{\cosh(y)-a}}dy
$$
the remaining integral is now pretty simple to estimate in the limit $n\rightarrow \infty$ since it is clearly dominated by the contribution of the branch point due to exponenital decay of the numerator. 
$$
C_n(a)\sim\frac{2}{(a^2-1)^{1/4}}\int_{\text{arccosh(a)}}^{\infty}dy\frac{e^{-ny}}{\sqrt{y-\text{arccosh(a)}}}=\\
\frac{4}{(a^2-1)^{1/4}}e^{-n\text{arccosh(a)}}\int_0^{\infty}dye^{-n y^2}=\\
\frac{2\sqrt{\pi}}{(a^2-1)^{1/4}}\frac{e^{-n\text{arccosh(a)}}}{\sqrt{n}}
$$
Taking the derivative w.r.t. $a$ gives us a mulitplicative factor of $n$ and OPs result follows by using the the connection between inverse hyperbolics and logarithms.

A similiar technique was used in the answer to this question
A: I don't think this is more elementary, but it may be a little more straightforward.  Start with Heine's toroidal identity,
$$ \frac{1}{\sqrt{z-\cos\,t}} = \frac{\sqrt{2}}{\pi} \sum_{m=-\infty}^\infty\,Q_{m-1/2}(z)\,\exp{(i \, m\, t}), \, \, |z|\ge1,  \, 0\le t \le 2\pi \, .$$
Splitting the summation into positive and negative $m$ and using the fact that for the toroidal function (Gradshteyn & Rhyzhik 8.737.4) 
$Q_{-(m-1/2)-1}(z) = Q_{m-1/2}(z)$ for integer $m$ then one obtains
$$ \frac{1}{\sqrt{z-\cos\,t}} = \frac{\sqrt{2}}{\pi}\Big(2 \sum_{m=1}^\infty\,Q_{m-1/2}(z)\,\cos{( m\, t)}+  Q_{-1/2}(z)\Big) .$$
Differentiate both sides with respect to $z$ and use G&R 8.706.1 (with the fact that $z>1$ to get the phase correct) to get
$$ (z-\cos\,t)^{-3/2} = \frac{-2\sqrt{2}}{\pi\,\sqrt{z^2-1}}\Big(2 \sum_{m=1}^\infty\,Q_{m-1/2}^{\,1}(z)\,\cos{( m\, t)}+  Q_{-1/2}^{\,1}(z)\Big) .$$
From the problem's statement, with $z=1/R,$ we see that the associated toroidal function $Q_{m-1/2}^{\,1}(z)$ must be asymptotically estimated.  G&R 8.777.2 has a hypergeometric representation that has an asymptotic nature for large $m,$
$$Q_{m-1/2}^{\,1}(z)=-2\sqrt{\pi}\,\frac{\Gamma(m+3/2)}{\Gamma(m+1)} 
\,\frac{\sqrt{z^2-1}}{\zeta^{\,m+3/2}} F\big(3/2,m+3/2; m+1; 1/\zeta^2\big)$$
where $\zeta = z+\sqrt{z^2-1} > 1.$  For convenience set $y=1/\zeta^2$ and expand the ration of Pochhamer symbols in the hypergeometric function as
$$\frac{(m+3/2)_k}{(m+1)_k} \sim 1 + \frac{k}{2m} + \, ...  \,\, m\to \infty.$$ Therefore
$$F\big(3/2,m+3/2; m+1; 1/\zeta^2\big)\sim \sum_{k=0}^\infty \dfrac{(3/2)_k}{k!} \big(1+\frac{k}{2m}\big)y^k= (1-y)^{-3/2}\Big(1+\frac{3y}{4m(1-y)} \Big).
$$
For this to be a useful asymptotic expression, $y$ needs to bounded away from 1, which implies $z$ must be bounded from 1, and $R$ must be bounded from 1.
Writing 
$$ (z-\cos\,t)^{-3/2} =d_0(z)+\sum_{m=1}^\infty d_m(z)\,\cos{(m\,t)} $$
then, using only the first term of the expression above and the first in the asymptotic expansion of the ratio of gamma functions, 
$$ d_m(z) \sim \frac{8\sqrt{2\, m}}{\sqrt{\pi}} (z+\sqrt{z^2-1})^{-3/2}(1-y)^{3/2}$$
In the proposer's notation it is clear to see that
$$c_m = R^{-3/2}d_m(z \to 1/R)$$ and, with some simplification one finds
$$c_m \sim \frac{4}{\sqrt{\pi}} \,\sqrt{m}\, (1-R^2)^{-3/4}
\Big(\frac{R}{1+\sqrt{1-R^2}}\Big)^m .$$
This expression was derived for $z>1$ implying $0\le R<1.$  To extend to
 $-1\le R \le 0,$ simply let $t \to t+\pi$ and a factor of $(-1)^m$ is introduced into the summand.  It is then clear that the $(-1)^m$ can be incorporated into the final expression as the proposer did, by making the unsigned $R$ an absolute value.
I have checked many of the intermediate steps in Mathematica, but to do so one must use $Q_{m-1/2}^{\,1}(z)=$ LegendreQ[m-1/2, 1, 3, z].
