There is another way to find the that may strike you as funny. Use homogeneous coordinates for the points in the place so the exterior point is located at $P = (-1,2,1)$.
The homogeneous coordinates of a circle located at $(c_x,c_y)$ with radius $r$ are
$$ C = \begin{vmatrix} 1 & 0 & -c_x \\ 0 & 1 & -c_y \\ -c_x & -c_y & c_x^2+c_y^2-r^2 \end{vmatrix} $$
such that the expression $$\begin{pmatrix} x & y & 1 \end{pmatrix} \begin{vmatrix} 1 & 0 & -c_x \\ 0 & 1 & -c_y \\ -c_x & -c_y & c_x^2+c_y^2-r^2 \end{vmatrix} \begin{pmatrix} x \\ y \\ 1 \end{pmatrix} = 0 $$ simplifies to the standard equation for a circle
$$ (x-c_x)^2 + (y-c_y)^2 - r^2 = 0 $$
Ok, back to the problem. If you take the product $L=C P$, then $L$ is a line with coordinates $L=[a,b,c]$ such as $a x+b y + c =0$ is the equation of the line. This line connects the two tangent points.
In your case
$$ C = \begin{vmatrix} 1 & 0 & -4 \\ 0 & 1 & -2 \\ -4 & -2 & 16 \end{vmatrix} $$
so
$$ L = \begin{vmatrix} 1 & 0 & -4 \\ 0 & 1 & -2 \\ -4 & -2 & 16 \end{vmatrix} \begin{pmatrix} -1 \\ 2 \\1 \end{pmatrix} = \begin{bmatrix} -5 \\ 0 \\ 16 \end{bmatrix} $$
or $$ \left. (x,y,1) \cdot [-5,0,16] = 0 \right\} \left. 16-5 x =0 \right\} x = \frac{16}{5} $$
Plug the line equation in the circle equation to get the other coordinate
$$ \left. \left( \frac{16}{5} \right)^2 - 8 \left( \frac{16}{5} \right) + y^2 -4 y + 16 =0 \right\} y =2 \pm \frac{2\sqrt{21}}{5} $$
The coordinates are thus
$$ \begin{pmatrix} \frac{16}{5} \\ 2 - \frac{2\sqrt{21}}{5} \end{pmatrix} \mbox{ or } \begin{pmatrix} \frac{16}{5} \\ 2 + \frac{2\sqrt{21}}{5} \end{pmatrix}$$
NOTE: The relationship between the point $P$ and the line $L$ is that of pole and polar
:-)It's just a problem for myself -- not a school assignment or anything like that. P.S. (I'll go to sleep now and come back tomorrow.) – Svend Tveskæg Mar 25 '16 at 02:46