Reduced row echelon form of matrix with trigonometric expressions I'm trying to solve for the eigenvalues of and the eigenvectors of a rotation matrix (about the z-axis):
$$A = \begin{bmatrix}
    \cos\theta & -\sin\theta & 0 \\
    \sin\theta & \cos\theta & 0 \\
    0 & 0 & 1
\end{bmatrix}$$
Using $\det (\lambda I - A) = 0$, I determined that $\lambda \in \{1, \cos\theta\}$ (though I'm not certain that that is 100% correct). To solve for the eigenvectors, I decided that it would be easiest to put the following matrix in reduced row echelon form.
$$\begin{bmatrix}
    1 - \cos\theta & \sin\theta & 0 \\
    -\sin\theta & 1 - \cos\theta & 0 \\
    0 & 0 & 0
\end{bmatrix}$$
My TI-89 suggests that this matrix in reduced row echelon form is as follows:
$$\begin{bmatrix}
    1 & 0 & 0 \\
    0 & 1 & 0 \\
    0 & 0 & 0 \\
\end{bmatrix}$$
But I'm not entirely sure which trigonometric identities that it's using to figure that out. Any ideas?

The characteristic polynomial is
$$
\det(\lambda I-A)=
\det\begin{bmatrix}
\lambda-\cos\theta & \sin\theta & 0 \\
-\sin\theta & \lambda-\cos\theta & 0 \\
0 & 0 & \lambda-1
\end{bmatrix}
=
(\lambda-1)
\det\begin{bmatrix}
\lambda-\cos\theta & \sin\theta\\
-\sin\theta & \lambda-\cos\theta
\end{bmatrix}
$$
and, finally,
$$
(\lambda-1)(\lambda^2-2\lambda\cos\theta+1)
$$
and the eigenvalues are $1$, $\cos\theta+i\sin\theta$ and $\cos\theta-i\sin\theta$.
An eigenvector for $1$ is readily computable as $[0\ 0\ 1]^T$, but if you're interested in finding a row echelon form of
$$
\begin{bmatrix}
1-\cos\theta & \sin\theta & 0 \\
-\sin\theta & 1-\cos\theta & 0 \\
0 & 0 & 0
\end{bmatrix}
$$
you can do as follows; set $\theta=2\varphi$, so $1-\cos\theta=2\sin^2\varphi$ and $\sin\theta=2\sin\varphi\cos\varphi$.
The case in which $\sin\varphi=0$ must be dealt with separately; assume $\sin\varphi\ne0$:
\begin{align}
\begin{bmatrix}
2\sin^2\varphi & 2\sin\varphi\cos\varphi & 0 \\
-2\sin\varphi\cos\varphi & 2\sin^2\varphi & 0 \\
0 & 0 & 0
\end{bmatrix}
&\to
\begin{bmatrix}
1 & \cot\varphi & 0 \\
-2\sin\varphi\cos\varphi & 2\sin^2\varphi & 0 \\
0 & 0 & 0
\end{bmatrix}
&&R_1\gets\frac{1}{2\sin^2\varphi}R_1
\\[6px]&\to
\begin{bmatrix}
1 & \cot\varphi & 0 \\
0 & 2 & 0 \\
0 & 0 & 0
\end{bmatrix}
&&R_2\gets R_2+2R_1\sin\varphi\cos\varphi
\\[6px]&\to
\begin{bmatrix}
1 & \cot\varphi & 0 \\
0 & 1 & 0 \\
0 & 0 & 0
\end{bmatrix}
&&R_2\gets \frac{1}{2}R_2
\\[6px]&\to
\begin{bmatrix}
1 & 0 & 0 \\
0 & 1 & 0 \\
0 & 0 & 0
\end{bmatrix}
&&R_1\gets R_1-R_2\cot\varphi
\end{align}
The case $\sin\varphi=0$ corresponds to $\theta=2k\pi$, when the original matrix is the identity.
