I am currently working to try and solve the integral:
∫π20√1−a2cos(x)2=√1−a2E(1+1−1+a2)
where E(m) is the complete elliptic integral of second kind with the parameter m=k2 and $0.
This is the soultion I found using WolframAlpha
In this case m is always smaller then 0, yet m is supposed to be the square of some k.
The implementation I found in NumericalRecipies seems to only accept k as argument, evident through a test, plotting values of the positive argument-branch and comparing them with WolframAlphas elliptic function.
However, the NR implementation doesn't use any sort of complex variables (therefor, won't accept them as input).
Do I need another implementation, am I completely off with this.
I am very confused and uncertain about this.
Help would be greatly appreciated, Thanks.
Answer
I found the solution to my problem here:
http://analyticphysics.com/Mathematical%20Methods/A%20Miscellany%20of%20Elliptic%20Integrals.htm
To quote:
E(−m)=√m+1E(m+1m)
Simply building an if clause into my driver to choose −k whenever k<0 and then use the formula indicated above gave me the result WolfrAmalpha shows in wolframalpha.com/input/?i=plot+EllipticE%5Bx%5D.
The discussion here was what got me to ask the right question. Thanks Jack D'Aurizio.
No comments:
Post a Comment