Derivation of the cartesian formula for an ellipse

From Freepedia

The derivation of the cartesian form for an ellipse is simple and instructive. An ellipse is defined as a the loci of points equidistant to two fixed points called the foci. Assuming that the foci are located at (-c,0) and (c,0) (ie. the ellipse is centered at (0,0)) then the sum of the distance between any point (x,y) and the two foci is constant.

If (x,y) is any point on the ellipse and if <math>d_1</math> is the distance between (x,y) and (-c,0) and <math>d_2</math> is the distance between (x,y) and (c,0), i.e.

Image:Ellipse derivation 1.jpg

<math>d_1 = \sqrt {(x+c)^2+y^2}</math>

<math>d_2 = \sqrt {(x-c)^2+y^2}</math>

then

<math>d_1 + d_2 = 2a</math>

where a is the semimajor axis. From this we can derive the cartesian equation. Substituting:

<math>\sqrt {(x+c)^2+y^2} + \sqrt {(x-c)^2+y^2} = 2a</math>

To simplify we rearrange and square both sides.

<math>\sqrt {(x+c)^2+y^2} = 2a - \sqrt {(x-c)^2+y^2}</math>

<math>(x+c)^2 + y^2 = \left ( 2a - \sqrt{(x-c)^2+y^2} \right )^2</math>

<math>(x+c)^2 + y^2 = 4a^2 - 4a\sqrt{(x-c)^2+y^2} + (x-c)^2 +y^2</math>

Solving for the root and simplifying:

<math>\sqrt{(x-c)^2+y^2} = -{1 \over 4a} ((x+c)^2+y^2-4a^2-(x-c)^2-y^2) </math>

<math>\sqrt{(x-c)^2+y^2} = -{1 \over 4a} (x^2 + 2xc + c^2 -4a^2 -x^2 +2xc -c^2)</math>

<math>\sqrt{(x-c)^2+y^2} = -{1 \over 4a} (4xc - 4a^2)</math>

<math>\sqrt{(x-c)^2+y^2} = a - {c \over a}x</math>

A final squaring

<math>(x-c)^2+y^2 = a^2 - 2cx + {c^2 \over a^2}x^2</math>

<math>x^2 - 2xc + c^2 + y^2 = a^2 -2xc + {c^2 \over a^2}x^2</math>

<math>x^2 + c^2 + y^2 = a^2 + {c^2 \over a^2}x^2</math>

Grouping the x-terms and dividing with <math>a^2-c^2</math>

<math>x^2 \left( 1 - {c^2 \over a^2} \right) + y^2 = a^2 - c^2</math>

<math>x^2 \left( {a^2 - c^2 \over a^2} \right) + y^2 = a^2 - c^2</math>

<math>{x^2 \over a^2} + {y^2 \over a^2-c^2} = 1</math>

If x = 0 then <math>d_1 = d_2 = a = \sqrt {c^2+b^2}</math>

Therefore we can substitute

<math>b^2 = a^2-c^2</math>

And we have our desired equation:

<math>{x^2 \over a^2} + {y^2 \over b^2} = 1</math>