Show that if $b\neq 0$, then the cross-product term can be eliminated from the quadratic $ax^2 + 2bxy + cy^2$ by rotating the coordinate axes through an angle $\theta$ that satisfies the equation $$ \cot{2\theta}=\frac{a-c}{2b}. $$
I know that $ax^2 + 2bxy + cy^2$ can be rewritten as $$ \begin{bmatrix}x&y\end{bmatrix} \begin{bmatrix}a&b\\b&c\end{bmatrix} \begin{bmatrix}x\\y\end{bmatrix}, $$ and I am familiar with how to rotate it using unitary matrices. However, it is very tedious to find the eigenvectorss of $$ \begin{bmatrix}a&b\\b&c\end{bmatrix}. $$ Is there a method for solving this problem that does not require finding the eigenvectors?
Working directly with the rotation substitution $x=cx'-sy',\ y=sx'+cy',$ where $c=\cos \theta,\ s=\sin \theta,$ we need only compute the resulting coefficient of $x'y'$ and set it to zero. From the form $ax^2+2bxy+cy^2$ this coefficient is $$a(-2sc)+2b(c^2-s^2)+c(2sc),$$ and then (assuming $b \neq 0$) we have $$\frac{a-c}{2b}=\frac{c^2-s^2}{2sc}=\frac{\cos 2\theta}{\sin 2\theta}=\cot(2 \theta),$$ the desired formula.
Note: maybe it's a bad choice to use $s,c$ for the sine and cosine of the rotation angle, since $c$ also appears in the form $ax^2+2bxy+cy^2.$ However I think it's clear which is which in the above so I'll leave it as it is.