1 Introduction
The isoperimetric problem is one of the oldest in geometry and it consists of finding a surface of minimum area that encloses a given volume . For the readers not familiar with this problem we briefly comment on a physical experiment to identify possible shapes of a soap bubble inside a box. See Fig. 1.
In Fig. 1 one must consider that tilting the box will not change the shape of the bubble, and gravity is negligible because the soap film weighs close to zero. On top of Fig. 1 we can easily see the three shapes that appear in a lab experiment: spherical, cylindrical and planar. Bottom right we have a surface named after the German mathematician Karl Hermann Amandus Schwarz (18431921) but this shape will never occur in the experiment as proved in [22].
Bottom left we see another surface, this one found by Herbert Blaine Lawson (1942). Though it has never appeared in the bubble experiment there is no mathematical proof that discards Lawson surface. Maybe this bubble could exist under special conditions but it has been an open question for two decades already. We shall resume it with details in Sect. 4.
The isoperimetric problem is particularly important in physics because of its strong relation with stability, and this also involves the study of phenomena in nonEuclidean spaces. For example, by considering homogeneous density we can work with either volume or mass . In [18] the mathematical physicist Roger Penrose (1931) combined several results and evidences regarding gravitational collapse in order to conjecture that, if is the total mass and the area of a black hole, then . In the theory of General Relativity this and other inequalities are called isoperimetric inequalities for black holes [12]. For this problem the nonEuclidean spaces of interest are the Schwarzschild and the ReissnerNordstrom, and here we cite [9] for the readers who want more details.
Differently from the example in Fig. 1, of course one cannot customize lab experiments with a nonEuclidean metric. But we can resort to computational simulations, and one of the mostly used softwares for this purpose is the Surface Evolver. Firstly introduced in 1989, now its most recent version is 2.70 [5] with several applications in many Areas of Knowledge like Aerodynamics [3], Fluid Dynamics [8], [10], [23], and Medicine [11]. It handles forces (contact, gravity, etc.), pressures, densities, dimensional spaces (including nonEuclidean), tracking of quantities and prescribed energies, among several other features.
In this paper we use Evolver to study the experiment of Fig. 1 in a compact hyperbolic cubic box obtained by a tesselation of the whole space. Here we represent the hyperbolic 3space by the unit ball in centred at the origin and endowed with the Poincaré metric. It is known that, up to isometries, there exists a unique tesselation of by nonideal cubes . This result is proved in Sect. 2, which also includes two weak assumptions on isoperimetric regions inside . The assumptions come from the fact that any torus , given by a group of translations in , will be too poor in space symmetries compared with a Euclidean torus , where is a group of translations in .
2 Preliminaries
This section is devoted to main theorems and propositions used throughout the text.
Proposition 1.
Up to isometries there exists a unique tesselation of by nonideal cubes . These cubes have dihedral angle , meaning that cubes meet at each edge.
Proof.
For the uniqueness pick a vertex of the tessellation incident to edges. Now consider the unit sphere of the tangent space . This space is isometric to , whose origin is now . Of the edges each pair that belongs to a face of the cubic tesselation will make this face intersect with . These intersections determine a triangulation of that is both equilateral and equiangular, hence one of only three possible types: tetrahedral (), octahedral (), and icosahedral (). However, any hyperbolic cube has an acute dihedral angle. Therefore, only the case is possible.
Now define a ¯cube of radius by placing vertices on the 12 geodesic rays with icosahedral symmetry emanating from . The dihedral angle can be computed in terms of hyperbolic trigonometry: it is a continuous and strictly decreasing function of . This function converges to for , and to for . The Intermediate Value Theorem gives a cube with dihedral angle . ∎
Lemma 1.
Let be the unit disk in centred at the origin. Inversion with respect to a circumference of radius centred in keeps if and only if .
Proof.
Fig. 2 illustrates in light grey and the circumference , . By taking as the complex plane we know that such inversion is
(1) 
Let us define a lattice of isometries in which keeps the cubic tessellation. Consider the three pairs of opposite faces of . For each such pair there is a mirror reflection , , in one of the faces, and there are also mirror reflections , , which flip only the pair of opposite faces (and fix the midpoint of ). The latters are represented by mirror reflections in the three coordinate planes in Fig. 3(a). From Lemma 2 the formers are represented by spherical inversions such as
where and . These are the values that keep the tiling depicted in Fig. 3(b). Now we compose to obtain a “translation”, namely an isometry of hyperbolic type. Finally we let and . Then is a fundamental domain for .
By means of the Alexandrov Reflection Principle [1] we can reduce our study of the isoperimetric problem to a symmetric eighth . In the Poincaré model if is as in Fig. 3(a) then . If denotes hyperbolic volume then and . These values were obtained through the Surface Evolver, whose details are given in Sect. 3.
The conformal metric in is given by
(2) 
where is the identity matrix. For instance, if we consider a Euclidean radius the area of and volume inside are given by
(3) 
respectively. Notice that both and are strictly increasing with . As mentioned at the Introduction, in this work we adopt two assumptions on isoperimetric regions in .
Assumption 1.
Any isoperimetric region in intersects the three coordinate planes.
As a matter of fact we believe that Assumption 1 is always true but were not able to prove this fact yet. Up to ambient isometries it is valid in Euclidean threedimensional boxes, as proved in [22]. For now here is an example that motivates Assumption 1:
Example 1.
The set encloses volume . Now let be a sphere centred at the upper right corner of the central square in Fig. 3(b), which in Euclidean coordinates corresponds to where
(4) 
If encloses the same volume then and therefore .
Now we present a strong result that will be widely used in our work. According to [24], from [2], [14], [16], [17] we have:
Theorem 1.
Suppose is compact and is either empty or piecewise smooth. Then for any there exists a compact domain such that minimizes area among regions of volume . Moreover, the boundary of any minimizing region is a smooth embedded surface with constant mean curvature and, if , then meets orthogonally.
Any that minimizes area under volume constraint has constant mean curvature (CMC). For instance, see [4] to check this wellknown property. In their turn CMC surfaces follow the maximum principle (see [13]), which is the key to conclude that whenever . The property of orthogonal intersections will be used extensively in this work. The following result is an important clue on Assumption 1:
Proposition 2.
Let be an isoperimetric region of and . Then cannot contain open subsets of only zero, one, or two incident faces of . In particular we have .
Proof.
Let us apply Theorem 1 to . By contradiction, consider a continuous family of hyperbolic isometries that leave invariant the zero, one, or two coordinate planes. Let be the identity and . Now if intersects only zero, one, or two incident faces, then we consider for the respective family the supremum of the values such that . Since keeps then is also minimizing. However, touches tangentially, and so violates . ∎
Because of Proposition 2 successive reflections of in the coordinate planes will result in a potential solution of the isoperimetric problem for in the torus (equivalently in its fundamental domain ).
Assumption 2.
Let be as in Proposition 2. Then its intersection with any coordinate plane is either empty or a twodimensional connected graph.
In this paper we analyse under Assumptions 1 and 2. A little reflection shows that for there are only four nonempty possible connected topological cases depicted in Fig. 4. Notice that Fig. 4(b) is a graph in both directions since our geometry is hyperbolic. Moreover, we shall see that the numerical candidates for this case are either unduloids or a pair of equidistant tori.
Not all combinations of the three axes are feasible. Table 1 shows the eight possible combinations and the name we attribute for each one. The empty graph is called “e”, the triad of letters corresponds to , and in this order, respectively. More details on Table 1 will be given in Sect. 4. There we explain why is called inverted Lawson, and why some cases like , and are not included. Notice that is always in alphabetical order because of congruence. For example, acb is congruent to .
The isoperimetric problem considers . We could also take the origin but for the sake of visibility some cases in Table 1 are studied for .
Combination  Surface 

sphere  
unduloid  
inverted Lawson  
Lawson  
pair of tori  
Schwarz 
3 Evolver Datafiles
In order to study each case numerically we make use of the Surface Evolver [6], [7]. For lengths and areas Evolver computes Riemannian values through the metric given in the datafile. However, volumes must be achieved as a declared quantity and the computation does not use the given metric any longer. Therefore, we implemented the submanifold quantity sbmvol in order to have
(5) 
where and are the hyperbolic and Euclidean elements of volume, respectively.
But Evolver only computes surface integrals
(line integrals are in fact sums of thin strips of surface). So we need to introduce a vector field
(q1,q2,q3) in such a way that its divergence is . By setting we may first compute(6) 
Hence, if we take q1, q2 as identically zero then q3 is given by
(7) 
In this way (q1,q2,q3) are not symmetric. A better choice is to compute q1 and q2 as done for q3, and then take one third of each expression. That is what we use in the datafiles.
Moreover, notice that Assumption 1 and Theorem 1 enable us to study the isoperimetric regions in such a way that the computation of the volume will not use . For instance, the face of on makes (7) vanish, and the analogous holds for the symmetric (q1,q2,q3) used in our datafiles. Hence we can work with , as described in Proposition 2.
Inside the sphere centred at the origin cannot have . Namely, for the first line of Table 1 our simulation takes . With Evolver 2.70 and Geomview 1.9.4 we obtained Fig. 6. It depicts the numerical surface, whereas Fig. 6(b) shows for (3) and Fig. 6(c) the corresponding , where is the difference between numerical and theoretical values of area.
In our simulation Evolver reached numerical values of volume that coincide with the theoretical ones up to the 8th decimal. Hence we consider them as identical. In Fig. 6(c) we have , and the abrupt change at just means that we refined the triangulation for that volume. This is to guarantee that , where is the number of triangles. The fact that can be seen in Fig. 6(b), namely lesser than , which means a great difference from the Euclidean case. Indeed, if that eighth of cube were Euclidean the extreme would have more than half of its volume.
Now back to Table 1 we consider . Numerically speaking this case is rather different from , in which the initial surface was just a tiny equilateral triangle with vertices on the coordinate axes. As usual, Evolver starts with a very simple polyhedral surface whose triangulation must be consistently refined, equalised and submitted to energy minimisation under geometrical and quantitative constraints. After some iterations the numerical surface may serve as evidence to help answer theoretical questions.
For we could try the initial surface depicted in Fig. 7(a). It consists of a rectangular blue face and a transparent triangular face on top. This one is used by (q1,q2,q3) to compute the hyperbolic volume but we ask Evolver to paint it in CLEAR for does not count it.
We end up with a numerical unduloid by starting from Fig. 7(a). But it presents two numerical problems: the convergence is slow and the resulting unduloid is just a local minimum of area under volume constraint. The convergence problem arises from the long thin shape and the initial contact angle of circa with the coordinate planes. It does converge to as expected from Theorem 1 but not as quickly as , since now we have that long thin face.
The nonisoperimetric unduloid arises from the metric (2), which makes the blue face in Fig. 7(a) broad on top and narrow on bottom. That is the feature we get for the resulting surface, which is just a local minimum under volume constraint. From the Euclidean perspective the isoperimetric unduloid must look broader on bottom than on top, which is the case if we start from Fig. 7(b). Notice that the initial surface is also orthogonal to , which speeds up convergence to a numerical answer compatible with Theorem 1.
These strategies may look rather tricky but when dealing with numerical optimisation either Evolver or any other software will strive for convergence around values that are potentially local minima. Global minima will be hardly found without applying theoretical framework.
That said, we generated unduloids of volume according to (3) for varying from to . This last one is depicted in Fig. 7(c). Finally, Fig. 7(d) shows both and for . There the and curves are red and blue, respectively (cf. Fig. 6(b)). Notice the turning point at .
Unduloids like in Fig. 7(c) are surfaces obtained by revolution of a stretch of hypercircle around . If we take , as the horizontal and vertical axes in Fig. 3(b), respectively, then hypercircles are all equated as
(8) 
in our case restricted to and in Euclidean values. This motivates the following definition:
Definition 1.
Consider the curve in for positive and . Then is the surface of revolution obtained by rotation of this curve around . We call the vertical unduloid with axis .
Back again to Table 1 we now take , a much easier case from the numerical pointofview. For the corresponding is one quarter of the area of the central square in Fig. 3(b). By means of polar coordinates we get
(9) 
where the value in (9) was computed via Cauchy’s Residue Theorem. Namely, we begin with and does not grow very much with . Of course, it is constant in the Euclidean case. In Fig. 8(a) we see for the extreme value , which already surpasses . The turning point between and is at , as shown in Fig. 8(b). It summarizes the three main cases for our discussions in the next section.
4 Results
We begin this section by recalling an important question rose in [24]: Could Lawson surface be a solution of the isoperimetric problem for the cubic lattice in ? More precisely, let us take the torus obtained by identifying the opposite faces of the unitary cube in . According to [24, p.11], with Evolver one can get a Lawson surface of area 1.017 that encloses the volume . Technically speaking, if the corresponding theoretical surface has in fact area slightly lesser than 1, then the answer will be yes. Fig. 9(a) summarizes the corresponding three cases for a cubic lattice of . Lawson hyperbolic surface is depicted in Fig. 9(b) with for , already above . For reasons that will soon be explained the minimum of the three graphs in Fig. 8(b) will be called the isopcurve. This curve is compared with the graph in magenta of our simulation of Lawson’s case (see Fig. 9(c)).
In Fig. 9(a) the turning points occur at circa and , respectively. Since our first and second turning points occur for the ratios and , respectively. Roughly speaking, they happen earlier than for the Euclidean case. This explains why the curve is visibly above the isopcurve in Fig. 9(c).
We used (3) to make the curve start at for . In Fig. 7(d) the unduloid started at for . In theory an unduloid can have arbitrarily small volume but numerical simulations will only be meaningful if we avoid extreme cases. With Lawson surface there is an additional problem: it simply does not exist when is too little. In our simulations, for slightly below () the initial surface converges to a degenerate surface consisting of an eighth of connected to and by tiny tubes.
The analogous problem occurs with . Topologically speaking, this case is obtained through a Euclidean reflection of bbd in the plane determined by the points , , and , where is defined in (4). This plane works as a mirror that inverts the way we see Lawson surface, hence the name inverted Lawson.
In the Euclidean cube of edgelength the inverted Lawson is again Lawson surface but translated by . However, in the hyperbolic cube we get a different graph for the inverted surface. This is due to the metric (2), which requires more area to comprise the same volume when leaves the origin (see Example 1 in Sect. 2). Fig. 10(a) shows an inverted Lawson close to the extreme case of collapsing its handle around with the origin itself. Fig. 10(b) shows the graph for in blue, together with the other graphs in Fig. 9(c).
Since we have just recalled Example 1, notice that its arguments also work for unduloids:
Example 2.
Take instead of , respectively, where is given by Definition 1 and is the “inverted ”, as explained right above. Namely, and therefore . In the case of a horizontal unduloid, namely with axis either or , area and volume remain unchanged for the “inverted” surface.
In Example 2 the assertion about horizontal unduloids comes from the fact that their “inversion” can be equated. For the axis it is the composition of two isometries in : rotation followed by the spherical inversion
For the axis it is followed by
(11) 
Because of Examples 1 and 2 we had already expected the graph to lie above the graph as depicted in Fig. 10(b). For the sake of concision we could then have skipped the numerical analysis of . Nevertheless, it is important to include it here for two reasons: it justifies why Table 1 omits even some feasible cases (like and ), and it helps check the reliability of our numerical simulations as we are going to do right now.
Indeed, above was already expected but we can look even closer. The factor is due to the fivefold symmetry on the corners instead of the fourfold symmetry at the centre of any square in Fig. 3(b). In the Euclidean case any 3Dmanifold whose dimensions increase by the same factor will get increased by and , respectively. For a cylinder of constant height we have the same factor for both if top and bottom do not count. Roughly speaking Lawson surface is close to a pair of horizontal unduloids connected by a very small piece of sphere. By means of Taylor expansion we rewrite (3) as
(12) 
Namely grows by when grows by , but its contribution to Lawson surface must be lesser than that. As we have already explained in Example 2, the pair remains unchanged for horizontal unduloids. Therefore must grow very little for the inverted Lawson, as we see in Fig. 10(b). The growth ratio varies from to according to our numerical tests.
A little reflection shows that is the “inversion” of bdb, this one congruent to Lawson surface, and now we have arguments to skip . Table 1 also omits because it is the “inverted” . One obtains by taking the “inversion” of with respect to the rectangle , , , , hence this one was omitted from Table 1 as well.
Now we are going to study , namely Schwarz surface in Table 1. For reasons that we have already explained for the case of Lawson surface, it is already expected that Schwarz surface does not exist for too little volume. Indeed, our numerical experiments show that the initial surface degenerates to a piece of sphere connected to , and by tiny tubes when is too little. We begin with for until when the volume surpasses . See Fig. 11(a) for an illustration and also Fig. 11(b) for its graph in red.
Now what is the inverted Schwarz? By looking at Figs. 4(a) and (c) we realise that it should be again , which was called inverted Lawson. However, by arguments already given in this section we know that “inversion” leads to a nonisoperimetric and for this reason we do not care about the apparent dubious meaning of . The inverted Schwarz can then be skipped, even if Fig. 4 leaves only as a way to codify it. For the same reason we shall not differ subcases when dubious triads appear again.
Back to Table 1 we now study the remaining cases , and . Topologically speaking Fig. 12(a) shows that these surfaces are Lawson’s example reflected by (11). But differently from the previous inverted cases now we must analyse these ones numerically. Fig. 12(b) compares with Fig 9(c).
By looking at Figs. 12(b) and 10(b) we realise that is quite similar to the inverted Lawson. But lies above almost as a parallel curve while rises more quickly and is closer to for small values of . Anyway, we still can argue that needs more area for the same volume of because the contact of with the origin is not as “large” as in the case of . Together with the metric (2) this explains why performs worse than .
Finally we are going to study , which is in fact a degenerate case. From Figs. 13(a) to (c) one promptly recognises that its topology was discarded as a possible CMC surface in a Euclidean threedimensional box, as proved in [22]. However, since we cannot adapt all of the arguments in [22] to the hyperbolic geometry, and similar CMC surfaces were already found in Euclidean boxes (see the gyroids of [15, Fig.2]), then we must include the numerical analysis of this case. The volume is fixed and the area decreases in Figs. 13(a) to (c). Namely, no CMC surface can be found in with that topology.
We conclude this section with a brief comment on the cases depicted in Fig. 5. One of them is combined with aa and presented in Fig. 14. It looks like a Lawson surface rotated by around . What is a congruence in the Euclidean case changes drastically in . Fig. 14 shows an example with and , namely well above the corresponding for the magenta curve in Fig. 9(c). The magenta curve starts at but now what happens for slightly below is that the surface degenerates to a pair of eighths of sphere centred at and and connected by a tiny tube.
A little reflection explains the bad performance of the cases in Fig. 5: they make have less contact around with the fundamental planes compared with any option in Table 1. Indeed, the example in Fig. 14 is above all others in that table if we look at in Figs. 11(b) and 12(b). That is why we have not included the cases of Fig. 5 in Table 1.
5 Conclusions
At the beginning of Sect. 4 we recalled a question rose in [24]. It remains open until the present day but the corresponding question in has great chance to be answered because of the evidence shown in Fig. 9(c). There the curve is at least 10% distant from the isopcurve, whereas in the Euclidean case it is only 1.7%. This is consistent with the fact that close objects in the Euclidean viewpoint turn out to be farther apart in the hyperbolic geometry. Curiously in both cases the closest approach happens at the turning point between the two surfaces of genus one.
The numerical gap of 10% makes it easier to prove that no Lawson surface can be isoperimetric in a hyperbolic space form: the fist step it to simplify the elliptic equations that define these surfaces by means of a theoretical and accurate evaluation. See some examples in [20], [21]. By the way, in [21] we get results whose numerical error is at most 0.05%. Optimistically speaking, the 1.7% could even be handled by our techniques.
Neither the hyperbolic nor the Euclidean geometry have positive sectional curvatures. Lawson surfaces have chance to be isoperimetric in some lens spaces, even though they are just the quotient of by a finite cyclic group, which however gives bidirectional translations. Indeed, in Sect. 2 could have been defined just as instead of . In this case in Fig. 6(a) will rise indefinitely in the vertical direction. The same holds for the examples discussed in [22], [24] but since these works inspired ours then we have considered a cubic tesselation.
Notice that such tesselation exists in as we did for but without a group that acts totally discontinuously. Hence we cannot have translations because of fixed points. In order to see this, first consider the tiling in Fig. 3(b) and the corresponding tiling of depicted in Fig. 15(a). Of course, compactness implies there are six spherical squares and the south pole is centred at the bottom square, the darkest one in Fig. 15(a). Any edge of the tiling belongs to a plane of given by , . Hence continuous displacements that keep the tiling must come from rotations about the fundamental axes, and any such rotation will have two fixed points.
For the corresponding cubic tesselation of consider as the north and south poles, respectively. The stereographic projection that keeps the equator will take its inside to the unit ball in , so that corresponds to the origin. This ball is shown by grid lines in Fig. 15(b), and without loss of generality our tesselation contains the spherical cube centred at and depicted there. It has dihedral angle , hence by spherical reflections we get eight cubes that cover the whole . Now take any spherical displacement that keeps the tesselation. The origin will slide along a fundamental axis but the spherical cubes that are crossed by the other axes will remain invariant. Hence we shall get fixed points in the threedimensional case as well.
As a matter of fact the unit ball can be viewed as a spherical cube of dihedral angle , so that together with its reflection in one gets a tesselation of with only two cubes. However, for a new space form we must take the quotient by the group (identity and antipodal maps). Hence the quotient will have a nonorientable boundary, which is the projective plane. However, the interest in the isoperimetric problem relies mostly on spaces whose geometry fits in the following sequence of generalisations: from Euclidean to spaces forms, from these to homogeneous spaces, and from these to Lie groups in general.
This has to do with the amount of results available for these geometries. For instance, in homogeneous spaces we still count on the maximum principle for CMC surfaces, and on the fact that they are minima of area under volume constraint. The Surface Evolver works in any finite dimension and it displays threedimensional projections specified by the user. But of course our understanding and our control of the numerical computation will always work better in the third dimension. There the richest homogeneous spaces have fourdimensional isometry group, and for them much about CMC surfaces is already known. See [19] for a good start.
Acknowledgements
We thank Prof K.GroßeBrauckmann at the Maths Dept of the University of Darmstadt for his valuable help in this work. This was during his stay in Brazil in 2010, supported by São Paulo Research Foundation (FAPESP) proc.09/154080. The first author was partially supported by FAPESP proc.16/237466.
References
 [1] (1962) Uniqueness theorems for surfaces in the large. I. Amer. Math. Soc. Transl.(2) 21, pp. 341–354. Cited by: §2.
 [2] (1976) Existence and regularity almost everywhere of solutions to elliptic variational problems with constraints. Vol. 165, American Mathematical Soc.. Cited by: §2.
 [3] (2012) Deployment analysis of pneumatic envelopes including ascending balloons and inflatable aerodynamic decelerators. Journal of Spacecraft and Rockets 49 (2), pp. 413–421. Cited by: §1.
 [4] (1984) Stability of hypersurfaces with constant mean curvature. Math. Z. 185, pp. 339–353. Cited by: §2.
 [5] (2013) The Surface Evolver. Susquehanna University. External Links: Link Cited by: §1.
 [6] (1992) The Surface Evolver. Experimental Mathematics 1 (2), pp. 141–165. Cited by: §3.
 [7] (1994) Surface Evolver manual. Mathematics Department, Susquehanna Univerisity, Selinsgrove, PA 17870 (2.24), pp. 20. External Links: Link Cited by: §3.
 [8] (2016) Towards the computation of viscous flow resistance of a liquid bridge. International Journal of Computational Methods and Experimental Measurements 4 (1), pp. 42–49. Cited by: §1.
 [9] (2007) On isoperimetric surfaces in general relativity. Pacific Journal of Mathematics 231 (1), pp. 63–84. Cited by: §1.
 [10] (2016) Effects of ligaments shape on radiative heat transfer in metal foams. International Journal of Numerical Methods for Heat & Fluid Flow 26 (2). Cited by: §1.
 [11] (2015) A software tool based on the surface evolver for precise location of tumours as a preoperative procedure to partial mastectomy. In Journal of Physics: Conference Series, Vol. 633, pp. 012132. Cited by: §1.
 [12] (1972) The time symmetric initial value problem for black holes. Communications in Mathematical Physics 27 (2), pp. 87–102. Cited by: §1.

[13]
(2015)
Elliptic partial differential equations of second order
. Springer. Cited by: §2.  [14] (1983) On the regularity of boundaries of sets minimizing perimeter with a volume constraint. Indiana University Mathematics Journal 32 (1), pp. 25–37. Cited by: §2.
 [15] (1997) Gyroids of constant mean curvature. Experimental Mathematics 6 (1), pp. 33–50. Cited by: §4.
 [16] (1987) Boundary regularity for solutions of a partitioning problem. Archive for Rational Mechanics and Analysis 97 (3), pp. 261–270. Cited by: §2.
 [17] (2003) Regularity of isoperimetric hypersurfaces in riemannian manifolds. Transactions of the American Mathematical Society 355 (12), pp. 5041–5052. Cited by: §2.
 [18] (1973) Naked singularities. Annals of the New York Academy of Sciences 224 (1), pp. 125–134. Cited by: §1.
 [19] (2012) Constant mean curvature surfaces in homogeneous manifolds. Logos Verlag Berlin GmbH. Cited by: §5.
 [20] (2002) Theoretical evaluation of elliptic integrals based on computer graphics. Technical report External Links: Link Cited by: §5.
 [21] (2003) Accurate evaluation of elliptic integrals. Technical report External Links: Link Cited by: §5.
 [22] (1994) Superficies con curvatura media constante. PhD thesis, Universidad de Granada (ES) 1, pp. 1–106. External Links: Link Cited by: §1, §2, §4, §5.
 [23] (2016) Computational modelling of a spanning drop in a wedge with variable angle. International Journal of Undergraduate Research and Creative Activities 6, pp. 3. Cited by: §1.
 [24] (2001) The isoperimetric problem. In Global theory of minimal surfaces, Clay Math. Proc., Vol. 2, pp. 175–209. External Links: Link Cited by: Figure 1, §2, Figure 9, §4, §5, §5.