The goal is a recurrence relation for the series coefficients; maybe the
method, like the Frobenius method, to solve $yy’’=y’$. First, substitute $y=\sum\limits_{k=0}^\infty a_k x^{k+r}$:
$$\sum_{n=0}^\infty\sum_{m=0}^\infty(n+r+1)(n+r+2)a_ma_{n+2}x^{n+r}x^{m+r}=\sum_{k=0}^\infty(k+r+1)a_{k+1}x^{k+r}$$
Equating coefficients with Kronecker $\delta_{u,v}$:
$$(k+r+1)a_{k+1}= \sum_{n=0}^\infty\sum_{m=0}^\infty(n+r+1)(n+r+2)a_ma_{n+2}\delta_{k+r,m+n+2r}$$ and make $n$’s upper bound finite since the double sum is $0$ when $k+r\ne m+n+2r\iff k\ne m+n+r$. Following Quantile Mechanics, erase $\sum\limits_{m=0}^\infty$, $m=k-n-r$, and the upper bound of $n$ being $k-r$ is just a guess:
$$\boxed{(k+r+1)a_{k+1}=\sum_{n=0}^{k-r}(n+r+1)(n+r+2)a_{k-n-r}a_{n+2}}$$
Is this recurrence relation true and if not, what is the recurrence relation for the series coefficients of $yy’’=y’$?