There is a solution that does not involve a substitution... if that's any help...
First, multiply through by [tex]x + y^2[/tex], to get
[tex]x y^{\prime} + y^2 y^{\prime} = y[/tex]
rearrange to get
[tex]x y^{\prime} - y = -y^2 y^{\prime}[/tex]
but
[tex]x y^{\prime} - y = y^2 ( \phi - \frac{x}{y})^{\prime}[/tex]
(where [tex]\phi[/tex] is a constant.) So,
[tex]( \phi - \frac{x}{y})^{\prime} = -y^{\prime}[/tex]
which you can integrate to get
[tex]\phi - \frac{x}{y} = - y[/tex]
which you can turn into a quadratic by multiplying through by [tex]y[/tex], leaving you with.
[tex]y(x) = \frac{-\phi \pm \sqrt{\phi^2 + 4x}}{2}[/tex]