Quadrature for quadrics

We make a systematic investigation of quadrature properties for quadrics, namely integration of holomorphic functions over planar domains bounded by second degree curves. A full understanding requires extending traditional settings by allowing domains which are multi-sheeted, in other words domains which really are branched covering surfaces of the Riemann sphere, and in addition usage of the spherical area measure instead of the Euclidean. The first part of the paper discusses two different points of view of real algebraic curves: traditionally they live in the real projective plane, which is non-orientable, but for their role for quadrature they have to be pushed to the Riemann sphere. The main results include clarifying a previous theorem (joint work with V. Tkachev), which says that a branched covering map produces a domain with the required quadrature properties if and only it extends to be meromorphic on the double of the parametrizing Riemann surface. In the second half of the paper domains bounded by ellipses, hyperbolas, parabolas and their inverses are studied in detail, with emphasis on the hyperbola case, for which some of the results appear to be new.


Introduction
This note is inspired by a question by Henrik Shahgholian on why exterior of ellipses and parabolas in the plane are null-quadrature domains, but not the exterior of hyperbolas. Partial answers, from the point of view of potential theory and PDE (partial differential equations), can be read off from papers of Shahgholian himself and other mathematicians, but here we shall try to give a more philosophical answer in terms of what we call multi-sheeted algebraic domains, or quadrature Riemann surfaces.
The paper is more oriented towards real algebraic geometry than to PDE, but the two research areas are linked via potential theory, specifically quadrature domains. The title of the paper is a little of a play with words: quadrature refers to numerical integration via the latin noun 'quadrature', making square-shaped. After having made a region square-shaped, or dived it into small squares, one can read off the area of any arising square by taking the power two of (squaring) its side length. For related reasons, 'quadrics' refers to curves described by equations of degree two. The word 'conic' is more or less synonymous, but has a more geometric slant as the intersection between a cone and a plane. The actual subject of this paper eventually becomes integration of holomorphic functions over domains bounded by algebraic curves of degree two, essentially ellipses, hyperbolas and parabolas.
While the subject of conic sections is very old, potential theory dates back "only" to Newton's theory of gravitation (1687). But already Newton himself realized that conic sections have a special role in gravitational theory. Besides the role of these curves in planetary motions there is the theorem that ellipsoidal homoeoids (elliptic homothetic shells) give zero gravitational attraction inside the hole. This was proved in three dimensions by Newton himself (at least to some part) and later extended by P.S Laplace and J. Ivory, see [30,39], and in particular the recent paper [27] which contains a wealth of geometric aspects. For the life of James Ivory, see [8]. Several papers produced by the V. Arnold school, for example [2,3,47,4] contain further results for hyperboloids, and [37,11,14,44,12,42,28,29,6,32] give links between null quadrature domains and regularity theory for free boundaries.
The property of ellipsoidal shells creating no interior gravity can be understood in terms of a Laplacian growth process (or Hele-Shaw moving boundary problem) which preserves harmonic moments and which is related to quadrature domains. A few general sources are [46,36,21]. See [20] for a an overview of the theory of quadrature domains. Topology of quadrature domains, partly in connection with Hilbert's sixteenth problem, is studied in [34].
In the present paper we make a systematic investigation of quadrature properties for planar domains bounded by second degree curves, including to some extent inversions of such domains. We get rid of difficulties caused by exterior domains having infinite area by working with spherical area measure instead of Euclidean area measure. This does not change the class of quadrature domains as a whole, in that sense the change is inessential, however some computations become a more involved.
Another difference compared to traditional treatments is that we allow the domains to be "multi-sheeted", i.e. to be covered several times by a branched covering map from some uniformizing surface. The latter is to be one half of a symmetric, or "real", Riemann surface of "dividing type". This simply means the Riemann surface (assumed compact) has an anti-conformal involution and that it becomes disconnected after removal of the fixed points of the involution.
The boundaries of so arising multi-sheeted quadrature domains, or quadrature Riemann surfaces in the terminology of Sakai [40], are algebraic curves, and essentially one half of all algebraic curves arise in this way. Therefore much of our treatment falls into the subject of real algebraic geometry. But here arises a kind a dichotomy: real algebraic curves have their natural loci in the real projective plane RP 2 while the potential theoretic considerations rather take place in the Riemann sphere CP 1 ∼ = C ∪ {∞}.
Looking with positive eyes this actually makes the theory more interesting and rich, and there are even some surprises to come up. While the theories for the ellipse and parabola are rather straightforward, the hyperbola cannot be discussed without running into multi-sheeted domains and dichotomies between RP 2 and CP 1 : in RP 2 the hyperbola is an everywhere smooth curve (even at points of infinity), it is an "oval" in the terminology of real algebraic geometry. The interior of the oval are the two pieces containing the focal points, while the exterior is the intermediate region, which topologically is a Möbius strip (recall that RP 2 is non-orientable). See Figure 3. The corresponding multi-sheeted region in CP 1 which has potential theoretic significance consists of one of the pieces containing the focal point, having multiplicity two and with the focal points as branch point, continued with multiplicity one across the curve, up to the other component of the curve. See Figure 6. And in CP 1 the hyperbola has a singularity at infinity, in fact it becomes a lemniscate upon inversion.
After an initial discussion of algebraic curves in projective spaces, partly in terms of corresponding projection maps (Sections 2 and 3), we turn in Section 4 to basic definitions and general properties of quadrature domains in our setting. This part extends, and partly clarifies, a previous treatment [23]. The single main result of the paper is Theorem 4.1. Some computational details for its proof are deferred to an appendix Section 9.
After these general matters in Sections 2-4 we discuss systematically in the remaining Sections 5-8 the different special cases arising for second degree curves: the ellipse, the hyperbola, the parabola, and inversions of these curves. These sections are to a large extent computational in nature, but the computations are not always easy. One of the main outcomes is that hyperbolas are indeed multi-sheeted quadrature domains. If one works with Euclidean area measure instead of the spherical, then they are multi-sheeted null quadrature domains The text is accompanied with hopefully helpful pictures, constructed in most cases by "TikZ". These are however not computationally exact, they are mainly for illustrational purpose.

The genus and Riemann-Hurwitz formulas
For a plane algebraic curve P (z, w) = 0 of degree d and genus g the genus formula says that g + number of singular points = (d − 1)(d − 2) 2 . (1) It is assumed in this formula that the curve is irreducible, and for our considerations it will in addition be real in the sense that P (z, w) is real-valued when w =z. The coordinate functions z and w can be viewed as meromorphic functions on that compact Riemann surface M which uniformizes the curve, and thereby these coordinate functions can also be viewed as branched covering maps of M onto the Riemann sphere. For such a map, let m denote the number of sheets and b the number of branch points (in general these will be different for z and w). Then the Riemann-Hurwitz formula says that We refer to [31,15,7,13] for details. One fundamental difference between the genus formula and the Riemann-Hurwitz formula is that the genus formula strictly speaking refers to the completion of the curve in two dimensional complex projective space CP 2 (for example, one has to count the singular points within that space), and from that point of view the corresponding real algebraic curve P (z,z) = 0 becomes a subset of the Figure 1: The conic sections. Picture copied from mathhints.com. real projective space RP 2 . On the other hand, (2) refers to a meromorphic function z : M → CP 1 , and this makes essentially the same curve become a subset of the Riemann sphere CP 1 . The finite parts of RP 2 and CP 1 can both be naturally be identified with R 2 ∼ = C, but when it comes to points of infinity they differ considerably: CP 1 has only one point of infinity whereas RP 2 has a whole projective line of infinities, representing asymptotic directions. In addition, RP 2 is non-orientable. In fact, denoting by S 2 the unit sphere in R 3 , there are homoemorphisms CP 1 ∼ = S 2 and RP 2 ∼ = S 2 /(antipodes identified).
In the case of quadrics (conics) d = 2, hence the genus formula says that g = 0 and that there are no singular points. Thus b = 2(m − 1), which gives b = 2 in most cases to be discussed below. One exception is the circle, for which m = 1 and so b = 0. We shall also discuss curves obtained by making inversions (such as the antipodal map z → −1/z) of quadrics. This will lead to curves like hippopedes (special cases of hypocycloids), lemniscates and cardioids. Such inversions will not change the genus, since the curve is still uniformized by the same Riemann surface, but the degree of the curve may change. Typically it will change from d = 2 to d = 4, and then the genus formula tells that there will be three singular points. These will turn out to have roles of being quadrature nodes, "special points", intersections between smooth branches, and cusps, depending on the curve. In general, some of the singular points may be invisible in the real.

General
We shall be concerned with real algebraic curves, the meaning of "real" being that the polynomial Q(x, y) = P (x + iy, x − iy) has real coefficients. Then the coefficients of P (z, w) has a hermitean symmetry, and the complex curve {(z, w) ∈ C 2 : P (z, w) = 0} admits the involution (z, w) → (w,z). In general we will denote all anticonformal involutions by the same letter J, the exact meaning then to be understood from the context. The real locus of P , namely is the fixed point set of the mentioned involution. The Riemann surface M associated to the curve P (z, w) = 0 has a corresponding anti-conformal involution J : M → M . The locus (3) can be obtained from the fixed point set Γ M of J, which can also be viewed as the symmetry line of M , essentially as the image z(Γ M ) under the coordinate function z, the latter considered as a meromorphic function on M . It can alternatively be obtained using w, as The above picture is actually somewhat simplified, for example one has to carefully compactify the above algebraic curves in projective spaces, and it also turns out that V may contain accidental points, not coming from Γ M . We elaborate these details below.
The two dimensional real projective space is the set of equivalence classes (or ratios) representing lines through the origin in R 3 . As indicated, only the ratio counts, i.e. [ t : x : y ] = [ λt : λx : λy ] for any λ = 0, λ ∈ R. The finite plane R 2 can be naturally embedded in RP 2 via (x, y) → [ 1 : x : y ]. The complex projective plane CP 2 is defined similarly, just that t, x, y, λ are allowed to be complex numbers. However, it will be more suitable in our context to name the coordinates t, z, w and identify the real subspace via z = x + iy, w = x − iy, where x, y ∈ R. We shall thus work with The set of fixed points for the involution is the "real diagonal" in CP 2 , for which t ∈ R and w =z. Thus it consists of the points In this version of the real projective space the third component, z = x − iy, is redundant and can be discarded. Therefore the real projective space can effectively be identified with the set of equivalence classes where two pairs represent the same point if they are related by a nonzero real factor. The subscript red stands for "reduced", meaning that an inessential component has been discarded. The definition of set RP 2 red reminds of that of the one dimensional complex projective space a model for the Riemann sphere. Here it is always possible to choose t to be real, in fact the two values t = 0 and t = 1 cover all possible equivalence classes. It follows that there is a natural and surjective map which takes a point [ t : z ] ∈ RP 2 red to the point with the same name in CP 1 . But in CP 1 space it may belong to a bigger equivalence class. This is more precisely the case when t = 0: in CP 1 all points with t = 0 are equivalent to [ 0 : 1 ], i.e. there is only one point of infinity, while for RP 2 red , two pairs [ 0 : z 1 ] and [ 0 : z 2 ] represent the same point of infinity if and only if z 2 = λz 1 for some λ ∈ R \ {0}. Thus one can take representatives of the form [ 0 : e iϕ ], observing then that two ϕ represent the same point if and only if they differ by an integer multiple of π (not 2π). The latter remark is related to the fact that RP 2 red is a non-orientable surface, topologically a "cross-cap". The Riemann sphere, CP 1 ∼ = S 2 is, on the other hand, orientable.
In the sequel we shall also semi-complexify the real spaces by identifying R 3 with R ⊕ C with coordinates (t, z) = (t, x + iy), and also identify RP 2 with RP 2 red as in (4).

Ray representations
In general, projective spaces like RP 2 , CP 2 , RP 1 and CP 1 can be defined as the set of lines through a fixed point in a vector space of one dimension higher than the index given. For example, or, choosing another point, the north pole N = (1, 0, 0), Each of the above two representations (6) and (7) has its own advantages. In the first case, each line K intersects the unit sphere S 2 at two antipodal points, p and −p. Switching to RP 2 red this gives a bijective map RP 2 red → S 2 /(antipodes identified).
Here one hemisphere can be discarded, say the "southern" hemisphere, which gives the description RP 2 red ∼ = (northern hemisphere)∪(equator with antipodes identified), presenting RP 2 red indeed as a "cross-cap". As a further concretization of (6), each line K which is not horizontal intersects the plane {t = 1} at exactly one point, (1, z) ∈ R ⊕ C. This point corresponds to the projective coordinates for the line when written as [ 1 : z ], hence corresponds directly to a complex number z, which we denote z K = z K (p) ∈ C, with p as in (8). Obviously z K (−p) = z K (p). Any horizontal line represents a points of infinity for RP 2 red , and such a line has a projective representation as [ 0 : e iϕ ] with ϕ ∈ R taken modulo π. See Figure 2 for an illustration.
The second choice (7) is connected to the stereographic projection, which rather identifies S 2 with C ∪ {∞}. Each non-horizontal line L through the north pole N = (1, 0) ∈ R ⊕ C intersects S 2 in one more point, say p, and it also intersects the plane {t = 0} at one point, say (0, z), where then z = z L (p) ∈ C. This defines the stereographic projection which extends to all of S 2 by setting z L (N ) = ∞. Thus the horizontal lines L give via the stereographic projection only one point of infinity, despite there as as many horizontal lines of L as there are of K. On the computational side, the relation between the components of p ∈ S 2 \ {N } and z = z L (p) ∈ C in the stereographic projection is In the other direction we have if we identify p as the point p = (x 0 , The Riemannian metric on S 2 as inherited from the ambient Euclidean space is ds 2 = dx 2 0 + dx 2 1 + dx 2 2 , and it is straightforward to show (and of course well-known) that this becomes when expressed in terms of z = z L (p). The fact that this is of the form ds 2 = λ(z) 2 |dz| 2 for some scalar factor λ(z) > 0 means that the map p → z L (p) is conformal as a map between S 2 and the Riemann sphere Thinking next of the Euclidean coordinates (9) for p ∈ S 2 as projective coordinates they are, in this role, This means that the complex number 2z L /(|z L | 2 − 1) represents the line K, provided |z L | = 1. Accordingly we have, with z K = z K (p) as above, the latter showing that 1/z K is a harmonic function of 1/z L . Expressed in terms of p = (x 0 , x 1 + ix 2 ) we have In contrast to z L (p), the map p → z K (p) is not conformal. Indeed, the relationship (11) between z L and z K is not analytic (or anti-analytic), hence z L and z K cannot both be conformal. When |z L | = 1 the line K is horizontal, which means that any such z L represents a point of infinity in RP 2 , then with z L and −z L representing the same point. We finally note (or recall) that We summarize everything as follows (see again Figure 2). Figure 2: Illustration of projections p → z L (p) and p → z K (p) of sphere onto complex and real planes with their infinities (two-dimensional view).
Lemma 3.1. Points p = (x 0 , x 1 + ix 2 ) on the unit sphere S 2 ⊂ R ⊕ C have two natural representations as (extended) complex numbers, namely as z L and z K defined below. 1) Via the stereographic projection, which is a conformal map Alternatively, considering CP 1 to be the target for stereographic projection, In the converse direction p can be recovered from z L = z L (p) by 2) As a map onto the real projective plane, bijective after identification of antipodes: This is defined by and ϕ ∈ R/πZ is defined by The maps in the lemma are, in summary, For finite points the associations are given by

Classical quadrature domains
The traditional definition of a quadrature domain [20] is that it is a bounded domain (or open set) Ω ⊂ C for which there exist finitely many points z 1 , . . . , z m ∈ Ω, integers m j ≥ 1, and coefficients a kj ∈ C such that the identity holds for all integrable analytic functions h in Ω. The simplest example is the unit disk D, for which the mean-value property for analytic functions gives the quadrature identity Quadrature domains are known to be bounded by algebraic curves. In fact, if Ω satisfies (15) then ∂Ω is a real algebraic variety as in (3), possibly minus finitely many points, which then are singular point of the variety. Solving the algebraic equation P (z,z) = 0 in (3) for z gives the Schwarz function S(z) for ∂Ω. Then P (z, S(z)) = 0 holds identically, and S(z) =z on ∂Ω.
The latter equation is the defining property of a Schwarz function, generally required to be holomorphic in an at least a one-sided neighborhood of ∂Ω. The main statement [1] in this context is that Ω is a quadrature domain if and only if ∂Ω has a Schwarz function which extends to be meromorphic in all of Ω. A few sources in book form for quadrature domains are [9,38,45,46]. In the latter text quadrature domains as in (15) are named an algebraic domains. This is an appropriate terminology even though quadrature domains only represent a small subclass of all domains bounded by algebraic curves. In order to diminish the gap attempts were made in [23] to extend the class of quadrature domains/algebraic domains by first of all working with the spherical area measure on the Riemann sphere instead of Euclidean measure in the plane, and secondly by allowing multiply covered domains. The so arising class of multi-sheeted algebraic domains, or by another name, used in [40], quadrature Riemann surfaces, covers essentially one half of all domains bounded by algebraic curves. This class consists of those multiply covered domains which are images of one side of a symmetric (or "real") Riemann surface of "dividing type" under a function which is meromorphic on the entire surface. A symmetric Riemann surface is said to be of dividing type if the symmetry line disconnects the surface onto two halves.
One simple example of a domain bounded by an algebraic curve, which is not a quadrature domain in the sense of (15), is the interior of any non-circular ellipse, for example the domain defined by On the other hand, the exterior (with respect to the complex plane) of the ellipse is an unbounded domain admitting an identity (15), in fact it is a null quadrature domain: holds for all integrable holomorphic functions h in the domain. See [37] for a complete result in this direction. And in a general sense the ellipse can be rescued to belong to the class of multisheeted algebraic domains.
A domain which cannot be rescued by any means is The two Riemann surfaces associated to this domain, namely that defined by the algebraic equation for the boundary and the Schottky double of the domain, are not canonically isomorphic in the way they need in order to be a quadrature domain. As we shall try to argue below, almost everything goes wrong with the curve in (19). Generic algebraic curves are smooth, i.e. have no singular points. However algebraic curves bounding quadrature domains as in (15) usually have many singular points when completed in complex projective space. As follows from a classification carried out in [18] the right member in (15) is a manifestation of some of these singular points. The functional in right member accounts more precisely for n(n − 1) of the singular points, n = m k=1 m k being the order of the quadrature identity. The curve bounding the domain in (19) can be seen to have no singular points at all, even when viewed as an algebraic curve in CP 2 . This fact is in itself a sign that it cannot bound any quadrature domain.
Using the genus formula (1), the absence of singular points in addition shows that the mentioned algebraic curve has genus three, while the real locus, defining the boundary in (19), consists of only one closed curve. One single closed curve cannot separate a surface of genus three into two equal halves, hence the Riemann surface for (19) cannot be of "dividing type". Finally, the Schottky double of the domain in (19) obviously has genus zero, not three.
Somewhat similar results as those for ellipses hold for parabolas, but hyperbolas are different, and indeed quite interesting. We start by considering just one example, namely the open set defined by As a subset of the complex plane, or of the Riemann sphere, it has two components. As seen from the point of infinity, for example after an antipodal inversion, the two components look like the two holes in the figure "∞", i.e. the curve becomes a lemniscate. However, when viewed as a subset of RP 2 red the set (20) becomes connected via all asymptotic lines with directions y = cx, |c| < 1. The curve x 2 −y 2 = 1 is an "oval" in a terminology used in real algebraic geometry, see [41,16,24], and the domain defined by (20) is in that terminology the interior of it.
The exterior of the oval, i.e. the complement of (20) in RP 2 red , is topologically a Möbius strip. This is completely consistent with RP 2 red being topologically a "cross-cap", which is non-orientable and can Möbius strip "oval" "oval" interior of "oval" (topological disk) interior of "oval" (to be glued "upside-down" with the right hand part) be obtained by sewing a disk and a Möbius strip along their unique boundary curves. One can see the topology directly by providing the two branches of the hyperbola with positive directions downwards (decreasing y). This does not seem correct for a boundary of a domain, but indeed it is in this case, as we try to clarify in Figure 3. For example, the interior of the oval consists of two pieces in the picture, but these come with opposite orientation because of the identification of diametrically opposite points infinitely far away.
The hyperbola x 2 − y 2 = 1 is also the conic obtained by intersecting of the cone with the plane t = 1. Here the cone has two roles: besides being a choice of standard cone as in Figure 1, it gives the equation for our special hyperbola expressed in homogenous coordinates. Thinking thus of t, x, y as homogenous coordinates and changing names to is to satisfy Figure 4: Unusual picture of hyperbola: the boundary curves, as in (21), of two disks become a hyperbola when antipodal points on the sphere are identified. The sphere is in this role a double cover of the projective plane.
This gives the two circles representing the intersection of S 2 with two parallel planes. See Figure 4. The interpretation of the above is simply that we have identified, via (21), the tracks ±p ∈ S 2 of the hyperbola x 2 − y 2 = 1 under the identification z K (p) = x + iy in Lemma 3.1. Note that the two circles become only one circle in S 2 /J and that the intermediate part,

Quadrature domains with respect to spherical measure
Leaving now the strict realm of (15) and aiming at being able to treat domains with infinite area without loss of test functions, we shall replace ordinary area measure with spherical area measure. We identify this with the two-form here scaled so that the entire space C ∪ {∞}, representing the sphere S 2 , has area one.
If h(z) is analytic in a bounded domain Ω ⊂ C and if we set S(z) = z on ∂Ω, as in (17), then partial integration gives From this identity it follows that Ω is a quadrature domain for the spherical measure if and only if the differential S(z)dz 1+zS(z) has a meromorphic continuation to all of Ω. This holds if and only if S(z) is itself meromorphic in Ω, hence Ω is a spherical quadrature domain if and only if it is a Euclidean quadrature domain. The spherical quadrature identity then becomes where the right-hand side is a functional, as acting on h, on the same form as the right member in (15). A different usage of partial integration gives .
Here one sees that if S(z) has no pole at z = 0, then the two terms in the final member both contribute with residues at z = 0, but that these contributions cancel each other. Thus one can, if S(z) is regular at z = 0, write the final quadrature formula as In order to discuss multiply covered domains the coordinate z has to be substituted by a function z = f (ζ), where ζ is a coordinate on a uniformizing Riemann surface M , assumed to be "real" and of "dividing type". These assumptions on M mean more precisely that M shall be compact and be provided with an anti-conformal involution J : M → M such that, denoting by Γ M the set of fixed points of J, M \ Γ M has exactly two components, to be denoted M + and M − . It is assumed that M itself is connected. One of the components, say M + , is selected to be a model for the multi-sheeted domain Ω.
With M as above, any non-constant meromorphic function f on M gives rise to a spherical quadrature domain. Such a function may be viewed as a holomorphic and possibly branched covering map M → CP 1 , and the resulting quadrature domain Ω is to be identified with f (M + ) with appropriate multiplicities. The multivalued analytic functions on Ω allowed as test functions in quadrature identities will be those which become single-valued when lifted to M + . Thus the test functions will effectively be integrable analytic functions on M + itself.
We shall sometimes refer to M as the 'parameter' plane (or ζplane) and the target space CP 1 for f as the 'physical' plane' (or zplane). In terms of the above notations we formalize the concept of a multi-sheeted algebraic domain exactly as in [23]: holds for all functions h holomorphic and integrable with respect to the indicated spherical measure in M + . Here α kj ∈ C, ζ k ∈ M + are fixed quadrature data, and the derivatives are to be taken with respect to suitable local coordinates around the quadrature nodes. One can also express the quadrature property (25) without specifying the coefficients α kj by saying that the left member shall vanish for all holomorphic functions h which vanish at the points ζ k to order at least m k (1 ≤ k ≤ m).
The "easy" implication in the above statement, saying that f being meromorphic on all of M implies a formula of type (25), is of course mainly a matter of computation (partial integration combined with a residue calculus). However, it actually is slightly more delicate than one may at first expect. We shall therefore rework this computation, which was not fully carried out in [23].

Residue considerations
A first remark is that one may avoid certain troubles by rotating the Riemann sphere. Möbius transformations which preserve the spherical metric are of the form These orientation preserving rigid transformations can be combined with complex conjugations so that inversions like w = 1/z and antipodal map w = −1/z are included. Post-composing meromorphic maps like f above with such Möbius transformation allows us to rotate the Riemann sphere into suitable positions.
In particular one can arrange that f has no poles on the symmetry line Γ M . Indeed, if f has such a pole one may consider instead any function without such poles. Then and one may work with g in place of f . The above also means that the spherical area form itself is not sensitive for poles, or any other specific values, of f . Since the spherical area form, being a 2-form, is automatically closed, it is locally exact. For example, away from poles of f we can write Similarly, away from zeros of f , However the spherical area form cannot be globally exact because Here m is the degree of f (the number of zeros or poles), and we have used that the unit sphere has total area one. Starting from the right member of (29) and computing the exterior derivative there, noting that 1/(1 + ff ) is always a smooth function and taking, in remaining factors, distributional contributions carefully into account and using that df ∧ df = 0, gives Thus, rearranging terms, In a similar manner, In the above expressions the final terms consist of pure distributional contributions (Dirac currents): The right member here is the 2-form current corresponding to the divisor of f , ω k denoting the zeros and ξ j the poles of f , both repeated according to multiplicities. Localized at a zero of f , for example, the identity expresses that 1/πz is a fundamental solution for the Cauchy-Riemann operator: Written in an invariant way the latter relation becomes indicating also that Dirac measures shall be considered as 2-form currents. Compare Lemma 1 in [22]. Inserting (33) into (31) and (32), taking into account that the factor in front of d df /f kills in each case half of the Dirac currents, gives the following formulas. Proposition 4.1. Let f be a meromorphic function on M of degree m, let ω k be its zeros and ξ j its poles, both repeated according to multiplicities. Then, in the sense of currents, Remark 4.1. The presence of the last terms in (34) and (35) was slightly slipped over in [23]. One the other hand, those same terms but with the opposite sign, usually appear also in the first terms and therefore eventually cancel. See Theorem 4.1 below.  (27), the η j are the poles of the function g in (27), and in terms of g we have

Remark 4.3.
Instead of treating poles within the framework of distributions or currents one may cut small holes around singularities and then get additional boundary integrals, like in the derivation of (24).
Integrating over M + in Proposition 4.1, using Stokes' theorem together withf = f * on ∂M + , gives the following versions of the quadrature identity (25).
It is assumed, in accordance with Remark 4.2, that the sphere has been rotated so that zeros and poles of f on ∂M + are avoided. The first sums in the right members above are taken over all singular points in M + of the differentials indicated. The second sums (with explicit zeros and poles) are generically contained, but with opposite signs, also in the first sums and therefore usually cancel with terms from these.
The last statement can be made precise under the assumption that f * does not vanish at the poles of f and vice versa. In that case it is enough to sum over those points ζ ∈ M + for which 1 + f (ζ)f * (ζ) = 0, and the complete identities then become Proof. The formulas (36), (37) follow immediately on using Stokes' theorem together with Proposition 4.1. Note that the factor h can be multiplied to (35) and then goes inside the differential in the first term of the right member. For convenience we repeat the computation for (36) in one sequence of equalities: The use of Stokes' theorem can be justified by using Ahlfors-Bers mollifiers [5] in a way which was introduced in [38] and which is now a standard tool in the theory of quadrature domains. The proof of (38), (39) under the stated assumption, i.e. confirmation of the presence of cancellations, consists of a lengthy but straightforward local analysis of the differentials near an arbitrary point. In order not to interrupt the main lines of presentation we defer these details to an appendix, see Section 9.
Remark 4.4. One may scale the spherical area measure with a parameter ε > 0 to approach Euclidean area measure in the limit ε → 0.
Using that which is a generalization of (32), one obtains the quadrature identity Like in the statements of the theorem, the contributions from the poles ξ j cancel under natural assumptions. In the limit ε → 0 the identity takes the traditional form There may still be contributions from poles of f , but these come from solutions ζ ∈ M + of 1 + ε 2 f (ζ)f * (ζ) = 0 which have become poles of f in the limit ε → 0. On the other hand, poles of f make fewer functions h be allowed in the left member. For example, a pole of order one for f requires h to have a zero of order two at the same point (unless compensation comes from f * ), and therefore contributions from poles of f will usually not show up in the formula. This is what happens for null quadrature domains.
When n = 1 or n = 2, f is univalent on M + and the result can be written directly in the image domains Ω 1 = {Im z > 0} and Ω 2 = C \ [0, +∞) as follows: Here h(z) is considered as a function of z = f (ζ). In the second formula, Ω 2 can be replaced by C\I where I is an arbitrary closed subset of R, and it then actually follows from the first formula together with the corresponding formula for the reflected domain J(Ω 1 ). One may notice that ∂Ω 2 , and ∂(C \ I) in general, is not a full algebraic curve. As remarked in [17] (Lemma 1 and discussions after it) this is an exceptional situation showing up in the present case because f * = f , having as a consequence that the pair f and f * does not generate the field of meromorphic functions on M when n ≥ 2.
With Euclidean area measure, in place of spherical, both Ω 1 and Ω 2 become null quadrature domains.

The disk
The most basic quadrature domain, the unit disk D = {|z| < 1}, is bounded by the real algebraic curve Γ : x 2 + y 2 = 1.
As a complex algebraic curve in CP 2 it can be identified with the Riemann sphere. The Schottky double of D is also a model of the Riemann sphere, thought of as C ∪ {∞} with anti-conformal involution J(z) = 1/z, which has Γ as its fixed point set. The exterior of Γ is (as for the conformal structure) just another copy of D, and can be identified with the backside in the Schottky double. The fact that the two mentioned spheres are canonically identical, as Riemann surfaces with an anticonformal involution, is crucial for the property of D being a quadrature domain. This is a general statement applicable to all classical quadrature domains as in (15). The quadrature identities for the disk and its exterior D e are immediate: These identities can be identified with (41) after suitable rotations of the sphere, as in (26). The curve Γ can be naturally viewed also as a curve in the real projective space RP 2 red , which topologically can be identified with the sphere S 2 with antipodes identified. What is outside the curve then becomes topologically a Möbius strip, hence a non-orientable surface. This is the generic situation in real algebraic geometry: a simple closed algebraic curve, an "oval", which divides the space into a topological disk and a Möbius strip. Hilbert's sixteenth problem concerns the possible scenarios as for ovals inside each other in the case of real algebraic curves having several components, see [41,48,16,24].

The true ellipse
The standard ellipse in the plane is, with a > b > 0, Γ : x 2 a 2 + and it bounds the domain In terms of the complex coordinate the equation (43) becomes and solving forz gives a multivalued function with two branches as Here ±c are the foci, with The branch S − (z) corresponding to the minus sign (taking the square root to be positive for large positive values of z) equalsz on Γ and is single-valued in the region outside the ellipse, and also inside Γ up to the slit [−c, c]. Hence this branch is the actual Schwarz function of the curve Γ. The other branch, S + (z), is also single-valued in C \ [−c, c], and it analytically continues through the slit, to join S − (z) on the other side, compare Figure 5. However it does not agree with z on Γ. As z → ∞ both branches tend to infinity linearly.
It follows that S(z) is an algebraic function which becomes singlevalued on a two-sheeted Riemann surface M over C ∪ {∞} with branch points over ±c, and with each branch having a simple pole over z = ∞. Moreover, it uniformizes the algebraic equation P (z, w) = 0, where in the sense that P (z, S(z)) = 0 holds identically. The Riemann surface M is compact, has genus zero and has no singular points. Therefore M can be directly identified with the locus P (z, w) = 0, provided one compactifies the latter. This is to be done in two dimensional projective space, i.e. one introduces the homogenization P (t, z, w) = t 2 P (z/t, w/t) of P and defines loc P = {[ t : z : w ] ∈ CP 2 : P (t, z, w) = 0}.
Having genus zero, M can also be identified with the Riemann sphere, and so we have several different models: This function is univalent inside the unit disk D and maps it onto the exterior Ω = E e of the ellipse. The holomorphically reflected function, defined in terms of f by f It is easy to interpret (46) on the boundary: on taking ζ = e iϕ ∈ ∂D it gives the parametrization f (e iϕ ) = a cos ϕ − ib sin ϕ of Γ. The minus sign appears because decreasing ϕ is to correspond to positive orientation of ∂Ω.
As mentioned, the function f (ζ) is univalent in the unit disk, and it remains so up to the radius R = a+b a−b = a+b c > 1. The annulus is covered twice, and the circle |ζ| = R is mapped onto the focal segment [−c, c], which is also a branch cut for the inverse map. On this circle |ζ| = R we have We notice that, in the notations of (33), the zeros ω j and poles ξ j of f are given by As for the mapping properties we can write, symbolically, f (D e ) = 2 · E + 1 · E e , see Figure 5.
The points ±ζ int and ±ζ ext have no particular geometric meaning themselves, it is rather the image points ±z int = f (±ζ int ) and ±z ext = f (±ζ ext ) which can be interpreted geometrically. They are most easily obtained by solving directly the equation 1 + zS(z) = 0 (relevant in (24)), which corresponds to 1+f (ζ)f * (ζ) = 0. That equation becomes and the solutions are They have symmetries similar to those for ζ, for example z int · z ext = ±1, and are located along the imaginary axis in the order In this physical z-plane the quadrature identity (54) becomes in terms of the above data and with h considered as a function of z = f (ζ). The identities, (54) and (61), are somewhat more precise versions of formulas appearing in [23].
6 The hyperbola

General
We consider next a general hyperbola on the standard form Γ : obtained from the equation (43) for the ellipse by replacing b with ib. Thus most of the computations will be the same as those for the ellipse, but the geometric considerations and certain other issues will be different, and for that reason we prefer to be fairly explicit with the details. In addition we consider the hyperbola to be the main new object of investigation for the present paper. For the ellipse we assumed a > b, mainly because one usually visualizes the ellipse as oriented along the x-axis. For the hyperbola we assume only a, b > 0.
In complex coordinates the equation (62) becomes P (z,z) = 0, where Solving forz gives the Schwarz function where now c = √ a 2 + b 2 > 0. The two branch points ±c are located on the real axis, and branch cuts, not intersecting the curve, can be chosen to consist of the two segments [c, +∞), (−∞, −c], which meet at the point of infinity of the Riemann sphere. The above means that S(z) has two single-valued branches near the origin. That branch which satisfies S(0) = 2ab/c analytically continues to Γ and satisfies S(z) =z on Γ. So this branch is the actual Schwarz function of Γ. One might then expect that a null quadrature identity similar to that for the exterior of an ellipse, for example (18), would hold in that component of C \ Γ which is free from branch points. But it does not, and one naturally wonders why. In order to produce an answer we need to make a fairly systematic study of quadrature properties for hyperbola domains. We name the components of C \ Γ according to Thus S(z) has single-valued branches in Ω 0 , but not in Ω ± .
When trying to complete the above components by adding points of infinity one really has to make a choice in which space to make the completion: in CP 1 or in RP 2 red . In the first case there is only one point of infinity, and after an inversion the hyperbola becomes a lemniscate and looks like the figure "∞". Then Ω 0 is mapped to the exterior component, and the Ω ± are mapped to the insides of the two small loops (in "∞"). See further Section 8.2.
In the second case the curve Γ itself completes into the set of those [ t : z ] ∈ RP 2 red satisfying Choosing t = 1 gives the finite part of the hyperbola, namely (63), while the additional points of infinity are of the form [ 0 : This equation defines the asymptotic directions of the hyperbola. Note that the left member has period π, which means that the two directions of an asymptotic line correspond to the same point of infinity in the limit. In particular this means that the two components of Γ do not intersect at the projective infinity, but are pieces of one and the same smooth curve in RP 2 red . If the component of Γ which bounds Ω + is oriented so that Ω + lies to the left, then this curve continues with the other component of Γ oriented so that Ω 0 lies to the left, see Figure 3.
This orientation of Γ, in the direction of decreasing values of y, is consistent with the idea that Γ is the boundary of the multi-sheeted domain with an accompanying function f * (ζ) as in (47), namely With P (z, w) as in (63) with notations as in (33). Thus ω j ∈ ∂D, ξ j ∈ ∂U. The zeros ζ ± = ±(a + ib)/c of the derivative f are mapped by f onto the focal points ±c, which hence are branch points for f considered as a covering map.
Remark 6.1. One may wonder why the involution changes when one replaces b by ib, i.e. when moving from ellipse to hyperbola. With P (z, w) as in (45) and f and f * as in (46) and (47), the algebraic identity P (f, f * ) = 0 holds independent of what type of values a and b take. But when insisting on having the relationship f * = f • J the meaning of J has to change because of the action of complex conjugation.
As a further step one may take both of a and b to be imaginary, i.e. to change them, in (43), (45), (46), (47), to ia and ib respectively. Then P (f, f * ) = 0 still holds, but now we can interpret f * = f • J only with J(ζ) = −1/ζ. This involution has no fixed points, as is consistent with the fact that the curve (43) in this case has no real locus.

Quadrature data
According to (38) the quadrature nodes come from the solutions ζ ∈ U of 1 + f (ζ)f * (ζ) = 0. This equation spells out to giving in total four solutions in C, branch cut branch cut quadrature node (on lower sheet) Ω 0 Ω + Ω − lower sheet of quadrature surface, with continuation underneath Ω + upper sheet Figure 6: The hyperbola x 2 −y 2 = 1 as part of the Riemann sphere, with the single quadrature node (of order two) located on the curve itself, however on a diffeent sheet (as interior point of Ω 0 and its continuation under Ω + ).
The parameter value b = 1 marks a phase change, and we have to consider three different regimes: • 0 < b < 1: In this regime the right member of (69) is a negative real number for both choices of signs. Indeed It follows that all roots ζ j in (52) are purely imaginary, and that we can order them according to This means that ζ 1 and ζ 3 correspond to choosing the plus sign in (69), ζ 2 and ζ 4 to choosing the minus sign, and that the outer square root in (70) is chosen so that ζ 1 , ζ 2 ∈ U.
With the above agreements the identity (38) becomes, when b = 1, This is in principle of the form The quadrature identity associated to the lower half plane is of the same form: In similarity with the ellipse case, insertion of (69) in (72) gives only two possible values, A + and A − , of the coefficients A j . In the hyperbola case these need not be real, and to exhibit the situation clearly we have to separate the cases b < 1 and b > 1. After some computations one gets In both cases A + + A − = 1, and, for the difference, with positive square roots in both cases. A detailed analysis of the above formulas leads to the following identifications, for (73) and (74): As a consequence, the main quadrature identity (73) takes the form when b > 1. The first terms in the left member are quite natural since the partly double-sheeted image domain of U is to occupy as much area as the Riemann sphere itself (the same for the image domain of L). The second term is purely imaginary in the case b < 1, and a negative real number in the case b > 1. Note that it be made arbitrarily big by choosing b to be sufficiently close to one. Also, b > 1 can be adapted so that the contribution from h(ζ 1 ) disappears, see further Section 6.4. Turning to the special case b = 1, we have ζ 1 = ζ 2 = i and the quadrature identity becomes One can show this by letting ζ 1 , ζ 2 → i (whereby b → 1) in the right member of (73) and using (70), (77), (78), (79) when passing to the limit:

Quadrature identities in the physical plane
Returning to the general hyperbola, and considering it now as a partly double-sheeted domain in the z-plane, the quadrature nodes z j = f (ζ j ) in this picture are most easily obtained by solving the equation 1 + zS(z) = 0, namely (58) with b replaced by ib. Thus (compare (59), (60)), The analysis of the square roots is similar to what we had in the previous section: the outer square root selects between the two pairs {z 1 , z 2 } and {z 3 , z 4 } and the inner square root makes selections within each of them. The result is the following classification (compare Figure 7).
• 0 < b < 1: The nodes z j consist of two complex conjugate pairs located on the unit circle with z 1 , z 2 , z 3 , z 4 in, respectively, the fourth, first, second and third quadrant.
• b > 1: All nodes are real and satisfy The quadrature identity is in all cases of the form with A ± as in (75), (76). The integration takes place on a two-sheeted and branched Riemann surface and the function h is allowed to take different values on the two sheets, these being joined by branch points at ±c. When b = 1, with z 1 = z 2 = 1 and a derivative of h appearing in the quadrature identity, see (80), one has to take the chain rule into account when interpreting that derivative as a derivative with respect to z: since d dζ h(f (ζ)) = f (ζ)h (z) a factor f (i) = a has to be included. The result is Specializing further to a = b = 1, the unique quadrature node z 1 = z 2 = 1 for Ω is in this case located on the curve itself. The node belongs actually to the interior of that sheet of domain Ω to which z 1 = z 2 belongs, but happens to simultaneously be on the boundary of the other sheet. See Figure 6.

Loss of weight
Considering in some more detail the case b > 1, a perhaps astonishing phenomenon is that the coefficient for h(z 1 ) in (81) may vanish. In fact, given any a > 0, the coefficient A − vanishes if ab = (a 2 + 1)(b 2 − 1), which occurs for b = √ 1 + a 2 . See equation (76). The quadrature nodes are then, in the parameter space, given by c ζ j = ±i √ 1 ∓ 2iab, the order of the inner signs being coupled to those in A ± . The nodes in the z-plane become Thus the node z 1 = c looses its weight and becomes only "virtual". The quadrature identity reduces to  An important observation is that c is a branch point. In other words, f (ζ 1 ) = 0 in the above notation, and since df ∧ df /(1 + |f | 2 ) 2 contains the factor |f | 2 this makes functions h with a light singularity be integrable when lifted to the covering surface. Slightly singular functions could therefore in principle be accepted as test functions.
In that sense the loss of a quadrature node can be explained by shortage of test functions, since we actually allow only those which are holomorphic. The above phenomenon, loss of weight for quadrature nodes at branch points, has been discussed also in [19] (see in particular Section 8.4.1 there), and it is responsible for crucial properties of contractive zero divisors in Bergman space [25,26].

Non-real coefficients and dipole densities
When b > 1, the coefficients A j = A ± in (81) are real numbers and the quadrature identity holds for the real and imaginary parts of h independently of each other. This means that (81) actually holds for harmonic test functions h. However, when b < 1 the coefficients A ± are non-real and the identity (81) mixes the real and imaginary parts of h. One can still obtain a quadrature identity for harmonic functions, but then one must extend the setting and allow more general quadrature functionals in the right member. To make this precise we rewrite (81), when b < 1 and with z 1 = x 1 + iy 1 , as On the latter form there is no mixing between real and imaginary parts, hence the equation holds for harmonic functions h without involving their harmonic conjugates. A potential theoretic interpretation of the equation is that it expresses gravi-equivalence between the spherical area measure restricted to Ω and two point sources at z 1 and z 1 plus a line dipole density on the segment fromz 1 to z 1 . As b → 1, whereby y 1 → 0, the formula agrees in the limit with (82). One has to keep in mind that y 1 < 0, and for that reason the final integral has size ≈ −2y 1 · ∂h(x 1 ) ∂x 1 .

Summary
In summary, we have seen that multi-sheeted domains bounded by hyperbolas are indeed quadrature domains of order two for the spherical area measure. If one of the quadrature nodes happens to coincide with a branch point the order actually decreases to one. Replacing spherical area measure by the Euclidean, the mentioned domains become null quadrature domains (still multi-sheeted), as indicated in Remark 4.4. The explanation is that the quadrature nodes move to infinity and that the holomorphic functions h used as test functions have to vanish of order two at infinity to qualify as test functions. Hence the quadrature nodes become invisible. As for the forms of the quadrature identities we obtained the following identities in the parameter space: when b < 1, and when b > 1. The quadrature nodes ζ j are given by (70). When b = 1, then ζ 1 = ζ 2 = i, and the identity has the form (80), i.e.
In all above identities h is assumed to be holomorphic in U and integrable with respect to the measure in the left member.

The parabola
We consider a general parabola on the form As a complex curve in CP 2 with coordinates [ t : z : w ] it is and in real projective space RP 2 red , with coordinates [ t : x + iy ], 4atx = y 2 .  (22), see [18]. For the general hyperbola (62) the points of infinity become [ 0 : a − ib : a + ib ] and [ 0 : a + ib : a − ib ].

Solving (84) forz gives the Schwarz function
which has one branch point at z = a, the other being z = ∞. Let D denote that component of C \ Γ which contains the branch cut [a, +∞), and let Ω the complementary component in C.
It is known [37] that Ω is a null quadrature domain for planar area measure. The corresponding result is also true in higher dimensions and can then be proved by exhausting the paraboloid by ellipsoids and using the main result in [14]. This technique of exhaustions has been developed by H. Shahgholian [43], L. Karp [29] and others.
For the spherical measure, Ω becomes a two point quadrature domain. The real parabola Γ can be parametrized by z = This function is meromorphic on the Riemann sphere and is univalent in the upper half plane U, with f (U) = Ω. Taking M = C ∪ {∞}, M + = U, we are thus in a similar situation as for the hyperbola. The holomorphically reflected function with respect to J(ζ) =ζ is and the quadrature identity becomes, as in (38), The equation 1 + f f * = 0 giving the quadrature nodes spells out to be 1 + 16a 2 ζ 2 (ζ 2 + 1) = 0, with roots These numbers ζ j (j = 1, 2, 3, 4) are non-real and come in complex conjugate pairs. The quadrature nodes are those two, say ζ 1 and ζ 2 , which are in the upper half plane U. For the further analysis there are three cases: • 0 < a < 1 2 : There is one root in each quadrant. • a = 1 2 : There are two double roots on the imaginary axis, The double roots have been split into conjugate pairs of simple roots located on the imaginary axis.
In the case of simple roots ζ j ∈ U we have Res which in view of (88) makes the quadrature formula (87) explicit. In the z-plane one gets but simplifications of this expression seem difficult in general. Therefore we shall be more explicit only in the special case of double roots, namely when a = 1 2 . Then and the quadrature formula for Ω becomes, by straight-forward computations using (24), (86) directly, Note that the coefficient 1 − 1/(2 √ 2) is slightly smaller than one (= the area of the full sphere), as one might expect from the exterior Ω of the parabola occupying a major part of the sphere. See Figure 8 (44) gives or, on real form, See Figure 9 for the shape. It can be viewed as the "smash sum" of two partially overlapping disks, compare [10,35] for the planar case.  We see that the equation now has degree four, and it describes a classical two point quadrature domain, known as the Neumann oval or, according to some sources, hippopede, which is a special case of a hypocycloid. In the planar case (Euclidean metric), the quadrature identity for the domain Ω inside (90) is quite simple: For the spherical metric the formula will be as in (61) along with (55), but with the nodes ±z int moved by the map which performs the inversion. In other words, the quadrature nodes will be ±1/z int with z int as in (59). The curve (90) has degree d = 4, but is still uniformized by the Riemann sphere, hence it has genus zero. Therefore the genus formula (1) says that there must three singular points. Some of the singular points come form the quadrature nodes: in the general setting of (15) the extension of the algebraic curve ∂Ω to CP 2 has, at each of the two points of infinity, [ 0 : 1 : 0 ] and [ 0 : 0 : 1 ], m simple cusps of multiplicities n 1 , . . . , n m with distinct tangent directions. As the detailed analysis in [18] shows, this gives a total contribution n(n − 1) to the term for singular points in the genus formula (1). Here n = m k=1 m k is the order of the quadrature identity. In the present case with the hippopede, n = 2 and so n(n − 1) = 2. Thus there is one more singular point to account for. This is what is called a "special point" [44,45], that is an isolated solution of S(z) = z. It is located at [ 1 : 0 : 0 ] ∈ CP 2 , in other words at the origin in the complex plane.

The lemniscate
The antipodal version of the hyperbola is a lemniscate. Indeed, the inversion z → −1/z makes the equation (63) transform into or a Bernoulli lemniscates. Since inversions are rigid transformations of the Riemann sphere the lemniscate is a quadrature domain in the same sense as the hyperbola. The curve has degree d = 4 and has an obvious singularity at the origin (intersection between two smooth branches). The special case is a = b = 1 is conspicuous. The lemniscate then becomes z 2 +z 2 = 2z 2z2 , or which is the inversion of the hyperbola x 2 − y 2 = 1 in Figure 6. The quadrature identity takes the form (82) with the node located on the curve, see Figure 10. While hippopedes and cardioids (next section) are quadrature domains in the traditional sense, lemniscate domains really require the Riemann surface setting to qualify. In the same vein, lemniscates do not behave well with respect to planar Laplacian growth processes, see [33], one would need to consider such processes on covering surfaces, like in [19], in order for lemniscates to fit in.

The cardioid
On antipodal inversion of the parabola (84) one gets an unbounded cardioid. To keep the cardioid bounded one needs to invert with respect to a point not located on the parabola itself. The best choice is the focal point z = a, and thus setting w = (z − a) −1 one gets in the w-plane a bounded cardioid having a cusp at w = 0. It can be parametrized from the unit disk by and it is a well-studied quadrature domain. In the planar case (Euclidean area measure) there is a quadrature node of order two at w = −1/4a. The equation for the boundary becomes, with w = u + iv, again of degree d = 4. In similarity with the case of the hippopede, the double quadrature node accounts for two of the three singular points dictated by the genus formula. The remaining singular point is the cusp at the origin.
9 Appendix: Some computational details for the proof of Theorem 4.1 To prove the final statements of Theorem 4.1, namely to confirm the presence of cancellations leading to (38), we shall make a detailed local analysis of the differentials near an arbitrary point, which we take to be the origin in a local coordinate ζ. We assume the following expansions as ζ → 0, where a, b, c = 0 and k, , j are integers: f (ζ) = aζ k + higher powers, f (ζ) = bζ + higher powers, f * (ζ) = cζ j + higher powers.
such that the right member takes the simpler form Res h 1 + f f * · df f appearing in (38). For this purpose we need to identify the residue at ζ = 0 of the first differential in (94), which is There are three main cases, and some subcases of these.
• k + j > 0: In this case there is a residue if and only if r + − k = −1.
-If k = 0 then ≥ 0, hence r + − k ≥ 0, and there is no residue. -If k > 0 then r + − k = r − 1 and there is a residue if r = 0. The size of the residue is then b/a = k. However, this residue cancels with k terms named h(ω k ) (the latter "k" being just an index) in the first formula in the theorem. -If k < 0 then j > 0. This is excluded by assumption.
-If 1 + ac = 0, then the O(ζ)-term is not the main term, and we again get a residue contribution if an only if r + − k = −1. The rest of the analysis is exactly the same as in the case k + j > 0. -If 1 + ac = 0, then 1 + f f * has a zero at ζ = 0, of order s ≥ 1 say: This means that the differential in (95) looks like hdf f (1 + f f * ) = b aq · ζ r+ −k−s dζ + higher powers, and there is a residue contribution if r = s + k − − 1. This case represents the surviving residues in equation (36), namely those coming from points where 1 + f f * = 0.
To prove (38) we investigate the second differential in (94): There are again several cases.
• k + j > 0: In this case there is a residue contribution at ζ = 0 if and only if r + j + = −1. The size of the residue is bc.
-If k = 0, hence j = 0, then r + j + ≥ 0 and there is no residue. -If k = 0, so that r + j + = r + j + k − 1 = r − 1, there is a residue if r = 0. The size of the residue is bc = kac.
and there is a residue if r = s + k − − 1. This case represents the surviving residues in equation (38), namely those coming from points where 1 + f f * = 0.
• k + j < 0: Then the leading power in (96) is ζ r+ −k dζ, and there is a residue contribution if and only if r + − k = −1.
-If k = 0, then r + − k = r − 1 and there is a residue if r = 0.
* If k > 0 then j < 0. This is excluded by assumption. * If k < 0 there is a contribution, by bc/ac = b/a = k = −|k|. This residue cancels together with |k| terms named h(ζ j ) in the second formula of the theorem.
The above analysis confirms the quadrature identities (38) under the stated assumptions. In fact, it shows that even without such assumptions (of non-coinciding poles and zeros) the terms with explicit zeros and poles in the first formas (36) in the theorem are cancelled by terms in the residue sums. But in general it becomes difficult to specify the ranges of these sums.