while studying a physics problem I found that asymptotically the incomplete elliptic integral of the third kind, (using the Mathematica conventions, where it is called EllipticPi),
$\Pi(n;\phi|m)=\int_0^{\sin\phi} \frac{dt}{(1-nt^2)\sqrt{(1-mt^2)(1-t^2)}},$
which seems to logarithmically diverge when $\sin^2 \phi$ approaches $1/n$, approaches the logarithm indeed in the following way:
$\lim_{\epsilon \rightarrow 0} \left[ 2(n-1)\Pi\left(n;\arcsin \sqrt{\frac{1-\epsilon}{n}}|(2-n)n\right) +\log\epsilon\right]=C(n)$
where $C(n)$ is a constant. This identity (could it be a new discovery?) I checked numerically and it might be independent of the value of $m$.
I need to know the constant $C(n)$ in terms of some known functions (or constants from integrations independent of $n$) if possible, for the value of $m$ given.
Mathematica will not help. Is there any chance to do it? I have no training in this apart from undergraduate analysis and use elliptic integrals for the first time here.
(by the way, for large $n$, I can from looking at the graph see $C(n)$ being of the form $a+log(n+b)$, but for $n$ around unity this is no exact fit while being similar (divergence to minus infinity for small $n$ included).
Recent Edit:
putting the elliptic integral and the logarithm in the same integral and doing what was suggested by user8268, splitting them into two nondiverging integrals, I have now problems with the more complicated, see my comment to the answer. Can anybody tell me if this is an elliptic integral of the third kind or if Mathematica is right when it gives me a mixture of them (including ugly roots, making evaluation of the limit impossible) and if so must there be a calculation mistake or might the advice given even not be correct?
notsoimportantdetails(
Background:
This comes from wanting to evaluate an improper integral (§) that I physically know and numerically see to converge. In turn this integral comes frome the difference of two diverging integrals. Mathematica cannot do the improper integral of the combination but it can do the indefinite integral. But then it splits up the integrand again and the result is the elliptic integral and the logarithm, each diverging when doing the limit. Mathematica cannot do the limit and as I said astonishingly does not even know that there is a divergence for the elliptic integral. When the elliptic integral and the logarithm are put together into one integral (use $-log\epsilon=\int_0^{\sqrt{\frac{1-\epsilon}{n}}}\frac{2tn}{1-t^2n}dt$ or transform the elliptic integral to obtain (o) below) we are more or less back where we started, I suppose - well, my version of Mathematica even hangs up for the indefinite evaluation now.
Edit: Since maybe all I have done so far is making it more complicated, the original problem is to evaluate the integral (§) of
$\frac{1-\sqrt{1-\frac{\text{r0} (-2 m+\text{r0})}{x (-2 m+x)}}}{\left(1-\frac{2 m}{x}\right) \sqrt{1-\frac{\text{r0} (-2 m+\text{r0})}{x (-2 m+x)}}}$
from $r0$ to infinity, where $r0>2m$. This has the obvious split into the two divergent integrals.
For your reference, the combined integral (o) I talked about in the end above reads
$C(n)=\lim_{\epsilon\rightarrow 0}\int_1^\epsilon\frac{1 - n + \sqrt{(1 + n (-1 + t) - 2 t) (-1 + t) (-1 + n +
t)}}{t \sqrt{(1 + n (-1 + t) - 2 t) (-1 + t) (-1 + n + t)}}$ where the integration will solve the problem by providing $C(n)$ ; $\epsilon$ can actually be replaced by $0$ right away since the limit exists.
notsoimportantdetails)
Answer
I finally after all this time found the answer by simply using the second identity on http://functions.wolfram.com/EllipticIntegrals/EllipticPi3/17/01/ which generally expresses the incomplete elliptic integral of the 3rd kind with one term of an incomplete integral of the 3rd kind, one term of an incomplete integral of the 1st kind and one logarithm, which cancels the logarithm in my limit.
The result has whence indeed the form suggested by user8268.
No comments:
Post a Comment