Matrix representation of conic sections

<h2 id="classification">Classification</h2>
Proper (non-degenerate) and <a href="/facts/Degenerate_conic/9Ek4ja1g">degenerate conic sections</a> can be distinguished<a class="footnote-ref" id="fnref:6" href="#fn:6">6</a><a class="footnote-ref" id="fnref:7" href="#fn:7">7</a> based on the <a href="/facts/Determinant/eGYqJ2TF">determinant</a> of 
 
 
 
 
 A
 
 Q
 
 
 =
 (
 A
 C
 −
 
 
 
 B
 
 2
 
 
 4
 
 
 )
 F
 +
 
 
 
 B
 D
 E
 −
 C
 
 
 D
 
 
 2
 
 
 −
 A
 
 
 E
 
 
 2
 
 
 
 4
 
 
 
 
 {\displaystyle A_{Q}=(AC-{\frac {B^{2}}{4}})F+{\frac {BDE-C{D}^{2}-A{E}^{2}}{4}}}
 
:
If 
 
 
 
 det
 
 A
 
 Q
 
 
 =
 0
 
 
 {\displaystyle \det A_{Q}=0}
 
, the conic is degenerate.
If 
 
 
 
 det
 
 A
 
 Q
 
 
 ≠
 0
 
 
 {\displaystyle \det A_{Q}\neq 0}
 
 so that 
 
 
 
 Q
 
 
 {\displaystyle Q}
 
 is not degenerate, we can see what type of conic section it is by computing the <a href="/facts/Minor_(mathematics)/CMudF5Eh">minor</a>, 
 
 
 
 det
 
 A
 
 33
 
 
 =
 A
 C
 −
 
 
 
 B
 
 2
 
 
 4
 
 
 
 
 {\displaystyle \det A_{33}=AC-{\frac {B^{2}}{4}}}
 
:

<ul><li>
 
 
 
 Q
 
 
 {\displaystyle Q}
 
 is a <a href="/facts/Hyperbola/6R6ER0mY">hyperbola</a> if and only if 
 
 
 
 det
 
 A
 
 33
 
 
 <
 0
 
 
 {\displaystyle \det A_{33}<0}
 
,</li>
<li>
 
 
 
 Q
 
 
 {\displaystyle Q}
 
 is a <a href="/facts/Parabola/Cweo3wgw">parabola</a> if and only if 
 
 
 
 det
 
 A
 
 33
 
 
 =
 0
 
 
 {\displaystyle \det A_{33}=0}
 
, and</li>
<li>
 
 
 
 Q
 
 
 {\displaystyle Q}
 
 is an <a href="/facts/Ellipse/1wUDnGxY">ellipse</a> if and only if 
 
 
 
 det
 
 A
 
 33
 
 
 >
 0
 
 
 {\displaystyle \det A_{33}>0}
 
.</li></ul>
In the case of an ellipse, we can distinguish the special case of a <a href="/facts/Circle/9fy5ggKn">circle</a> by comparing the last two diagonal elements corresponding to the coefficients of 
 
 
 
 
 x
 
 2
 
 
 
 
 {\displaystyle x^{2}}
 
 , 
 
 
 
 x
 y
 
 
 {\displaystyle xy}
 
 and 
 
 
 
 
 y
 
 2
 
 
 
 
 {\displaystyle y^{2}}
 
:

<ul><li>If 
 
 
 
 A
 =
 C
 
 
 {\displaystyle A=C}
 
 and 
 
 
 
 B
 =
 0
 
 
 {\displaystyle B=0}
 
, then 
 
 
 
 Q
 
 
 {\displaystyle Q}
 
 is a circle.</li></ul>
Moreover, in the case of a non-degenerate ellipse (with 
 
 
 
 det
 
 A
 
 33
 
 
 >
 0
 
 
 {\displaystyle \det A_{33}>0}
 
 and 
 
 
 
 det
 
 A
 
 Q
 
 
 ≠
 0
 
 
 {\displaystyle \det A_{Q}\neq 0}
 
), we have a <a href="/facts/Real_number/R02gw5Pb">real</a> ellipse if 
 
 
 
 (
 A
 +
 C
 )
 det
 
 A
 
 Q
 
 
 <
 0
 
 
 {\displaystyle (A+C)\det A_{Q}<0}
 
 but an <a href="/facts/Imaginary_number/6Xh9W3hq">imaginary</a> ellipse if 
 
 
 
 (
 A
 +
 C
 )
 det
 
 A
 
 Q
 
 
 >
 0
 
 
 {\displaystyle (A+C)\det A_{Q}>0}
 
. An example of the latter is 
 
 
 
 
 x
 
 2
 
 
 +
 
 y
 
 2
 
 
 +
 10
 =
 0
 
 
 {\displaystyle x^{2}+y^{2}+10=0}
 
, which has no real-valued solutions.
If the conic section is <a href="/facts/Degenerate_conic/9Ek4ja1g">degenerate</a> (
 
 
 
 det
 
 A
 
 Q
 
 
 =
 0
 
 
 {\displaystyle \det A_{Q}=0}
 
), 
 
 
 
 det
 
 A
 
 33
 
 
 
 
 {\displaystyle \det A_{33}}
 
 still allows us to distinguish its form:

<ul><li>Two intersecting lines (a hyperbola degenerated to its two asymptotes) if and only if 
 
 
 
 det
 
 A
 
 33
 
 
 <
 0
 
 
 {\displaystyle \det A_{33}<0}
 
.</li>
<li>Two parallel straight lines (a degenerate parabola) if and only if 
 
 
 
 det
 
 A
 
 33
 
 
 =
 0
 
 
 {\displaystyle \det A_{33}=0}
 
. These lines are distinct and real if 
 
 
 
 
 D
 
 2
 
 
 +
 
 E
 
 2
 
 
 >
 4
 (
 A
 +
 C
 )
 F
 
 
 {\displaystyle D^{2}+E^{2}>4(A+C)F}
 
, coincident if 
 
 
 
 
 D
 
 2
 
 
 +
 
 E
 
 2
 
 
 =
 4
 (
 A
 +
 C
 )
 F
 
 
 {\displaystyle D^{2}+E^{2}=4(A+C)F}
 
, and non-existent in the real plane if 
 
 
 
 
 D
 
 2
 
 
 +
 
 E
 
 2
 
 
 <
 4
 (
 A
 +
 C
 )
 F
 
 
 {\displaystyle D^{2}+E^{2}<4(A+C)F}
 
.</li>
<li>A single point (a degenerate ellipse) if and only if 
 
 
 
 det
 
 A
 
 33
 
 
 >
 0
 
 
 {\displaystyle \det A_{33}>0}
 
.</li></ul>
The case of coincident lines occurs if and only if the <a href="/facts/Rank_of_a_matrix/Yp9b5LK3">rank</a> of the 3 × 3 matrix 
 
 
 
 
 A
 
 Q
 
 
 
 
 {\displaystyle A_{Q}}
 
 is 1; in all other degenerate cases its rank is 2.<a class="footnote-ref" id="fnref:8" href="#fn:8">8</a>

<h2 id="central-conics">Central conics</h2>
When 
 
 
 
 det
 
 A
 
 33
 
 
 ≠
 0
 
 
 {\displaystyle \det A_{33}\neq 0}
 
 a geometric center of the conic section exists and such conic sections (ellipses and hyperbolas) are called central conics.<a class="footnote-ref" id="fnref:9" href="#fn:9">9</a>

<h3>Center</h3>
The center of a conic, if it exists, is a point that bisects all the chords of the conic that pass through it. This property can be used to calculate the coordinates of the center, which can be shown to be the point where the <a href="/facts/Gradient/TsR7uukj">gradient</a> of the quadratic function Q vanishes—that is,<a class="footnote-ref" id="fnref:10" href="#fn:10">10</a>

∇
        Q
        =
        
          [
          
            
              
                
                  ∂
                  Q
                
                
                  ∂
                  x
                
              
            
            ,
            
              
                
                  ∂
                  Q
                
                
                  ∂
                  y
                
              
            
          
          ]
        
        =
        [
        0
        ,
        0
        ]
        .
      
    
    {\displaystyle \nabla Q=\left[{\frac {\partial Q}{\partial x}},{\frac {\partial Q}{\partial y}}\right]=[0,0].}

This yields the center as given below.
An alternative approach that uses the matrix form of the quadratic equation is based on the fact that when the center is the origin of the coordinate system, there are no linear terms in the equation. Any translation to a coordinate origin (x0, y0), using x* = x – x0, y* = y − y0 gives rise to

 
 
 
 
 
 (
 
 
 
 
 x
 
 ∗
 
 
 +
 
 x
 
 0
 
 
 
 
 
 y
 
 ∗
 
 
 +
 
 y
 
 0
 
 
 
 
 
 )
 
 
 
 
 (
 
 
 
 A
 
 
 B
 
 /
 
 2
 
 
 
 
 B
 
 /
 
 2
 
 
 C
 
 
 
 )
 
 
 
 
 (
 
 
 
 
 x
 
 ∗
 
 
 +
 
 x
 
 0
 
 
 
 
 
 
 
 y
 
 ∗
 
 
 +
 
 y
 
 0
 
 
 
 
 
 )
 
 
 +
 
 (
 
 
 
 
 D
 
 
 E
 
 
 
 
 )
 
 
 (
 
 
 
 
 
 x
 
 ∗
 
 
 +
 
 x
 
 0
 
 
 
 
 
 
 
 y
 
 ∗
 
 
 +
 
 y
 
 0
 
 
 
 
 
 
 )
 
 +
 F
 =
 0.
 
 
 {\displaystyle {\begin{pmatrix}x^{*}+x_{0}&y^{*}+y_{0}\end{pmatrix}}{\begin{pmatrix}A&B/2\\B/2&C\end{pmatrix}}{\begin{pmatrix}x^{*}+x_{0}\\y^{*}+y_{0}\end{pmatrix}}+\left({\begin{matrix}D&E\end{matrix}}\right)\left({\begin{matrix}x^{*}+x_{0}\\y^{*}+y_{0}\end{matrix}}\right)+F=0.}

The condition for (x0, y0) to be the conic's center (xc, yc) is that the coefficients of the linear x* and y* terms, when this equation is multiplied out, are zero. This condition produces the coordinates of the center:

(
            
              
                
                  
                    x
                    
                      c
                    
                  
                
              
              
                
                  
                    y
                    
                      c
                    
                  
                
              
            
            )
          
        
        =
        
          
            
              (
              
                
                  
                    A
                  
                  
                    B
                    
                      /
                    
                    2
                  
                
                
                  
                    B
                    
                      /
                    
                    2
                  
                  
                    C
                  
                
              
              )
            
          
          
            
            −
            1
          
        
        
          
            (
            
              
                
                  −
                  D
                  
                    /
                  
                  2
                
              
              
                
                  −
                  E
                  
                    /
                  
                  2
                
              
            
            )
          
        
        =
        
          
            (
            
              
                
                  (
                  B
                  E
                  −
                  2
                  C
                  D
                  )
                  
                    /
                  
                  (
                  4
                  A
                  C
                  −
                  
                    B
                    
                      2
                    
                  
                  )
                
              
              
                
                  (
                  D
                  B
                  −
                  2
                  A
                  E
                  )
                  
                    /
                  
                  (
                  4
                  A
                  C
                  −
                  
                    B
                    
                      2
                    
                  
                  )
                
              
            
            )
          
        
        .
      
    
    {\displaystyle {\begin{pmatrix}x_{c}\\y_{c}\end{pmatrix}}={\begin{pmatrix}A&B/2\\B/2&C\end{pmatrix}}^{\!-1}{\begin{pmatrix}-D/2\\-E/2\end{pmatrix}}={\begin{pmatrix}(BE-2CD)/(4AC-B^{2})\\(DB-2AE)/(4AC-B^{2})\end{pmatrix}}.}

This calculation can also be accomplished by taking the first two rows of the associated 
matrix AQ, multiplying each by (x, y, 1)⊤ and setting both inner products equal to 0, obtaining the following system:

 
 
 
 
 
 
 
 A
 x
 +
 (
 B
 
 /
 
 2
 )
 y
 +
 D
 
 /
 
 2
 
 
 
 =
 0
 ,
 
 
 
 
 (
 B
 
 /
 
 2
 )
 x
 +
 C
 y
 +
 E
 
 /
 
 2
 
 
 
 =
 0.
 
 
 
 
 
 
 {\displaystyle {\begin{aligned}Ax+(B/2)y+D/2&=0,\\(B/2)x+Cy+E/2&=0.\end{aligned}}}

This yields the above center point.
In the case of a parabola, that is, when 4AC − B2 = 0, there is no center since the above denominators become zero (or, interpreted <a href="/facts/Projective_geometry/rHKSaEbx">projectively</a>, the center is on the <a href="/facts/Line_at_infinity/RFnjNo29">line at infinity</a>.)

<h4>Centered matrix equation</h4>
A central (non-parabola) conic 
 
 
 
 A
 
 x
 
 2
 
 
 +
 B
 x
 y
 +
 C
 
 y
 
 2
 
 
 +
 D
 x
 +
 E
 y
 +
 F
 =
 0
 
 
 {\displaystyle Ax^{2}+Bxy+Cy^{2}+Dx+Ey+F=0}
 
 can be rewritten in centered matrix form as

(
            
              
                
                  x
                  −
                  
                    x
                    
                      c
                    
                  
                
                
                  y
                  −
                  
                    y
                    
                      c
                    
                  
                
              
            
            )
          
        
        
          
            (
            
              
                
                  A
                
                
                  B
                  
                    /
                  
                  2
                
              
              
                
                  B
                  
                    /
                  
                  2
                
                
                  C
                
              
            
            )
          
        
        
          
            (
            
              
                
                  x
                  −
                  
                    x
                    
                      c
                    
                  
                
              
              
                
                  y
                  −
                  
                    y
                    
                      c
                    
                  
                
              
            
            )
          
        
        =
        K
        ,
      
    
    {\displaystyle {\begin{pmatrix}x-x_{c}&y-y_{c}\end{pmatrix}}{\begin{pmatrix}A&B/2\\B/2&C\end{pmatrix}}{\begin{pmatrix}x-x_{c}\\y-y_{c}\end{pmatrix}}=K,}

where

K
        =
        −
        
          
            
              det
              (
              
                A
                
                  Q
                
              
              )
            
            
              A
              C
              −
              (
              B
              
                /
              
              2
              
                )
                
                  2
                
              
            
          
        
        =
        −
        
          
            
              det
              (
              
                A
                
                  Q
                
              
              )
            
            
              det
              (
              
                A
                
                  33
                
              
              )
            
          
        
        .
      
    
    {\displaystyle K=-{\frac {\det(A_{Q})}{AC-(B/2)^{2}}}=-{\frac {\det(A_{Q})}{\det(A_{33})}}.}

Then for the ellipse case of AC > (B/2)2, the ellipse is real if the sign of K equals the sign of (A + C) (that is, the sign of each of A and C), imaginary if they have opposite signs, and a degenerate point ellipse if K = 0. In the hyperbola case of AC < (B/2)2, the hyperbola is degenerate if and only if K = 0.

<h3>Standard form of a central conic</h3>
Main articles: <a href="/facts/Conic_section/teqgLptX">Conic section § Standard forms in Cartesian coordinates</a>, and <a href="/facts/Conic_section/teqgLptX">Conic section § Conversion to canonical form</a>
The standard form of the equation of a central conic section is obtained when the conic section is translated and rotated so that its center lies at the center of the coordinate system and its axes coincide with the coordinate axes. This is equivalent to saying that the coordinate system's center is moved and the coordinate axes are rotated to satisfy these properties. In the diagram, the original xy-coordinate system with origin O is moved to the x'y'-coordinate system with origin O'.

The translation is by the vector 
 
 
 
 
 t
 
 =
 
 
 (
 
 
 
 
 x
 
 c
 
 
 
 
 
 
 
 y
 
 c
 
 
 
 
 
 )
 
 
 .
 
 
 {\displaystyle \mathbf {t} ={\begin{pmatrix}x_{c}\\y_{c}\end{pmatrix}}.}

The rotation by <a href="/facts/Angle/gFUTEQYq">angle</a> α can be carried out by <a href="/facts/Matrix_diagonalization/y3MkyTCy">diagonalizing</a> the matrix A33.
Thus, if 
 
 
 
 
 λ
 
 1
 
 
 
 
 {\displaystyle \lambda _{1}}
 
 and 
 
 
 
 
 λ
 
 2
 
 
 
 
 {\displaystyle \lambda _{2}}
 
 are the <a href="/facts/Eigenvalue/8TjEoT8u">eigenvalues</a> of the matrix A33, the centered equation can be rewritten in new variables x' and y' as<a class="footnote-ref" id="fnref:11" href="#fn:11">11</a>

 
 
 
 
 λ
 
 1
 
 
 
 x
 
 ′
 
 2
 
 
 
 +
 
 λ
 
 2
 
 
 
 y
 
 ′
 
 2
 
 
 
 =
 −
 
 
 
 det
 
 A
 
 Q
 
 
 
 
 det
 
 A
 
 33
 
 
 
 
 
 .
 
 
 {\displaystyle \lambda _{1}x'^{2}+\lambda _{2}y'^{2}=-{\frac {\det A_{Q}}{\det A_{33}}}.}

Dividing by 
 
 
 
 K
 =
 −
 
 
 
 det
 
 A
 
 Q
 
 
 
 
 det
 
 A
 
 33
 
 
 
 
 
 
 
 {\displaystyle K=-{\frac {\det A_{Q}}{\det A_{33}}}}
 
 we obtain a standard canonical form.
For example, for an ellipse this form is

x
                  ′
                
              
              
                2
              
            
            
              a
              
                2
              
            
          
        
        +
        
          
            
              
                
                  y
                  ′
                
              
              
                2
              
            
            
              b
              
                2
              
            
          
        
        =
        1.
      
    
    {\displaystyle {\frac {{x'}^{2}}{a^{2}}}+{\frac {{y'}^{2}}{b^{2}}}=1.}

From here we get a and b, the lengths of the semi-major and semi-minor axes in conventional notation.
For central conics, both eigenvalues are non-zero and the classification of the conic sections can be obtained by examining them.<a class="footnote-ref" id="fnref:12" href="#fn:12">12</a> 

<ul><li>If λ1 and λ2 have the same algebraic sign, then Q is a real ellipse, imaginary ellipse or real point if K has the same sign, has the opposite sign or is zero, respectively.</li>
<li>If λ1 and λ2 have opposite algebraic signs, then Q is a hyperbola or two intersecting lines depending on whether K is nonzero or zero, respectively.</li></ul>
<h3>Axes</h3>
By the <a href="/facts/Principal_axis_theorem/cr8VMued">principal axis theorem</a>, the two <a href="/facts/Eigenvectors/8TjEoT8u">eigenvectors</a> of the matrix of the quadratic form of a central conic section (ellipse or hyperbola) are <a href="/facts/Perpendicular/u0k8xqx4">perpendicular</a> (<a href="/facts/Orthogonality/JVIjYPog">orthogonal</a> to each other) and each is parallel to (in the same direction as) either the <a href="/facts/Major_axis/l7JES2wo">major or minor axis</a> of the conic. The eigenvector having the smallest eigenvalue (in <a href="/facts/Absolute_value/dwJBpEs5">absolute value</a>) corresponds to the major axis.<a class="footnote-ref" id="fnref:13" href="#fn:13">13</a>
Specifically, if a central conic section has center (xc, yc) and an eigenvector of A33 is given by v(v1, v2) then the principal axis (major or minor) corresponding to that eigenvector has equation,

x
              −
              
                x
                
                  c
                
              
            
            
              v
              
                1
              
            
          
        
        =
        
          
            
              y
              −
              
                y
                
                  c
                
              
            
            
              v
              
                2
              
            
          
        
        .
      
    
    {\displaystyle {\frac {x-x_{c}}{v_{1}}}={\frac {y-y_{c}}{v_{2}}}.}

<h3>Vertices</h3>
The <a href="/facts/Vertex_(curve)/FkoCoHhV">vertices</a> of a central conic can be determined by calculating the intersections of the conic and its axes — in other words, by solving the system consisting of the quadratic conic equation and the linear equation for alternately one or the other of the axes. Two or no vertices are obtained for each axis, since, in the case of the hyperbola, the minor axis does not intersect the hyperbola at a point with real coordinates. However, from the broader view of the <a href="/facts/Complex_plane/vE0QrxJp">complex plane</a>, the minor axis of an hyperbola does intersect the hyperbola, but at points with complex coordinates.<a class="footnote-ref" id="fnref:14" href="#fn:14">14</a>

<h2 id="poles-and-polars">Poles and polars</h2>
Main article: <a href="/facts/Pole_and_polar/cqg8T8qM">Pole and polar</a>
Using <a href="/facts/Homogeneous_coordinates/TG9ebKN4">homogeneous coordinates</a>,<a class="footnote-ref" id="fnref:15" href="#fn:15">15</a> the points<a class="footnote-ref" id="fnref:16" href="#fn:16">16</a>

p
        
        =
        
          
            (
            
              
                
                  
                    p
                    
                      0
                    
                  
                
              
              
                
                  
                    p
                    
                      1
                    
                  
                
              
              
                
                  
                    p
                    
                      2
                    
                  
                
              
            
            )
          
        
      
    
    {\displaystyle \mathbf {p} ={\begin{pmatrix}p_{0}\\p_{1}\\p_{2}\end{pmatrix}}}
  
 and 
  
    
      
        
          r
        
        =
        
          
            (
            
              
                
                  
                    r
                    
                      0
                    
                  
                
              
              
                
                  
                    r
                    
                      1
                    
                  
                
              
              
                
                  
                    r
                    
                      2
                    
                  
                
              
            
            )
          
        
      
    
    {\displaystyle \mathbf {r} ={\begin{pmatrix}r_{0}\\r_{1}\\r_{2}\end{pmatrix}}}

are conjugate with respect to the conic Q provided

p
          
          
            
              T
            
          
        
        
          A
          
            Q
          
        
        
          r
        
        =
        0.
      
    
    {\displaystyle \mathbf {p} ^{\mathsf {T}}A_{Q}\mathbf {r} =0.}

The conjugates of a fixed point p either form a line or consist of all the points in the plane of the conic. When the conjugates of p form a line, the line is called the polar of p and the point p is called the pole of the line, with respect to the conic. This relationship between points and lines is called a polarity.
If the conic is non-degenerate, the conjugates of a point always form a line and the polarity defined by the conic is a <a href="/facts/Bijection/j4gTuTmW">bijection</a> between the points and lines of the extended plane containing the conic (that is, the plane together with the <a href="/facts/Point_at_infinity/GAuS0wHQ">points</a> and <a href="/facts/Line_at_infinity/RFnjNo29">line at infinity</a>).
If the point p lies on the conic Q, the polar line of p is the tangent line to Q at p.
The equation, in homogeneous coordinates, of the polar line of the point p with respect to the non-degenerate conic Q is given by

p
          
          
            T
          
        
        
          A
          
            Q
          
        
        
          
            (
            
              
                
                  x
                
              
              
                
                  y
                
              
              
                
                  z
                
              
            
            )
          
        
        =
        0.
      
    
    {\displaystyle \mathbf {p} ^{T}A_{Q}{\begin{pmatrix}x\\y\\z\end{pmatrix}}=0.}

Just as p uniquely determines its polar line (with respect to a given conic), so each line determines a unique pole p. Furthermore, a point p is on a line L which is the polar of a point r, if and only if the polar of p passes through the point r (<a href="/facts/Philippe_de_La_Hire/q3WSd0eH">La Hire</a>'s theorem).<a class="footnote-ref" id="fnref:17" href="#fn:17">17</a> Thus, this relationship is an expression of geometric <a href="/facts/Duality_(projective_geometry)/dxrjhRbB">duality</a> between points and lines in the plane.
Several familiar concepts concerning conic sections are directly related to this polarity. The center of a non-degenerate conic can be identified as the pole of the line at infinity. A parabola, being tangent to the line at infinity, would have its center being a point on the line at infinity. Hyperbolas intersect the line at infinity in two distinct points and the polar lines of these points are the asymptotes of the hyperbola and are the tangent lines to the hyperbola at these points of infinity. Also, the polar line of a focus of the conic is its corresponding directrix.<a class="footnote-ref" id="fnref:18" href="#fn:18">18</a>

<h2 id="tangents">Tangents</h2>
Let line L be the polar line of point p with respect to the non-degenerate conic Q. By La Hire's theorem, every line passing through p has its pole on L. If L intersects Q in two points (the maximum possible) then the polars of those points are tangent lines that pass through p and such a point is called an exterior or outer point of Q. If L intersects Q in only one point, then it is a tangent line and p is the point of tangency. Finally, if L does not intersect Q then p has no tangent lines passing through it and it is called an interior or inner point.<a class="footnote-ref" id="fnref:19" href="#fn:19">19</a>
The equation of the tangent line (in homogeneous coordinates) at a point p on the non-degenerate conic Q is given by,

 
 
 
 
 
 p
 
 
 
 T
 
 
 
 
 A
 
 Q
 
 
 
 
 (
 
 
 
 x
 
 
 
 
 y
 
 
 
 
 z
 
 
 
 )
 
 
 =
 0.
 
 
 {\displaystyle \mathbf {p} ^{\mathsf {T}}A_{Q}{\begin{pmatrix}x\\y\\z\end{pmatrix}}=0.}

If p is an exterior point, first find the equation of its polar (the above equation) and then the intersections of that line with the conic, say at points s and t. The polars of s and t will be the tangents through p.
Using the theory of poles and polars, the problem of finding the four mutual tangents of two conics reduces to finding the <a href="/facts/Conic_section/teqgLptX">intersection of two conics</a>.

<h2 id="see-also">See also</h2>
<ul><li><a href="/facts/Conic_section/teqgLptX">Conic section § General Cartesian form</a></li>
<li><a href="/facts/Quadratic_form_(statistics)/qH3sK2VG">Quadratic form (statistics)</a></li></ul>
<h2 id="notes">Notes</h2>

<ul><li>Ayoub, A. B. (1993), "The central conic sections revisited", <a href="/facts/Mathematics_Magazine/A55MtD1G">Mathematics Magazine</a>, 66 (5): 322–325, <a href="/facts/Doi_(identifier)/muM9Etpq">doi</a>:<a href="https://doi.org/10.1080%2F0025570x.1993.11996157">10.1080/0025570x.1993.11996157</a></li>
<li>Brannan, David A.; Esplen, Matthew F.; Gray, Jeremy J. (1999), Geometry, Cambridge University Press, <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 978-0-521-59787-6</li>
<li>Lawrence, J. Dennis (1972), A Catalog of Special Plane Curves, Dover</li>
<li>Ostermann, Alexander; Wanner, Gerhard (2012), Geometry by its History, Springer, <a href="/facts/Doi_(identifier)/muM9Etpq">doi</a>:<a href="https://doi.org/10.1007%2F978-3-642-29163-0">10.1007/978-3-642-29163-0</a>, <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 978-3-642-29163-0</li>
<li>Pettofrezzo, Anthony (1978) [1966], <a href="https://archive.org/details/matricestransfor0000pett">Matrices and Transformations</a>, Dover, <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 978-0-486-63634-4</li>
<li>Spain, Barry (2007) [1957], Analytical Conics, Dover, <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 978-0-486-45773-4</li></ul>

<h2 id="references">References</h2>

<ol>
<li id="fn:1">Brannan, Esplen & Gray 1999, p. 30 - Brannan, David A.; Esplen, Matthew F.; Gray, Jeremy J. (1999), Geometry, Cambridge University Press, ISBN 978-0-521-59787-6 <a href="#fnref:1" class="footnote-back-ref">↩</a></li>
<li id="fn:2">Pettofrezzo 1978, p. 110 - Pettofrezzo, Anthony (1978) [1966], Matrices and Transformations, Dover, ISBN 978-0-486-63634-4 <a href="https://archive.org/details/matricestransfor0000pett" target="_blank">https://archive.org/details/matricestransfor0000pett</a> <a href="#fnref:2" class="footnote-back-ref">↩</a></li>
<li id="fn:3">Spain 2007, pp. 59–62 - Spain, Barry (2007) [1957], Analytical Conics, Dover, ISBN 978-0-486-45773-4 <a href="#fnref:3" class="footnote-back-ref">↩</a></li>
<li id="fn:4">It is also a matrix of a quadratic form, but this form has three variables and is 
 
 
 
 A
 
 x
 
 2
 
 
 +
 B
 x
 y
 +
 C
 
 y
 
 2
 
 
 +
 D
 x
 z
 +
 E
 y
 z
 +
 F
 
 z
 
 2
 
 
 
 
 {\displaystyle Ax^{2}+Bxy+Cy^{2}+Dxz+Eyz+Fz^{2}}
 
. <a href="#fnref:4" class="footnote-back-ref">↩</a></li>
<li id="fn:5">Spain 2007, pp. 59–62 - Spain, Barry (2007) [1957], Analytical Conics, Dover, ISBN 978-0-486-45773-4 <a href="#fnref:5" class="footnote-back-ref">↩</a></li>
<li id="fn:6">Lawrence 1972, p. 63 - Lawrence, J. Dennis (1972), A Catalog of Special Plane Curves, Dover <a href="#fnref:6" class="footnote-back-ref">↩</a></li>
<li id="fn:7">Spain 2007, p. 70 - Spain, Barry (2007) [1957], Analytical Conics, Dover, ISBN 978-0-486-45773-4 <a href="#fnref:7" class="footnote-back-ref">↩</a></li>
<li id="fn:8">Pettofrezzo 1978, p. 110 - Pettofrezzo, Anthony (1978) [1966], Matrices and Transformations, Dover, ISBN 978-0-486-63634-4 <a href="https://archive.org/details/matricestransfor0000pett" target="_blank">https://archive.org/details/matricestransfor0000pett</a> <a href="#fnref:8" class="footnote-back-ref">↩</a></li>
<li id="fn:9">Pettofrezzo 1978, p. 105 - Pettofrezzo, Anthony (1978) [1966], Matrices and Transformations, Dover, ISBN 978-0-486-63634-4 <a href="https://archive.org/details/matricestransfor0000pett" target="_blank">https://archive.org/details/matricestransfor0000pett</a> <a href="#fnref:9" class="footnote-back-ref">↩</a></li>
<li id="fn:10">Ayoub 1993, p. 322 - Ayoub, A. B. (1993), "The central conic sections revisited", Mathematics Magazine, 66 (5): 322–325, doi:10.1080/0025570x.1993.11996157 <a href="https://doi.org/10.1080%2F0025570x.1993.11996157" target="_blank">https://doi.org/10.1080%2F0025570x.1993.11996157</a> <a href="#fnref:10" class="footnote-back-ref">↩</a></li>
<li id="fn:11">Ayoub 1993, p. 324 - Ayoub, A. B. (1993), "The central conic sections revisited", Mathematics Magazine, 66 (5): 322–325, doi:10.1080/0025570x.1993.11996157 <a href="https://doi.org/10.1080%2F0025570x.1993.11996157" target="_blank">https://doi.org/10.1080%2F0025570x.1993.11996157</a> <a href="#fnref:11" class="footnote-back-ref">↩</a></li>
<li id="fn:12">Pettofrezzo 1978, p. 108 - Pettofrezzo, Anthony (1978) [1966], Matrices and Transformations, Dover, ISBN 978-0-486-63634-4 <a href="https://archive.org/details/matricestransfor0000pett" target="_blank">https://archive.org/details/matricestransfor0000pett</a> <a href="#fnref:12" class="footnote-back-ref">↩</a></li>
<li id="fn:13">Ostermann & Wanner 2012, p. 311 - Ostermann, Alexander; Wanner, Gerhard (2012), Geometry by its History, Springer, doi:10.1007/978-3-642-29163-0, ISBN 978-3-642-29163-0 <a href="https://doi.org/10.1007%2F978-3-642-29163-0" target="_blank">https://doi.org/10.1007%2F978-3-642-29163-0</a> <a href="#fnref:13" class="footnote-back-ref">↩</a></li>
<li id="fn:14">Kendig, Keith (2005), Conics, The Mathematical Association of America, pp. 89–102, ISBN 978-0-88385-335-1 <a href="978-0-88385-335-1" target="_blank">978-0-88385-335-1</a> <a href="#fnref:14" class="footnote-back-ref">↩</a></li>
<li id="fn:15">This permits the algebraic inclusion of infinite points and a line at infinity which are necessary to have for some of the following results <a href="#fnref:15" class="footnote-back-ref">↩</a></li>
<li id="fn:16">This section follows Fishback, W.T. (1969), Projective and Euclidean Geometry (2nd ed.), Wiley, pp. 167–172 <a href="#fnref:16" class="footnote-back-ref">↩</a></li>
<li id="fn:17">Brannan, Esplen & Gray 1999, p. 189 - Brannan, David A.; Esplen, Matthew F.; Gray, Jeremy J. (1999), Geometry, Cambridge University Press, ISBN 978-0-521-59787-6 <a href="#fnref:17" class="footnote-back-ref">↩</a></li>
<li id="fn:18">Akopyan, A.V.; Zaslavsky, A.A. (2007), Geometry of Conics, American Mathematical Society, p. 72, ISBN 978-0-8218-4323-9 <a href="978-0-8218-4323-9" target="_blank">978-0-8218-4323-9</a> <a href="#fnref:18" class="footnote-back-ref">↩</a></li>
<li id="fn:19">Interpreted in the complex plane such a point is on two complex tangent lines that meet Q in complex points. <a href="#fnref:19" class="footnote-back-ref">↩</a></li>
</ol>

Matrix representation of conic sections open-in-new

Matrix representation of conic sections