7
$\begingroup$

I am trying to manipulate a conformal map from the half-plane to a square $z \rightarrow w(z)$ defined by:

$$ w(z) = \int \limits^{z} dx \frac{1}{\sqrt{(1-x^2)x}} = \sqrt{2} \; F\left(\sqrt{z+1},\frac{\sqrt{2}}{2}\right)$$

where $F$ is the incomplete elliptic integral of the first kind according to this.

My problem is that the incomplete elliptic integral built into Mathematica as EllipticF does not seem to behave like this function, as its derivative is not singular at $z=1$ as the formula is telling us.

Proof that the Mathematica derivative is finite:

D[EllipticF[Sqrt[x + 1], Sqrt[2]/2], x] /. x -> 1 // N

gives me

0.63491

I am very unfamiliar with these functions; is it a convention difference between the two?

$\endgroup$
11
$\begingroup$

(I just knew someone would ask this someday...)

I had talked (ranted?) about this issue at some length here, so I'd like you to read that first. I'll just give an executive summary here: Mathematica uses the parameter convention, while the formula you found on Wikipedia is based on the modulus convention (quickly betrayed by the comma separating the two arguments).

In fact, even the convention in the first argument is different, so EllipticF[] is not exactly the Mathematica function you want. What you need here, in Mathematica's convention, is denoted as InverseJacobiSN[]:

D[InverseJacobiSN[Sqrt[x + 1], 1/2], x] /. x -> 1 // FullSimplify
   ComplexInfinity

After looking through my Schwarz-Christoffel notes, I'm not sure where the Wikipedia formula was obtained. In any case, allow me to present a short demonstration of the conventional conformal mapping between the upper half plane and a rectangle (though I cheat here and use the inverse map, JacobiSN[], in the following code):

a = 1/2; b = 1; (* 2 a × b rectangle *)
m = ModularLambda[I b/a];

{ParametricPlot[With[{w = JacobiSN[EllipticK[m] (u + I v)/a, m]},
                     {Re[w], Im[w]}], {u, -a, a}, {v, 0, 99 b/100},
                BoundaryStyle -> AbsoluteThickness[2], Mesh -> 15, 
                MeshShading -> {{Transparent}}, 
                MeshStyle -> {Directive[AbsoluteThickness[2], ColorData[97, 1]], 
                              Directive[AbsoluteThickness[2], ColorData[97, 2]]}, 
                PlotRange -> {{-5, 5}, {0, 10}}], 
 ParametricPlot[{u, v}, {u, -a, a}, {v, 0, 99 b/100}, 
                BoundaryStyle -> AbsoluteThickness[2], Mesh -> 15, 
                MeshShading -> {{Transparent}}, 
                MeshStyle -> {Directive[AbsoluteThickness[2], ColorData[97, 1]], 
                              Directive[AbsoluteThickness[2], ColorData[97, 2]]}]}
// GraphicsRow

upper half-plane to rectangle map

As a check of correctness, try restricting the range of v to {0, b/2}. The mapping on the left should yield a semicircular region with radius m^(-1/4).

$\endgroup$
  • 1
    $\begingroup$ P.S. You'll want to use ParametricPlot[] to visualize conformal mappings. $\endgroup$ – J. M. is away Nov 1 '15 at 20:25
  • $\begingroup$ @ J.M. how do you relate the integral w to the function InverseJacobiSN? Plotting them shows they are rather different. $\endgroup$ – Dr. Wolfgang Hintze Nov 2 '15 at 7:26
  • $\begingroup$ @Dr. Hintze, due to branch cuts, things are a little murkier. In fact, the actual mapping function is something like Sqrt[2] (EllipticK[1/2] - InverseJacobiSN[Sqrt[1 - x], 1/2]). $\endgroup$ – J. M. is away Nov 2 '15 at 7:37
  • $\begingroup$ @ J.M. In fact we have -Im[Sqrt[2] (EllipticK[1/2]-InverseJacobiSN[Sqrt[1-z],1/2])] = f for z>=1. Hence my function f gives the complete solution for the intergal w in terms of incomplete (F) and complete (K) elliptic integrals - as I said in my answer. $\endgroup$ – Dr. Wolfgang Hintze Nov 2 '15 at 7:50
  • $\begingroup$ Well, the OP uses $1-x^2$, while you use $x^2-1$; maybe that explains the discrepancy, @Dr. Hintze? $\endgroup$ – J. M. is away Nov 2 '15 at 7:52
3
$\begingroup$

I would prefer to have consistency within the framework of Mathematica and hence start off with the original definition of your function as an integral taken here between 1 and z:

f = 
 Integrate[1/Sqrt[x (x^2 - 1)], {x, 1, z}, Assumptions -> z > 1]

(*
Out[35]= -2 EllipticF[ArcSin[1/Sqrt[z]], -1] + 2 EllipticK[-1]
*)

This is obviously different from the result you quote from Wikipedia. Hence the validity of the latter should be clarified first.

$\endgroup$
  • $\begingroup$ It seems that there's a confusion in the definitions, as J.M. has just pointed out. $\endgroup$ – Dr. Wolfgang Hintze Nov 1 '15 at 20:17
  • $\begingroup$ This might be of interest. $\endgroup$ – J. M. is away Nov 2 '15 at 8:05

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.