Solving an initial value problem involving the nonlinear ordinary second order differential equation \$ y & # 39; & # 39; (x) = y (x) cdot y # (x) + (y # (x)) ^ 2 \$

I have the following equation $$y & # 39; (x) = y (x) cdot y # (x) + (y # (x)) ^ 2$$ with the initial values $$y (1) = 0$$ and $$y (1) = -1$$.

I am looking for advice on the best way to tackle this particular problem.

This article suggests that a substitution v can reduce the problem in a certain way. So I leave $$v = y (x)$$ as following:

begin {align *} y & # 39; & # 39; (x) & = y (x) cdot y (x) + (y (x)) ^ 2 & (1) \ y & # 39; & # 39; (x) & = y (x) cdot v + v ^ 2 & (2) \ end {align *}

When I enter (2) in wolfram, the engine says that the equation is second-rate linear ordinary differential equation.

Can we use substitution from here to solve the problem? Or is there a better way to do that?