Area enclosed by an equipotential curve for an electric dipole on the plane I am currently teaching Physics in an Italian junior high school. Today, while talking about the electric dipole generated by two equal charges in the plane, I was wondering about the following problem:

Assume that two equal charges are placed in $(-1,0)$ and $(1,0)$.
There is an equipotential curve through the origin, whose equation is
given by:
$$\frac{1}{\sqrt{(x-1)^2+y^2}}+\frac{1}{\sqrt{(x+1)^2+y^2}}=2 $$ and
whose shape is very lemniscate-like:



Is there a fast&tricky way to compute the area enclosed by such a curve?

Numerically, it is $\approx 3.09404630427286$.

HINT:
Almost sure area can be evaluated between limits $(u_1,u_2),(v_1,v_2)$ for holomorphic functions of complex variables in conformal maps.
The equipotentials and lines of force form an orthogonal net  of curvilinear rectangle curves (resemble Cassinian Ovals, but different.. one out of which is shown above), and force lines which are all concurrent hyperbolas passing through the pole of same polarity / charge and asymptotic to y-axis. For these $(u,v)$ parameter plots that satisfy Laplace equation all we need before evaluating integral of  area is to identify limits of u and v for any oval chosen. (shall upload equation , if it is found in sites e.g. 2d curves.com)
EDIT1:
The closest that I can get to the function is sketched by Mark McClure in another context for  $ f(z)= z^2 (z-1)^3 $:
Orth Traj to Hyperbolae
I think we should choose another complex variable part so that our hyperbolae pass through $ (\pm 1,0)$ instead of $ ( 0,0),(1,0)$ as given above. We are interested in a particular orthogonal trajectory doubly intersecting at origin. Calculation of area may be easier.
