# Ellipsoid

Tri-axial ellipsoid with distinct semi-axis lengths ${\displaystyle c>b>a}$
Tri-axial ellipsoid with distinct semi-axes a, b and c
Ellipsoids of revolution (spheroid) with a pair of equal semi-axes (a) and a distinct third semi-axis (c) which is an axis of symmetry. The ellipsoid is prolate (top) or oblate (bottom) as c is greater than or less than a.
Circular sections of an ellipsoid

An ellipsoid is a quadric surface such that every planar cross section is either an ellipse, or is empty, or is reduced to a single point. This explains the name, which means "ellipse like". Equivalently, an ellipsoid is a quadric surface that is bounded, which means that it may be enclosed in a sufficiently large sphere.

An ellipsoid has three pairwise perpendicular axes of symmetry which intersect at a center of symmetry, called the center of the ellipsoid. The line segments that are delimited on the axes of symmetry by the ellipsoid are called the principal axes, or simply axes of the ellipsoid. If the three axes have different lengths, the ellipsoid is said to be tri-axial or (rarely) scalene, and the axes are uniquely defined.

If two of the axes have the same length, then the ellipsoid is an "ellipsoid of revolution", also called a spheroid. In this case, the ellipsoid is invariant under a rotation around the third axis, and there are thus infinitely many ways of choosing the two perpendicular axes of the same length. If the third axis is shorter, the ellipsoid is an oblate spheroid, if it is longer, it is prolate spheroid. If the three axes have the same length, the ellipsoid is a sphere.

## Standard equation

When choosing a Cartesian coordinate system, such that the origin is the center of the ellipsoid, and the coordinate axes are axes of the ellipsoid, the implicit equation of the ellipsoid has the standard form

${\displaystyle {x^{2} \over a^{2}}+{y^{2} \over b^{2}}+{z^{2} \over c^{2}}=1,}$

where a, b, c are positive real numbers.

The points (a, 0, 0), (0, b, 0) and (0, 0, c) lie on the surface. The line segments from the origin to these points are called the semi-principal axes of the ellipsoid, because a, b, c are half the length of the principal axes. They correspond to the semi-major axis and semi-minor axis of an ellipse.

One can assume without loss of generality that ${\displaystyle a\geq b\geq c}$. If ${\displaystyle a=b>c,}$ one has an oblate spheroid; if ${\displaystyle a>b=c,}$ one has a prolate spheroid; if ${\displaystyle a=b=c,}$ one has a sphere.

## Parameterization

The ellipsoid may be parameterized in several ways, which are simpler to express when the ellipsoid axes coincide with coordinate axes. A common choice is

{\displaystyle {\begin{aligned}x&=a\cos(\theta )\cos(\varphi ),\\y&=b\cos(\theta )\sin(\varphi ),\\z&=c\sin(\theta ),\end{aligned}}\,\!}

where

${\displaystyle -{\frac {\pi }{2}}\leq \theta \leq {\frac {\pi }{2}},\qquad -\pi \leq \varphi \leq \pi .}$

These parameters may be interpreted as spherical coordinates. More precisely, ${\displaystyle \pi /2-\theta }$ is the polar angle, and ${\displaystyle \varphi }$ is the azimuth angle of the point (x, y, z) of the ellipsoid.

## In general position

More generally, an arbitrarily oriented ellipsoid, centered at v, is defined by the solutions x to the equation

${\displaystyle (\mathbf {x-v} )^{\mathrm {T} }\!A\,(\mathbf {x-v} )=1,}$

where A is a positive definite matrix and x, v are vectors.

The eigenvectors of A define the principal axes of the ellipsoid and the eigenvalues of A are the reciprocals of the squares of the semi-axes: ${\displaystyle a^{-2}}$, ${\displaystyle b^{-2}}$ and ${\displaystyle c^{-2}}$.[1] An invertible linear transformation applied to a sphere produces an ellipsoid, which can be brought into the above standard form by a suitable rotation, a consequence of the polar decomposition (also, see spectral theorem). If the linear transformation is represented by a symmetric 3-by-3 matrix, then the eigenvectors of the matrix are orthogonal (due to the spectral theorem) and represent the directions of the axes of the ellipsoid: the lengths of the semiaxes are given by the eigenvalues. The singular value decomposition and polar decomposition are matrix decompositions closely related to these geometric observations.

### Parametric representation

ellipsoid as an affine image of the unit sphere

The key to a parametric representation of an ellipsoid in general position is the alternative definition:

• An ellipsoid is the affine image of the unit sphere.

An affine transformation can be represented by a translation with a vector ${\displaystyle {\vec {f}}_{0}}$ and a regular 3×3-matrix ${\displaystyle A}$:

${\displaystyle {\vec {x}}\mapsto {\vec {f}}_{0}+A{\vec {x}}={\vec {f}}_{0}+x{\vec {f}}_{1}+y{\vec {f}}_{2}+z{\vec {f}}_{3}}$,

where ${\displaystyle {\vec {f}}_{1},{\vec {f}}_{2},{\vec {f}}_{3}}$ are the column vectors of matrix ${\displaystyle A}$.

A parametric representation of an ellipsoid in general position can be obtained by the parametric representation of a unit sphere (see above) and an affine transformation:

• ${\displaystyle {\vec {x}}(\theta ,\varphi )={\vec {f}}_{0}+{\vec {f}}_{1}\cos \theta \cos \varphi +{\vec {f}}_{2}\cos \theta \sin \varphi +{\vec {f}}_{3}\sin \theta ,\quad -\pi /2\leq \theta \leq \pi /2,\ 0\leq \varphi <2\pi }$.

If the vectors ${\displaystyle {\vec {f}}_{1},{\vec {f}}_{2},{\vec {f}}_{3}}$ form an orthogonal system, the points with vectors ${\displaystyle {\vec {f}}_{0}\pm {\vec {f}}_{i},\ i=1,2,3,}$ are the vertices of the ellipsoid and ${\displaystyle |{\vec {f}}_{1}|,|{\vec {f}}_{2}|,|{\vec {f}}_{3}|}$ are the semi pricipal axes.

A surface normal vector at point ${\displaystyle {\vec {x}}(\theta ,\varphi )}$ is

${\displaystyle {\vec {n}}(\theta ,\varphi )={\vec {f}}_{2}\times {\vec {f}}_{3}\cos \theta \cos \varphi +{\vec {f}}_{3}\times {\vec {f}}_{1}\cos \theta \sin \varphi +{\vec {f}}_{1}\times {\vec {f}}_{2}\sin \theta .}$

For any ellipsoid there exists an implicit representation ${\displaystyle F(x,y,z)=0}$. If for simplicity the center of the ellipsoid is the origin, i.e. ${\displaystyle {\vec {f}}_{0}=(0,0,0)^{T}}$, the following equation describes the ellipsoid above [2]:

${\displaystyle F(x,y,z)=\operatorname {det} ({\vec {x}},{\vec {f}}_{2},{\vec {f}}_{3})^{2}\;+\;\operatorname {det} ({\vec {f}}_{1},{\vec {x}},{\vec {f}}_{3})^{2}\;+\;\operatorname {det} ({\vec {f}}_{1},{\vec {f}}_{2},{\vec {x}})^{2}\;-\;\operatorname {det} ({\vec {f}}_{1},{\vec {f}}_{2},{\vec {f}}_{3})^{2}=0}$

## Volume and surface area

### Volume

The volume bounded by the ellipsoid is

${\displaystyle V={\frac {4}{3}}\pi abc.}$

Note that this equation reduces to that of the volume of a sphere when all three elliptic radii are equal, and to that of an oblate or prolate spheroid when two of them are equal.

The volume of an ellipsoid is ${\displaystyle {\frac {2}{3}}}$ the volume of a circumscribed elliptic cylinder, and ${\displaystyle {\frac {\pi }{6}}}$ the volume of the circumscribed box.

The volumes of the inscribed and circumscribed boxes are respectively:

${\displaystyle V_{\text{inscribed}}={\frac {8}{3{\sqrt {3}}}}abc,\qquad V_{\text{circumscribed}}=8abc.}$

### Surface area

The surface area of a general (tri-axial) ellipsoid is[3][4]

${\displaystyle S=2\pi c^{2}+{\frac {2\pi ab}{\sin(\varphi )}}\left(E(\varphi ,k)\,\sin ^{2}(\varphi )+F(\varphi ,k)\,\cos ^{2}(\varphi )\right),}$
where
${\displaystyle \cos(\varphi )={\frac {c}{a}},\qquad k^{2}={\frac {a^{2}(b^{2}-c^{2})}{b^{2}(a^{2}-c^{2})}},\qquad a\geq b\geq c,}$

and where F(φ,k) and E(φ,k) are incomplete elliptic integrals of the first and second kind respectively.[1]

The surface area of an ellipsoid of revolution (or spheroid) may be expressed in terms of elementary functions:

${\displaystyle S_{\rm {oblate}}=2\pi a^{2}\left(1+{\frac {(1-e^{2})}{e}}\cdot {\mbox{artanh}}(e)\right)\quad {\mbox{where}}\quad e^{2}=1-{\frac {c^{2}}{a^{2}}}\quad (c
${\displaystyle S_{\rm {prolate}}=2\pi a^{2}\left(1+{\frac {c}{ae}}\cdot \arcsin(e)\right)\quad \qquad {\mbox{where}}\;\quad e^{2}=1-{\frac {a^{2}}{c^{2}}}\quad (c>a),}$

which, as follows from basic trigonometric identities, are equivalent expressions (i.e. the formula for ${\displaystyle S_{\rm {oblate}}}$ can be used to calculate the surface area of a prolate ellipsoid and vice versa). In both cases e may again be identified as the eccentricity of the ellipse formed by the cross section through the symmetry axis. (See ellipse). Derivations of these results may be found in standard sources, for example Mathworld.[5]

#### Approximate formula

${\displaystyle S\approx 4\pi {\sqrt[{p}]{\frac {a^{p}b^{p}+a^{p}c^{p}+b^{p}c^{p}}{3}}}.\,\!}$

Here p ≈ 1.6075 yields a relative error of at most 1.061%;[6] a value of p = 8/5 = 1.6 is optimal for nearly spherical ellipsoids, with a relative error of at most 1.178%.

In the "flat" limit of c much smaller than a, b, the area is approximately 2πab.

## Applications

The ellipsoidal shape finds many practical applications:

Geodesy
• Earth ellipsoid, a mathematical figure approximating the shape of the Earth.
• Reference ellipsoid, a mathematical figure approximating the shape planetary bodies in general.
Mechanics
Crystallography
• Index ellipsoid, a diagram of an ellipsoid that depicts the orientation and relative magnitude of refractive indices in a crystal.
• Thermal ellipsoid, ellipsoids used in crystallography to indicate the magnitudes and directions of the thermal vibration of atoms in crystal structures.
Lighting

### Dynamical properties

The mass of an ellipsoid of uniform density ρ is:

${\displaystyle m=\rho V=\rho {\frac {4}{3}}\pi abc.\,\!}$

The moments of inertia of an ellipsoid of uniform density are:

${\displaystyle I_{\mathrm {xx} }={\frac {m(b^{2}+c^{2})}{5}},\qquad I_{\mathrm {yy} }={\frac {m(c^{2}+a^{2})}{5}},\qquad I_{\mathrm {zz} }={\frac {m(a^{2}+b^{2})}{5}},}$
${\displaystyle I_{\mathrm {xy} }=I_{\mathrm {yz} }=I_{\mathrm {zx} }=0.\,\!}$

For a=b=c these moments of inertia reduce to those for a sphere of uniform density.

Artist's conception of Haumea, a Jacobi-ellipsoid dwarf planet, with its two moons

Ellipsoids and cuboids rotate stably along their major or minor axes, but not along their median axis. This can be seen experimentally by throwing an eraser with some spin. In addition, moment of inertia considerations mean that rotation along the major axis is more easily perturbed than rotation along the minor axis.[7]

One practical effect of this is that scalene astronomical bodies such as Haumea generally rotate along their minor axes (as does Earth, which is merely oblate); in addition, because of tidal locking, moons in synchronous orbit such as Mimas orbit with their major axis aligned radially to their planet.

A relaxed ellipsoid, that is, one in hydrostatic equilibrium, has an oblateness a − c directly proportional to its mean density and mean radius. Ellipsoids with a differentiated interior—that is, a denser core than mantle—have a lower oblateness than a homogeneous body. Over all, the ratio (bc)/(ac) is approximately 0.25, though this drops for rapidly rotating bodies.[8]

The terminology typically used for bodies physically rotating on their minor axis – not necessarily ellipsoids of revolution, in their mathematical construction – and whose shape is determined by their gravitational field is Maclaurin spheroid (oblate spheroid) and Jacobi ellipsoid (scalene ellipsoid). At faster rotations, piriform or oviform shapes can be expected, but these are not stable.

### Fluid properties

The ellipsoid is the most general shape for which it has been possible to calculate the creeping flow of fluid around the solid shape. The calculations include the force required to translate through a fluid and to rotate within it. Applications include determining the size and shape of large molecules, the sinking rate of small particles, and the swimming abilities of microorganisms.[9]

### In probability and statistics

The elliptical distributions, which generalize the multivariate normal distribution and are used in finance, can be defined in terms of their density functions. When they exist, the density functions f have the structure:

${\displaystyle f(x)=k\cdot g((x-\mu )\Sigma ^{-1}(x-\mu )^{\mathrm {T} })}$

where ${\displaystyle k}$ is a scale factor, ${\displaystyle x}$ is an ${\displaystyle n}$-dimensional random row vector with median vector ${\displaystyle \mu }$ (which is also the mean vector if the latter exists), ${\displaystyle \Sigma }$ is a positive definite matrix which is proportional to the covariance matrix if the latter exists, and ${\displaystyle g}$ is a function mapping from the non-negative reals to the non-negative reals giving a finite area under the curve.[10] The multivariate normal distribution is the special case in which ${\displaystyle g(z)=e^{-z/2}}$ for quadratic form ${\displaystyle z}$.

Thus the density function is a scalar-to-scalar transformation of a quadric expression. Moreover, the equation for any iso-density surface states that the quadric expression equals some constant specific to that value of the density, and the iso-density surface is an ellipsoid.

## In higher dimensions

The volume of a higher-dimensional ellipsoid (a hyperellipsoid) can be calculated using the dimensional constant given for the volume of a hypersphere. One can also define hyperellipsoids as the images of spheres under invertible linear transformations. The spectral theorem can again be used to obtain a standard equation akin to the one given above.