Ellipse

DOWNLOAD Mathematica Notebook EXPLORE THIS TOPIC IN the MathWorld Classroom
EllipseBipolarEllipse construction

An ellipse is a curve that is the locus of all points in the plane the sum of whose distances r_1 and r_2 from two fixed points F_1 and F_2 (the foci) separated by a distance of 2c is a given positive constant 2a (Hilbert and Cohn-Vossen 1999, p. 2). This results in the two-center bipolar coordinate equation

 r_1+r_2=2a,
(1)

where a is the semimajor axis and the origin of the coordinate system is at one of the foci. The corresponding parameter b is known as the semiminor axis.

The ellipse is a conic section and a Lissajous curve.

An ellipse can be specified in the Wolfram Language using Circle[{x, y}, {a, b}].

If the endpoints of a segment are moved along two intersecting lines, a fixed point on the segment (or on the line that prolongs it) describes an arc of an ellipse. This is known as the trammel construction of an ellipse (Eves 1965, p. 177).

Elliptic gears

It is possible to construct elliptical gears that rotate smoothly against one another (Brown 1871, pp. 14-15; Reuleaux and Kennedy 1876, p. 70; Clark and Downward 1930; KMODDL).

The ellipse was first studied by Menaechmus, investigated by Euclid, and named by Apollonius. The focus and conic section directrix of an ellipse were considered by Pappus. In 1602, Kepler believed that the orbit of Mars was oval; he later discovered that it was an ellipse with the Sun at one focus. In fact, Kepler introduced the word "focus" and published his discovery in 1609. In 1705 Halley showed that the comet now named after him moved in an elliptical orbit around the Sun (MacTutor Archive). An ellipse rotated about its minor axis gives an oblate spheroid, while an ellipse rotated about its major axis gives a prolate spheroid.

A ray of light passing through a focus will pass through the other focus after a single bounce (Hilbert and Cohn-Vossen 1999, p. 3). Reflections not passing through a focus will be tangent to a confocal hyperbola or ellipse, depending on whether the ray passes between the foci or not.

Let an ellipse lie along the x-axis and find the equation of the figure (1) where F_1 and F_2 are at (-c,0) and (c,0). In Cartesian coordinates,

 sqrt((x+c)^2+y^2)+sqrt((x-c)^2+y^2)=2a.
(2)

Bring the second term to the right side and square both sides,

 (x+c)^2+y^2=4a^2-4asqrt((x-c)^2+y^2)+(x-c)^2+y^2.
(3)

Now solve for the square root term and simplify

sqrt((x-c)^2+y^2)=-1/(4a)(x^2+2xc+c^2+y^2-4a^2-x^2+2xc-c^2-y^2)
(4)
=-1/(4a)(4xc-4a^2)
(5)
=a-c/ax.
(6)

Square one final time to clear the remaining square root,

 x^2-2xc+c^2+y^2=a^2-2cx+(c^2)/(a^2)x^2.
(7)

Grouping the x terms then gives

 x^2(a^2-c^2)/(a^2)+y^2=a^2-c^2,
(8)

which can be written in the simple form

 (x^2)/(a^2)+(y^2)/(a^2-c^2)=1.
(9)

Defining a new constant

 b^2=a^2-c^2
(10)

puts the equation in the particularly simple form

 (x^2)/(a^2)+(y^2)/(b^2)=1.
(11)

The parameter b is called the semiminor axis by analogy with the parameter a, which is called the semimajor axis (assuming b<a). The fact that b as defined above is actually the semiminor axis is easily shown by letting r_1 and r_2 be equal. Then two right triangles are produced, each with hypotenuse a, base c, and height b=sqrt(a^2-c^2). Since the largest distance along the minor axis will be achieved at this point, b is indeed the semiminor axis.

If, instead of being centered at (0, 0), the center of the ellipse is at (x_0, y_0), equation (◇) becomes

 ((x-x_0)^2)/(a^2)+((y-y_0)^2)/(b^2)=1.
(12)

As can be seen from the Cartesian equation for the ellipse, the curve can also be given by a simple parametric form analogous to that of a circle, but with the x and y coordinates having different scalings,

x=acost
(13)
y=bsint.
(14)

The general quadratic curve

 ax^2+2bxy+cy^2+2dx+2fy+g=0
(15)

is an ellipse when, after defining

Delta=|a b d; b c f; d f g|
(16)
J=|a b; b c|
(17)
I=a+c,
(18)

Delta!=0, J>0, and Delta/I<0. Also assume the ellipse is nondegenerate (i.e., it is not a circle, so a!=c, and we have already established is not a point, since J=ac-b^2!=0). In that case, the center of the ellipse (x_0,y_0) is given by

x_0=(cd-bf)/(b^2-ac)
(19)
y_0=(af-bd)/(b^2-ac),
(20)

the semi-axes lengths are

a^'=sqrt((2(af^2+cd^2+gb^2-2bdf-acg))/((b^2-ac)[sqrt((a-c)^2+4b^2)-(a+c)]))
(21)
b^'=sqrt((2(af^2+cd^2+gb^2-2bdf-acg))/((b^2-ac)[-sqrt((a-c)^2+4b^2)-(a+c)])).
(22)

and the counterclockwise angle of rotation from the x-axis to the major axis of the ellipse is

 phi={0   for b=0 and a<c; 1/2pi   for b=0 and a>c; 1/2cot^(-1)((a-c)/(2b))   for b!=0 and a<c; pi/2+1/2cot^(-1)((a-c)/(2b))   for b!=0 and a>c.
(23)
EllipseDirectrix

The ellipse can also be defined as the locus of points whose distance from the focus is proportional to the horizontal distance from a vertical line known as the conic section directrix, where the ratio is <1. Letting r be the ratio and d the distance from the center at which the directrix lies, then in order for this to be true, it must hold at the extremes of the major and minor axes, so

 r=(a-c)/(d-a)=(sqrt(b^2+c^2))/d.
(24)

Solving gives

d=(a^2)/(sqrt(a^2-b^2))=(a^2)/c
(25)
r=(sqrt(a^2-b^2))/a=c/a.
(26)

The focal parameter of the ellipse is

p=(b^2)/(sqrt(a^2-b^2))
(27)
=(a^2-c^2)/c
(28)
=(a(1-e^2))/e,
(29)

where e is a characteristic of the ellipse known as the eccentricity, to be defined shortly.

Ellipse4Points

An ellipse whose axes are parallel to the coordinate axes is uniquely determined by any four non-concyclic points on it, and the ellipse passing through the four points (x_1,y_1), (x_2,y_2), (x_3,y_3), and (x_4,y_4) has equation

 |x^2 y^2 x y 1; x_1^2 y_1^2 x_1 y_1 1; x_2^2 y_2^2 x_2 y_2 1; x_3^2 y_3^2 x_3 y_3 1; x_4^2 y_4^2 x_4 y_4 1|=0.
(30)

Let four points on an ellipse with axes parallel to the coordinate axes have angular coordinates t_i for i=1, 2, 3, and 4. Such points are concyclic when

 s_1s_2s_3+s_1s_2s_4+s_1s_3s_4+s_2s_3s_4-(s_1+s_2+s_3+s_4)=0,
(31)

where the intermediate variable s_i=tan(t_i/2) has been defined (Berger et al. 1984; Trott 2006, pp. 39-40). Rather surprisingly, this same relationship results after simplification of the above where s_i is now interpreted as s_i=sin(t_i). An equivalent, but more complicated, condition is given by

 (s_1^4+s_2^4+s_3^4+s_4^4)+(4s_2^2s_3^2s_1^2+s_2^2s_4^2s_1^2+s_3^2s_4^2s_1^2+s_2^2s_3^2s_4^2)-4(s_2s_3s_4s_1^3+s_2s_3s_4^3s_1+s_2s_3^3s_4s_1+s_2^3s_3s_4s_1)-2(s_2^2s_1^2+s_3^2s_1^2+s_4^2s_1^2+s_2^2s_3^2+s_2^2s_4^2-2s_3^2s_4^2)+8s_2s_3s_4s_1. 
=0
(32)

Like hyperbolas, noncircular ellipses have two distinct foci and two associated directrices, each conic section directrix being perpendicular to the line joining the two foci (Eves 1965, p. 275).

Define a new constant 0<=e<1 called the eccentricity (where e=0 is the case of a circle) to replace b

 e=sqrt(1-(b^2)/(a^2)),
(33)

from which it follows that

b=asqrt(1-e^2)
(34)
=sqrt(a^2-c^2)
(35)
c=sqrt(a^2-b^2)
(36)
=ae
(37)
e=sqrt(1-(b^2)/(a^2))
(38)
=c/a.
(39)

The eccentricity can therefore be interpreted as the position of the focus as a fraction of the semimajor axis.

EllipseFocus

If r and theta are measured from a focus F instead of from the center C (as they commonly are in orbital mechanics) then the equations of the ellipse are

x=c+rcostheta
(40)
y=rsintheta,
(41)

and (◇) becomes

 ((c+rcostheta)^2)/(a^2)+(r^2sin^2theta)/(b^2)=1.
(42)

Clearing the denominators gives

 b^2(c^2+2crcostheta+r^2cos^2theta)+a^2r^2sin^2theta=a^2b^2.
(43)

Substituting in sin^2theta=1-cos^2theta gives

 b^2c^2+2rcb^2costheta+b^2r^2cos^2theta+a^2r^2-a^2r^2cos^2theta=a^2b^2.
(44)

Plugging in to re-express b and c in terms of a and e,

 a^2(1-e^2)a^2e^2+2aea^2(1-e^2)rcostheta+a^2(1-e^2)r^2cos^2theta 
 +a^2r^2-a^2r^2cos^2theta=a^2[a^2(1-e^2)].
(45)

Dividing by -a^2 and simplifying gives

 -r^2+[ercostheta-a(1-e^2)]^2=0,
(46)

which can be solved for r to obtain

 r=+/-[ercostheta-a(1-e^2)].
(47)

The sign can be determined by requiring that r must be positive. When e=0, (47) becomes r=+/-(-a), but since a is always positive, we must take the negative sign, so (47) becomes

 r=a(1-e^2)-ercostheta
(48)
 r(1+ecostheta)=a(1-e^2)
(49)
 r=(a(1-e^2))/(1+ecostheta).
(50)
Ellipsex

The distance from a focus to a point with horizontal coordinate x (where the origin is taken to lie at the center of the ellipse) is found from

 costheta=(x-c)/r.
(51)

Plugging this into (50) yields

 r+e(x-c)=a(1-e^2)
(52)
 r=a(1-e^2)-e(x-c).
(53)

In pedal coordinates with the pedal point at the focus, the equation of the ellipse is

 (b^2)/(p^2)=(2a)/r-1.
(54)

The arc length of the ellipse is

s(t)=aE(t,e)
(55)
=aE(t,sqrt(1-(b^2)/(a^2)))
(56)
=bE(t,sqrt(1-(a^2)/(b^2))),
(57)

where E(t,e) is an incomplete elliptic integral of the second kind with elliptic modulus e (the eccentricity).

The relationship between the polar angle from the ellipse center theta and the parameter t follows from

 theta=tan^(-1)(y/x)=tan^(-1)(b/atant).
(58)
EllipseFunction

This function is illustrated above with theta shown as the solid curve and t as the dashed, with b/a=0.6. Care must be taken to make sure that the correct branch of the inverse tangent function is used. As can be seen, theta weaves back and forth around t, with crossings occurring at multiples of pi/2. The curvature and tangential angle of the ellipse are given by

kappa(t)=(ab)/((b^2cos^2t+a^2sin^2t)^(3/2))
(59)
phi(t)=tan^(-1)(a/btant).
(60)

The entire perimeter p of the ellipse is given by setting t=2pi (corresponding to theta=2pi), which is equivalent to four times the length of one of the ellipse's quadrants,

p=aE(2pi,e)
(61)
=4aE(1/2pi,e)
(62)
=4aE(e),
(63)

where E(e) is a complete elliptic integral of the second kind with elliptic modulus e (the eccentricity). The perimeter can be computed using the rapidly converging Gauss-Kummer series as

p=pi(a+b)sum_(n=0)^(infty)(1/2; n)^2h^(2n)
(64)
=pi(a+b)(1+1/4h^2+1/(64)h^4+1/(256)h^6+...)
(65)

(OEIS A056981 and A056982), where (n; k) is a binomial coefficient and

 h=(a-b)/(a+b).
(66)

This can also be written analytically as

p=pi(a+b)_2F_1(-1/2,-1/2;1;h^2)
(67)
=2(a+b)[2E(h)+2(h^2-1)K(h)],
(68)

where _2F_1(a,b;c;z) is a hypergeometric function, K(k) is a complete elliptic integral of the first kind.

EllipsePerimeter

Approximations to the perimeter include

p approx pisqrt(2(a^2+b^2))
(69)
 approx pi[3(a+b)-sqrt((a+3b)(3a+b))]
(70)
 approx pi(a+b)(1+(3h)/(10+sqrt(4-3h))),
(71)

where the last two are due to Ramanujan (1913-1914), and (71) has a relative error of ∼3·2^(-17)h^5 for small values of h. The error surfaces are illustrated above for these functions.

The maximum and minimum distances from the focus are called the apoapsis and periapsis, and are given by

r_+=r_(apoapsis)=a(1+e)
(72)
r_-=r_(periapsis)=a(1-e).
(73)

The area of an ellipse may be found by direct integration

A=int_(-a)^aint_(-bsqrt(a^2-x^2)/a)^(bsqrt(a^2-x^2)/a)dydx
(74)
=int_(-a)^a(2b)/asqrt(a^2-x^2)dx
(75)
=(2b)/a{1/2[xsqrt(a^2-x^2)+a^2sin^(-1)(x/(|a|))]}_(x=-a)^a
(76)
=ab[sin^(-1)1-sin^(-1)(-1)]
(77)
=ab[pi/2-(-pi/2)]
(78)
=piab.
(79)

The area can also be computed more simply by making the change of coordinates x^'=(b/a)x and y^'=y from the elliptical region R to the new region R^'. Then the equation becomes

 1/(a^2)(a/bx^')^2+(y^('2))/(b^2)=1,
(80)

or x^('2)+y^('2)=b^2, so R^' is a circle of radius b. Since

 (partialx)/(partialx^')=((partialx^')/(partialx))^(-1)=(b/a)^(-1)=a/b,
(81)

the Jacobian is

|(partial(x,y))/(partial(x^',y^'))|=|(partialx)/(partialx^') (partialy)/(partialx^'); (partialx)/(partialy^') (partialy)/(partialy^')|
(82)
=|a/b 0; 0 1|=a/b.
(83)

The area is therefore

intint_(R)dxdy=intint_(R^')|(partial(x,y))/(partial(x^',y^'))|dx^'dy^'
(84)
=a/bintint_(R^')dx^'dy^'
(85)
=a/b(pib^2)
(86)
=piab,
(87)

as before. The area of an arbitrary ellipse given by the quadratic equation

 ax^2+bxy+cy^2=1
(88)

is

 A=(2pi)/(sqrt(4ac-b^2)).
(89)

The area of an ellipse with semiaxes a and b with respect to a pedal point P is

 A=1/2pi(a^2+b^2+|OP|^2).
(90)
EllipseNormalTangent

The unit tangent vector of the ellipse so parameterized is

x_T(t)=-(asint)/(sqrt(b^2cos^2t+a^2sin^2t))
(91)
y_T(t)=(bcost)/(sqrt(b^2cos^2t+a^2sin^2t)).
(92)

A sequence of normal and tangent vectors are plotted above for the ellipse.

The locus of the apex of a variable cone containing an ellipse fixed in three-space is a hyperbola through the foci of the ellipse. In addition, the locus of the apex of a cone containing that hyperbola is the original ellipse. Furthermore, the eccentricities of the ellipse and hyperbola are reciprocals. The locus of centers of a Pappus chain of circles is an ellipse. Surprisingly, the locus of the end of a garage door mounted on rollers along a vertical track but extending beyond the track is a quadrant of an ellipse (Wells 1991, p. 66). (The envelope of the door's positions is an astroid.)

Wolfram Web Resources

Mathematica »

The #1 tool for creating Demonstrations and anything technical.

Wolfram|Alpha »

Explore anything with the first computational knowledge engine.

Wolfram Demonstrations Project »

Explore thousands of free applications across science, mathematics, engineering, technology, business, art, finance, social sciences, and more.

Computerbasedmath.org »

Join the initiative for modernizing math education.

Online Integral Calculator »

Solve integrals with Wolfram|Alpha.

Step-by-step Solutions »

Walk through homework problems step-by-step from beginning to end. Hints help you try the next step on your own.

Wolfram Problem Generator »

Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.

Wolfram Education Portal »

Collection of teaching and learning tools built by Wolfram education experts: dynamic textbook, lesson plans, widgets, interactive Demonstrations, and more.

Wolfram Language »

Knowledge-based programming for everyone.