Reescribamos la EDO como un sistema de primer orden:
$$\begin{align*} y_1 &= y\\ y_2 &= y' \\ y'_1 &= y_2 \\ y'_2 &= 2y_1. \end{align*}$$
Este es un sistema lineal, así que dejemos que $\mathbf{y} = \begin{pmatrix} y_1 \\ y_2 \end{pmatrix}$ . Entonces,
$$\mathbf{y}' = \underbrace{\begin{pmatrix} 0 & 1 \\ -2 & 0 \end{pmatrix}}_{=A} \mathbf{y}.$$
Ahora, podemos utilizar la exponencial de la matriz para resolver esta EDO de la manera habitual, pero más bien vamos a explorar la ecuación característica de la matriz $A$ :
$$A-\lambda I = \begin{pmatrix} -\lambda & 1 \\ -2 & -\lambda\end{pmatrix}\\ \det A-\lambda I = (-\lambda)^2 + 2.$$
Si se iguala a cero, se obtiene
$$\lambda^2 +2 = 0$$
que es precisamente la ecuación característica obtenida mediante el método "tradicional" de conversión de derivadas de $y$ a los poderes de $\lambda$ .