Since there does not appear to be a complete derivation of the surface area of an ellipsoid readily available on the web, consider the manual evaluation in all its lovely detail. The concept of volume element indicates that surface area can be calculated much like volume, as an integral over the square root of the determinant of the metric, and that is how the evaluation will proceed.
A typical parametrization of an ellipsoid is
Forming the space-like interval on this surface, the corresponding metric tensor is
which looks rather complicated, but the determinant simplifies significantly:
Evaluating an integral of the square root of this expression can be done with expansions of binomials. The overall factor of a sine in the square root of this determinant makes the choice natural for reducing factors in the integrand, and the other angular variable can be left in place. With the abbreviations
and taking advantage of symmetries in the integrand, the surface area integral is>
Given the generalized binomial expansion
as well as a trigonometric integral, the integration over φ is
As an aside, this is the expansion of a complete elliptic integral of the second kind. The evaluation of the surface area continues with further applications of the binomial expansion:
This looks rather hopeless until one realizes the sum over l does not touch the two parameter combinations and can be performed independently of them. The factors involved are
where the last line uses the Gauss summation formula. Putting this into the previous sum and reordering arguments slightly gives
which is a single Appell double hypergeometric function. This single Appell function form for the surface of an ellipsoid does not seem to be widely known. While a should be less than both b and c in accordance with the intermediate elliptic integral expansion, this form does not impose an ordering on b and c due to symmetry in the second and third parameters.
For a symmetric ellipsoid with the surface area reduces to a Gauss hypergeometric function that can be related to circular functions,
where the second step applies a contiguous relation and the third identifies two special cases. The two options for parameters are
using the fact that the inverse sine or inverse hyperbolic sine of an imaginary argument is the other function multiplied by the imaginary unit. For both expressions simplify to
since both the inverse sine and and inverse hyperbolic sine have a leading linear term in their expansions.
To relate the single Appell function form to Legendre elliptic integrals, recall the expansions near the beginning of this presentation:
The result for the surface area can be expressed as a linear superposition of the Appell functions appearing here, along with a third function to be determined:
Since the only one of these three Appell functions not symmetric in the second and third parameters is the third, its minus sign will be associated with the argument β for final consistency. The simplest way to establish the linear superposition appears to via an integral representation of the Appell function:
While this representation does not converge at the lower limit for the cases of parameters considered here, one can assume the infinities cancel from both sides of the equation and proceed anyway. With identifications
the linear superposition takes the form
Moving integrals to the left-hand side of the equation, the combined integrand is
Since derivatives of linear square roots have one less power algebraically, a combination to consider is
For the term inverse in x here to match that in the previous bracket, an overall factor of one half is needed and the other constants follow immediately:
The remaining function, ignoring the singular lower limit, is now simply
and the linear relation among the three Appell functions is
With the identifications
the surface area of the three-dimensional ellipsoid can be written
which is the expected result with a and c interchanged. This form reduces to the result for the prolate ellipsoid for , which for the traditional Legendre form means .
The single Appell function result above is patently simpler than the traditional form. It is equivalent to expressing the surface area in terms of a Carlson completely symmetric integral of the second kind.
Uploaded 2020.06.21 analyticphysics.com