Mixing
Inverse of elliptic integral of second kind
The Wikipedia articles on elliptic integral and elliptic functions state that “elliptic functions were discovered as inverse functions of elliptic integrals.” Some elliptic functions have names and are thus well-known special functions, and the same holds for some elliptic integrals. But what is the relation between the named elliptic functions and the named elliptic integrals?
It seems that the Jacobi amplitude $varphi=operatorname{am}(u,k)$ is the inverse of the elliptic integral of the first kind, $u=F(varphi,k)$. Or related to this, $x=operatorname{sn}(u,k)$ is the inverse of $u=F(x;k)$. It looks to me as if all of Jacobi's elliptic functions relate to the elliptic integral of the first kind. For other named elliptic functions listed by Wikipedia, like Jacobi's $vartheta$ function or Weierstrass's $wp$ function, it is even harder to see a relation to Legendre's integrals.
Is there a way to express the inverse of $E$, the elliptic integral of the second kind, in terms of some named elliptic functions? I.e. given $E(varphi,k)=u$, can you write a closed form expression for $varphi$ in terms of $k$ and $u$ using well-known special functions and elementary arithmetic operations?
In this post the author uses the Mathematica function FindRoot to do this kind of inversion, but while reading that post, I couldn't help wondering whether there is an easier formulation. Even though the computation behind the scenes might in fact boil down to root-finding in any case, it feels like this task should be common enough that someone has come up with a name for the core of this computation.
special-functions inverse elliptic-integrals elliptic-functions
edited Apr 13 '17 at 12:21
Community♦
1
asked Jan 28 '15 at 13:31
MvG
30.7k449101
add a comment |
The Wikipedia articles on elliptic integral and elliptic functions state that “elliptic functions were discovered as inverse functions of elliptic integrals.” Some elliptic functions have names and are thus well-known special functions, and the same holds for some elliptic integrals. But what is the relation between the named elliptic functions and the named elliptic integrals?
It seems that the Jacobi amplitude $varphi=operatorname{am}(u,k)$ is the inverse of the elliptic integral of the first kind, $u=F(varphi,k)$. Or related to this, $x=operatorname{sn}(u,k)$ is the inverse of $u=F(x;k)$. It looks to me as if all of Jacobi's elliptic functions relate to the elliptic integral of the first kind. For other named elliptic functions listed by Wikipedia, like Jacobi's $vartheta$ function or Weierstrass's $wp$ function, it is even harder to see a relation to Legendre's integrals.
Is there a way to express the inverse of $E$, the elliptic integral of the second kind, in terms of some named elliptic functions? I.e. given $E(varphi,k)=u$, can you write a closed form expression for $varphi$ in terms of $k$ and $u$ using well-known special functions and elementary arithmetic operations?
In this post the author uses the Mathematica function FindRoot to do this kind of inversion, but while reading that post, I couldn't help wondering whether there is an easier formulation. Even though the computation behind the scenes might in fact boil down to root-finding in any case, it feels like this task should be common enough that someone has come up with a name for the core of this computation.
special-functions inverse elliptic-integrals elliptic-functions
edited Apr 13 '17 at 12:21
Community♦
1
asked Jan 28 '15 at 13:31
MvG
30.7k449101
This paper might be interesting for you: Numerical computation of incomplete elliptic integrals of a general form by T. Fukushima and H. Ishizaki. Also the paper Numerical computation of inverse complete elliptic integrals of first and second kinds by T Fukushima.
– Loreno Heer
Jan 28 '15 at 16:07
yes inverse calculations are in the second paper
– Loreno Heer
Jan 28 '15 at 16:08
@sanjab: Thanks for these pointers. But as far as I can see, the first paper only does forward computation, and the second does inverse (using notation like $m_E$), but of the complete integral and solving for the parameter $m$ not the amplitude $varphi$. So at least as far as I can see, this is related but won't immediately help answering my question.
– MvG
Jan 28 '15 at 16:29
1
FWIW, Wolfram is looking for something nice, too.
– J. M. is not a mathematician
May 17 '16 at 3:40
N.B. the inverse of the Weierstrass elliptic function can of course be expressed in terms of $F$, or Carlson's $R_F$ if you prefer the symmetric integrals. Similarly, the Weierstrass and Jacobi functions are expressible in terms of each other (corresponding to the algebro-geometric notion of a quartic being equivalent to a cubic under suitable rational transformations.)
– J. M. is not a mathematician
May 17 '16 at 3:43
add a comment |
The Wikipedia articles on elliptic integral and elliptic functions state that “elliptic functions were discovered as inverse functions of elliptic integrals.” Some elliptic functions have names and are thus well-known special functions, and the same holds for some elliptic integrals. But what is the relation between the named elliptic functions and the named elliptic integrals?
It seems that the Jacobi amplitude $varphi=operatorname{am}(u,k)$ is the inverse of the elliptic integral of the first kind, $u=F(varphi,k)$. Or related to this, $x=operatorname{sn}(u,k)$ is the inverse of $u=F(x;k)$. It looks to me as if all of Jacobi's elliptic functions relate to the elliptic integral of the first kind. For other named elliptic functions listed by Wikipedia, like Jacobi's $vartheta$ function or Weierstrass's $wp$ function, it is even harder to see a relation to Legendre's integrals.
Is there a way to express the inverse of $E$, the elliptic integral of the second kind, in terms of some named elliptic functions? I.e. given $E(varphi,k)=u$, can you write a closed form expression for $varphi$ in terms of $k$ and $u$ using well-known special functions and elementary arithmetic operations?
In this post the author uses the Mathematica function FindRoot to do this kind of inversion, but while reading that post, I couldn't help wondering whether there is an easier formulation. Even though the computation behind the scenes might in fact boil down to root-finding in any case, it feels like this task should be common enough that someone has come up with a name for the core of this computation.
special-functions inverse elliptic-integrals elliptic-functions
edited Apr 13 '17 at 12:21
Community♦
1
asked Jan 28 '15 at 13:31
MvG
30.7k449101
The Wikipedia articles on elliptic integral and elliptic functions state that “elliptic functions were discovered as inverse functions of elliptic integrals.” Some elliptic functions have names and are thus well-known special functions, and the same holds for some elliptic integrals. But what is the relation between the named elliptic functions and the named elliptic integrals?
It seems that the Jacobi amplitude $varphi=operatorname{am}(u,k)$ is the inverse of the elliptic integral of the first kind, $u=F(varphi,k)$. Or related to this, $x=operatorname{sn}(u,k)$ is the inverse of $u=F(x;k)$. It looks to me as if all of Jacobi's elliptic functions relate to the elliptic integral of the first kind. For other named elliptic functions listed by Wikipedia, like Jacobi's $vartheta$ function or Weierstrass's $wp$ function, it is even harder to see a relation to Legendre's integrals.
Is there a way to express the inverse of $E$, the elliptic integral of the second kind, in terms of some named elliptic functions? I.e. given $E(varphi,k)=u$, can you write a closed form expression for $varphi$ in terms of $k$ and $u$ using well-known special functions and elementary arithmetic operations?
In this post the author uses the Mathematica function FindRoot to do this kind of inversion, but while reading that post, I couldn't help wondering whether there is an easier formulation. Even though the computation behind the scenes might in fact boil down to root-finding in any case, it feels like this task should be common enough that someone has come up with a name for the core of this computation.
special-functions inverse elliptic-integrals elliptic-functions
special-functions inverse elliptic-integrals elliptic-functions
edited Apr 13 '17 at 12:21
Community♦
1
asked Jan 28 '15 at 13:31
MvG
30.7k449101
edited Apr 13 '17 at 12:21
Community♦
1
asked Jan 28 '15 at 13:31
MvG
30.7k449101
edited Apr 13 '17 at 12:21
Community♦
1
edited Apr 13 '17 at 12:21
Community♦
1
edited Apr 13 '17 at 12:21
Community♦
1
1
asked Jan 28 '15 at 13:31
MvG
30.7k449101
asked Jan 28 '15 at 13:31
MvG
30.7k449101
asked Jan 28 '15 at 13:31
MvG
30.7k449101
30.7k449101
This paper might be interesting for you: Numerical computation of incomplete elliptic integrals of a general form by T. Fukushima and H. Ishizaki. Also the paper Numerical computation of inverse complete elliptic integrals of first and second kinds by T Fukushima.
– Loreno Heer
Jan 28 '15 at 16:07
yes inverse calculations are in the second paper
– Loreno Heer
Jan 28 '15 at 16:08
@sanjab: Thanks for these pointers. But as far as I can see, the first paper only does forward computation, and the second does inverse (using notation like $m_E$), but of the complete integral and solving for the parameter $m$ not the amplitude $varphi$. So at least as far as I can see, this is related but won't immediately help answering my question.
– MvG
Jan 28 '15 at 16:29
1
FWIW, Wolfram is looking for something nice, too.
– J. M. is not a mathematician
May 17 '16 at 3:40
N.B. the inverse of the Weierstrass elliptic function can of course be expressed in terms of $F$, or Carlson's $R_F$ if you prefer the symmetric integrals. Similarly, the Weierstrass and Jacobi functions are expressible in terms of each other (corresponding to the algebro-geometric notion of a quartic being equivalent to a cubic under suitable rational transformations.)
– J. M. is not a mathematician
May 17 '16 at 3:43
add a comment |
This paper might be interesting for you: Numerical computation of incomplete elliptic integrals of a general form by T. Fukushima and H. Ishizaki. Also the paper Numerical computation of inverse complete elliptic integrals of first and second kinds by T Fukushima.
– Loreno Heer
Jan 28 '15 at 16:07
yes inverse calculations are in the second paper
– Loreno Heer
Jan 28 '15 at 16:08
@sanjab: Thanks for these pointers. But as far as I can see, the first paper only does forward computation, and the second does inverse (using notation like $m_E$), but of the complete integral and solving for the parameter $m$ not the amplitude $varphi$. So at least as far as I can see, this is related but won't immediately help answering my question.
– MvG
Jan 28 '15 at 16:29
1
FWIW, Wolfram is looking for something nice, too.
– J. M. is not a mathematician
May 17 '16 at 3:40
N.B. the inverse of the Weierstrass elliptic function can of course be expressed in terms of $F$, or Carlson's $R_F$ if you prefer the symmetric integrals. Similarly, the Weierstrass and Jacobi functions are expressible in terms of each other (corresponding to the algebro-geometric notion of a quartic being equivalent to a cubic under suitable rational transformations.)
– J. M. is not a mathematician
May 17 '16 at 3:43
This paper might be interesting for you: Numerical computation of incomplete elliptic integrals of a general form by T. Fukushima and H. Ishizaki. Also the paper Numerical computation of inverse complete elliptic integrals of first and second kinds by T Fukushima.
– Loreno Heer
Jan 28 '15 at 16:07
This paper might be interesting for you: Numerical computation of incomplete elliptic integrals of a general form by T. Fukushima and H. Ishizaki. Also the paper Numerical computation of inverse complete elliptic integrals of first and second kinds by T Fukushima.
– Loreno Heer
Jan 28 '15 at 16:07
yes inverse calculations are in the second paper
– Loreno Heer
Jan 28 '15 at 16:08
yes inverse calculations are in the second paper
– Loreno Heer
Jan 28 '15 at 16:08
@sanjab: Thanks for these pointers. But as far as I can see, the first paper only does forward computation, and the second does inverse (using notation like $m_E$), but of the complete integral and solving for the parameter $m$ not the amplitude $varphi$. So at least as far as I can see, this is related but won't immediately help answering my question.
– MvG
Jan 28 '15 at 16:29
@sanjab: Thanks for these pointers. But as far as I can see, the first paper only does forward computation, and the second does inverse (using notation like $m_E$), but of the complete integral and solving for the parameter $m$ not the amplitude $varphi$. So at least as far as I can see, this is related but won't immediately help answering my question.
– MvG
Jan 28 '15 at 16:29
1
1
FWIW, Wolfram is looking for something nice, too.
– J. M. is not a mathematician
May 17 '16 at 3:40
FWIW, Wolfram is looking for something nice, too.
– J. M. is not a mathematician
May 17 '16 at 3:40
N.B. the inverse of the Weierstrass elliptic function can of course be expressed in terms of $F$, or Carlson's $R_F$ if you prefer the symmetric integrals. Similarly, the Weierstrass and Jacobi functions are expressible in terms of each other (corresponding to the algebro-geometric notion of a quartic being equivalent to a cubic under suitable rational transformations.)
– J. M. is not a mathematician
May 17 '16 at 3:43
N.B. the inverse of the Weierstrass elliptic function can of course be expressed in terms of $F$, or Carlson's $R_F$ if you prefer the symmetric integrals. Similarly, the Weierstrass and Jacobi functions are expressible in terms of each other (corresponding to the algebro-geometric notion of a quartic being equivalent to a cubic under suitable rational transformations.)
– J. M. is not a mathematician
May 17 '16 at 3:43
add a comment |
2 Answers
2
active
oldest
votes
I've just found a physical problem (in classical mechanics) involving the trajectory of a particle in which I had to take the inverse of $E(phi,k)$, the incomplete elliptic integral of the 2nd kind. For the elliptic integral of the 1st kind $F(x,k)$ this is an easy task beacuse the Jacobi elliptic function $sn(x,k)$ (or JacobiSN(x,k) in mathematical software) is just $F^{-1}(x,k)$. However, currently there is no built-in function for the inverse of $E(phi,k)$. My computational solution was then to build a procedure for the inverse, using FindRoot (in Mathematica9.0) or fsolve (in Maple 2015). F.M.S. Lima (University of Brasilia).
answered Feb 4 '16 at 21:54
Dr. Fabio M. S. Lima
765
add a comment |
I know this isn't a closed form, but I was interested in this question and have found one can relate the two functions together as series representations. I have written a short article here but this is the crux of it
Write begin{equation}
E(phi,k) = sum_{i=0}^infty frac{Q_i(k)}{(2i-1)!}phi^{2i-1}
end{equation}
where $Q_i(k)$ are polynomials in $k$, from the series expansion here we can get a finite form for these polynomials as begin{equation}
Q_n(m) = 2(-4)^nsum_{k=1}^n frac{(2k-3)!!}{k!}left(frac{-m}{8}right)^k sum_{j=0}^{k-1} binom{2k}{j}(-1)^{1-j}(j-k)^{2n}, ;;; n>0
end{equation}
with $Q_0(k)$ defined as $1$. Then you can write the inverse series using series reversion in a very similar manner to $E(phi,k)$ begin{equation}
phi(E,k) = sum_{i=0}^infty frac{R_i(k)}{(2i-1)!}E^{2i-1}
end{equation}
where the relation between the new polynomials $R_i(k)$ is given by the explicit reversion formula found at the bottom of this link, giving begin{equation}
R_n(k) = (2n)! sum_{tau_n=n}(-1)^{sigma_n} frac{prod_{i=1}^{sigma_n}2n+i}{prod_{j=1}^n k_j!}prod_{l=1}^n left(frac{Q_l(k)}{(2l+1)!}right)^{k_l}
end{equation}
where $sigma_n=k_1+k_2+k_3+cdots+k_n$, and $tau_n=k_1+2k_2+3k+3+cdots + nk_n$ and the sum is over all sets of indices $k_i$ that meet the requirement $tau_n=n$. I don't know if any nice simplifications or tricks can be made to reduce this to a functional form. The only numerical element here is converging the series to the desired accuracy.
answered May 25 '17 at 14:32
Benedict W. J. Irwin
1,463529
The incomplete elliptic integral of the second kind is known to be expressible in terms of a (bivariate) hypergeometric function, so the appearance of factorials and binomial coefficients is not that surprising. The inverse obtained from Lagrangian inversion, however, does not seem to show any known pattern in the coefficients.
– J. M. is not a mathematician
Jul 31 '17 at 6:10
add a comment |
Your Answer
StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");
StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});
function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f1123360%2finverse-of-elliptic-integral-of-second-kind%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
2 Answers
2
active
oldest
votes
2 Answers
2
active
oldest
votes
active
oldest
votes
active
oldest
votes
I've just found a physical problem (in classical mechanics) involving the trajectory of a particle in which I had to take the inverse of $E(phi,k)$, the incomplete elliptic integral of the 2nd kind. For the elliptic integral of the 1st kind $F(x,k)$ this is an easy task beacuse the Jacobi elliptic function $sn(x,k)$ (or JacobiSN(x,k) in mathematical software) is just $F^{-1}(x,k)$. However, currently there is no built-in function for the inverse of $E(phi,k)$. My computational solution was then to build a procedure for the inverse, using FindRoot (in Mathematica9.0) or fsolve (in Maple 2015). F.M.S. Lima (University of Brasilia).
answered Feb 4 '16 at 21:54
Dr. Fabio M. S. Lima
765
add a comment |
I've just found a physical problem (in classical mechanics) involving the trajectory of a particle in which I had to take the inverse of $E(phi,k)$, the incomplete elliptic integral of the 2nd kind. For the elliptic integral of the 1st kind $F(x,k)$ this is an easy task beacuse the Jacobi elliptic function $sn(x,k)$ (or JacobiSN(x,k) in mathematical software) is just $F^{-1}(x,k)$. However, currently there is no built-in function for the inverse of $E(phi,k)$. My computational solution was then to build a procedure for the inverse, using FindRoot (in Mathematica9.0) or fsolve (in Maple 2015). F.M.S. Lima (University of Brasilia).
answered Feb 4 '16 at 21:54
Dr. Fabio M. S. Lima
765
add a comment |
I've just found a physical problem (in classical mechanics) involving the trajectory of a particle in which I had to take the inverse of $E(phi,k)$, the incomplete elliptic integral of the 2nd kind. For the elliptic integral of the 1st kind $F(x,k)$ this is an easy task beacuse the Jacobi elliptic function $sn(x,k)$ (or JacobiSN(x,k) in mathematical software) is just $F^{-1}(x,k)$. However, currently there is no built-in function for the inverse of $E(phi,k)$. My computational solution was then to build a procedure for the inverse, using FindRoot (in Mathematica9.0) or fsolve (in Maple 2015). F.M.S. Lima (University of Brasilia).
answered Feb 4 '16 at 21:54
Dr. Fabio M. S. Lima
765
I've just found a physical problem (in classical mechanics) involving the trajectory of a particle in which I had to take the inverse of $E(phi,k)$, the incomplete elliptic integral of the 2nd kind. For the elliptic integral of the 1st kind $F(x,k)$ this is an easy task beacuse the Jacobi elliptic function $sn(x,k)$ (or JacobiSN(x,k) in mathematical software) is just $F^{-1}(x,k)$. However, currently there is no built-in function for the inverse of $E(phi,k)$. My computational solution was then to build a procedure for the inverse, using FindRoot (in Mathematica9.0) or fsolve (in Maple 2015). F.M.S. Lima (University of Brasilia).
answered Feb 4 '16 at 21:54
Dr. Fabio M. S. Lima
765
answered Feb 4 '16 at 21:54
Dr. Fabio M. S. Lima
765
answered Feb 4 '16 at 21:54
Dr. Fabio M. S. Lima
765
answered Feb 4 '16 at 21:54
Dr. Fabio M. S. Lima
765
765
add a comment |
add a comment |
I know this isn't a closed form, but I was interested in this question and have found one can relate the two functions together as series representations. I have written a short article here but this is the crux of it
Write begin{equation}
E(phi,k) = sum_{i=0}^infty frac{Q_i(k)}{(2i-1)!}phi^{2i-1}
end{equation}
where $Q_i(k)$ are polynomials in $k$, from the series expansion here we can get a finite form for these polynomials as begin{equation}
Q_n(m) = 2(-4)^nsum_{k=1}^n frac{(2k-3)!!}{k!}left(frac{-m}{8}right)^k sum_{j=0}^{k-1} binom{2k}{j}(-1)^{1-j}(j-k)^{2n}, ;;; n>0
end{equation}
with $Q_0(k)$ defined as $1$. Then you can write the inverse series using series reversion in a very similar manner to $E(phi,k)$ begin{equation}
phi(E,k) = sum_{i=0}^infty frac{R_i(k)}{(2i-1)!}E^{2i-1}
end{equation}
where the relation between the new polynomials $R_i(k)$ is given by the explicit reversion formula found at the bottom of this link, giving begin{equation}
R_n(k) = (2n)! sum_{tau_n=n}(-1)^{sigma_n} frac{prod_{i=1}^{sigma_n}2n+i}{prod_{j=1}^n k_j!}prod_{l=1}^n left(frac{Q_l(k)}{(2l+1)!}right)^{k_l}
end{equation}
where $sigma_n=k_1+k_2+k_3+cdots+k_n$, and $tau_n=k_1+2k_2+3k+3+cdots + nk_n$ and the sum is over all sets of indices $k_i$ that meet the requirement $tau_n=n$. I don't know if any nice simplifications or tricks can be made to reduce this to a functional form. The only numerical element here is converging the series to the desired accuracy.
answered May 25 '17 at 14:32
Benedict W. J. Irwin
1,463529
The incomplete elliptic integral of the second kind is known to be expressible in terms of a (bivariate) hypergeometric function, so the appearance of factorials and binomial coefficients is not that surprising. The inverse obtained from Lagrangian inversion, however, does not seem to show any known pattern in the coefficients.
– J. M. is not a mathematician
Jul 31 '17 at 6:10
add a comment |
I know this isn't a closed form, but I was interested in this question and have found one can relate the two functions together as series representations. I have written a short article here but this is the crux of it
Write begin{equation}
E(phi,k) = sum_{i=0}^infty frac{Q_i(k)}{(2i-1)!}phi^{2i-1}
end{equation}
where $Q_i(k)$ are polynomials in $k$, from the series expansion here we can get a finite form for these polynomials as begin{equation}
Q_n(m) = 2(-4)^nsum_{k=1}^n frac{(2k-3)!!}{k!}left(frac{-m}{8}right)^k sum_{j=0}^{k-1} binom{2k}{j}(-1)^{1-j}(j-k)^{2n}, ;;; n>0
end{equation}
with $Q_0(k)$ defined as $1$. Then you can write the inverse series using series reversion in a very similar manner to $E(phi,k)$ begin{equation}
phi(E,k) = sum_{i=0}^infty frac{R_i(k)}{(2i-1)!}E^{2i-1}
end{equation}
where the relation between the new polynomials $R_i(k)$ is given by the explicit reversion formula found at the bottom of this link, giving begin{equation}
R_n(k) = (2n)! sum_{tau_n=n}(-1)^{sigma_n} frac{prod_{i=1}^{sigma_n}2n+i}{prod_{j=1}^n k_j!}prod_{l=1}^n left(frac{Q_l(k)}{(2l+1)!}right)^{k_l}
end{equation}
where $sigma_n=k_1+k_2+k_3+cdots+k_n$, and $tau_n=k_1+2k_2+3k+3+cdots + nk_n$ and the sum is over all sets of indices $k_i$ that meet the requirement $tau_n=n$. I don't know if any nice simplifications or tricks can be made to reduce this to a functional form. The only numerical element here is converging the series to the desired accuracy.
answered May 25 '17 at 14:32
Benedict W. J. Irwin
1,463529
The incomplete elliptic integral of the second kind is known to be expressible in terms of a (bivariate) hypergeometric function, so the appearance of factorials and binomial coefficients is not that surprising. The inverse obtained from Lagrangian inversion, however, does not seem to show any known pattern in the coefficients.
– J. M. is not a mathematician
Jul 31 '17 at 6:10
add a comment |
I know this isn't a closed form, but I was interested in this question and have found one can relate the two functions together as series representations. I have written a short article here but this is the crux of it
Write begin{equation}
E(phi,k) = sum_{i=0}^infty frac{Q_i(k)}{(2i-1)!}phi^{2i-1}
end{equation}
where $Q_i(k)$ are polynomials in $k$, from the series expansion here we can get a finite form for these polynomials as begin{equation}
Q_n(m) = 2(-4)^nsum_{k=1}^n frac{(2k-3)!!}{k!}left(frac{-m}{8}right)^k sum_{j=0}^{k-1} binom{2k}{j}(-1)^{1-j}(j-k)^{2n}, ;;; n>0
end{equation}
with $Q_0(k)$ defined as $1$. Then you can write the inverse series using series reversion in a very similar manner to $E(phi,k)$ begin{equation}
phi(E,k) = sum_{i=0}^infty frac{R_i(k)}{(2i-1)!}E^{2i-1}
end{equation}
where the relation between the new polynomials $R_i(k)$ is given by the explicit reversion formula found at the bottom of this link, giving begin{equation}
R_n(k) = (2n)! sum_{tau_n=n}(-1)^{sigma_n} frac{prod_{i=1}^{sigma_n}2n+i}{prod_{j=1}^n k_j!}prod_{l=1}^n left(frac{Q_l(k)}{(2l+1)!}right)^{k_l}
end{equation}
where $sigma_n=k_1+k_2+k_3+cdots+k_n$, and $tau_n=k_1+2k_2+3k+3+cdots + nk_n$ and the sum is over all sets of indices $k_i$ that meet the requirement $tau_n=n$. I don't know if any nice simplifications or tricks can be made to reduce this to a functional form. The only numerical element here is converging the series to the desired accuracy.
answered May 25 '17 at 14:32
Benedict W. J. Irwin
1,463529
I know this isn't a closed form, but I was interested in this question and have found one can relate the two functions together as series representations. I have written a short article here but this is the crux of it
Write begin{equation}
E(phi,k) = sum_{i=0}^infty frac{Q_i(k)}{(2i-1)!}phi^{2i-1}
end{equation}
where $Q_i(k)$ are polynomials in $k$, from the series expansion here we can get a finite form for these polynomials as begin{equation}
Q_n(m) = 2(-4)^nsum_{k=1}^n frac{(2k-3)!!}{k!}left(frac{-m}{8}right)^k sum_{j=0}^{k-1} binom{2k}{j}(-1)^{1-j}(j-k)^{2n}, ;;; n>0
end{equation}
with $Q_0(k)$ defined as $1$. Then you can write the inverse series using series reversion in a very similar manner to $E(phi,k)$ begin{equation}
phi(E,k) = sum_{i=0}^infty frac{R_i(k)}{(2i-1)!}E^{2i-1}
end{equation}
where the relation between the new polynomials $R_i(k)$ is given by the explicit reversion formula found at the bottom of this link, giving begin{equation}
R_n(k) = (2n)! sum_{tau_n=n}(-1)^{sigma_n} frac{prod_{i=1}^{sigma_n}2n+i}{prod_{j=1}^n k_j!}prod_{l=1}^n left(frac{Q_l(k)}{(2l+1)!}right)^{k_l}
end{equation}
where $sigma_n=k_1+k_2+k_3+cdots+k_n$, and $tau_n=k_1+2k_2+3k+3+cdots + nk_n$ and the sum is over all sets of indices $k_i$ that meet the requirement $tau_n=n$. I don't know if any nice simplifications or tricks can be made to reduce this to a functional form. The only numerical element here is converging the series to the desired accuracy.
answered May 25 '17 at 14:32
Benedict W. J. Irwin
1,463529
answered May 25 '17 at 14:32
Benedict W. J. Irwin
1,463529
answered May 25 '17 at 14:32
Benedict W. J. Irwin
1,463529
answered May 25 '17 at 14:32
Benedict W. J. Irwin
1,463529
1,463529
The incomplete elliptic integral of the second kind is known to be expressible in terms of a (bivariate) hypergeometric function, so the appearance of factorials and binomial coefficients is not that surprising. The inverse obtained from Lagrangian inversion, however, does not seem to show any known pattern in the coefficients.
– J. M. is not a mathematician
Jul 31 '17 at 6:10
add a comment |
The incomplete elliptic integral of the second kind is known to be expressible in terms of a (bivariate) hypergeometric function, so the appearance of factorials and binomial coefficients is not that surprising. The inverse obtained from Lagrangian inversion, however, does not seem to show any known pattern in the coefficients.
– J. M. is not a mathematician
Jul 31 '17 at 6:10
The incomplete elliptic integral of the second kind is known to be expressible in terms of a (bivariate) hypergeometric function, so the appearance of factorials and binomial coefficients is not that surprising. The inverse obtained from Lagrangian inversion, however, does not seem to show any known pattern in the coefficients.
– J. M. is not a mathematician
Jul 31 '17 at 6:10
The incomplete elliptic integral of the second kind is known to be expressible in terms of a (bivariate) hypergeometric function, so the appearance of factorials and binomial coefficients is not that surprising. The inverse obtained from Lagrangian inversion, however, does not seem to show any known pattern in the coefficients.
– J. M. is not a mathematician
Jul 31 '17 at 6:10
add a comment |
Thanks for contributing an answer to Mathematics Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
Use MathJax to format equations. MathJax reference.
To learn more, see our tips on writing great answers.
Some of your past answers have not been well-received, and you're in danger of being blocked from answering.
Please pay close attention to the following guidance:
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f1123360%2finverse-of-elliptic-integral-of-second-kind%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
This paper might be interesting for you: Numerical computation of incomplete elliptic integrals of a general form by T. Fukushima and H. Ishizaki. Also the paper Numerical computation of inverse complete elliptic integrals of first and second kinds by T Fukushima.
– Loreno Heer
Jan 28 '15 at 16:07
yes inverse calculations are in the second paper
– Loreno Heer
Jan 28 '15 at 16:08
@sanjab: Thanks for these pointers. But as far as I can see, the first paper only does forward computation, and the second does inverse (using notation like $m_E$), but of the complete integral and solving for the parameter $m$ not the amplitude $varphi$. So at least as far as I can see, this is related but won't immediately help answering my question.
– MvG
Jan 28 '15 at 16:29
1
FWIW, Wolfram is looking for something nice, too.
– J. M. is not a mathematician
May 17 '16 at 3:40
N.B. the inverse of the Weierstrass elliptic function can of course be expressed in terms of $F$, or Carlson's $R_F$ if you prefer the symmetric integrals. Similarly, the Weierstrass and Jacobi functions are expressible in terms of each other (corresponding to the algebro-geometric notion of a quartic being equivalent to a cubic under suitable rational transformations.)
– J. M. is not a mathematician
May 17 '16 at 3:43