The conic is $\Gamma(x,y):8x^2+8xy+2y^2-5y=0$.
Since
\begin{equation}
det(A)=\begin{pmatrix} 8&4&0\\4&2&-\frac{5}{2}\\0&-\frac{5}{2}&0\end{pmatrix}=-50\neq 0
\end{equation}
\begin{equation}
det(Q)=\begin{pmatrix} 8&4\\4&2\end{pmatrix}=0
\end{equation}
$\Gamma$ is an parabola.
By searching the eigenvalues of $Q$ I get:
\begin{equation}
det(Q-\lambda I)=\begin{vmatrix} 8-\lambda&4\\4&2-\lambda\end{vmatrix}=\lambda(\lambda-10)=0
\end{equation}
i.e. $\lambda_{1}=0$ and$\lambda_{2}=10$.
The eigenspaces are $E(\lambda_{1})=\mathcal{L}(-\frac{1}{2},1)$ while $E(\lambda_{2})=\mathcal{L}(2,1)$. The basis of such spaces are evidently orthogonal basis of their respective eigenspace. The first one, $v=(-\frac{1}{2},1)$ can be normalized by dividing it by $\frac{\sqrt{5}}{2}$ while the second, $w=(2,1)$, by $\sqrt{5}$. Hence, the rotation matrix $R$ is:
\begin{equation}
\tilde{R}=\begin{pmatrix} -\frac{1}{2}&2\\1&1\end{pmatrix} \begin{pmatrix}\frac{2}{\sqrt{5}}\\\frac{1}{\sqrt{5}} \end{pmatrix} =\begin{pmatrix} -\frac{1}{\sqrt{5}}&\frac{2}{\sqrt{5}}\\\frac{2}{\sqrt{5}}&\frac{1}{\sqrt{5}}\end{pmatrix}
\end{equation}
This matrix has determinant $-1$ while it has to be $det(R)=1$, thus I exchange sign:
\begin{equation}
R=\begin{pmatrix} \frac{1}{\sqrt{5}}&\frac{2}{\sqrt{5}}\\-\frac{2}{\sqrt{5}}&\frac{1}{\sqrt{5}}\end{pmatrix}
\end{equation}
It follows that:
\begin{equation}
\begin{cases} x=\frac{1}{\sqrt{5}}(\tilde{x}+2\tilde{y})\\y=\frac{1}{\sqrt{5}}(-2\tilde{x}+\tilde{y})\end{cases}
\end{equation}
and, by substituting in the equation of the conic and simplifying, I get:
\begin{equation}
\Gamma(\tilde{x},\tilde{y}):10\tilde{y}^2+2\sqrt{5}\tilde{x}-\sqrt{5}\tilde{y}=0
\end{equation}
Then I should translate to obtain the canonic form (by completing the squares). Even if I'm pretty convinced of the procedure and that I checked the calculations with the PC, I do not understand why I have a term in $x$ instead that in $x^2$. Have I made some mistakes?
I'd really appreciate a check of the procedure and/or any advise, thank you!