\[ y'(x)=\frac {F(-(x-y(x)) (y(x)+x))+x}{y(x)} \] ✓ Mathematica : cpu = 30.5687 (sec), leaf count = 110
\[\text {Solve}\left [\int _1^{y(x)} \left (-\int _1^x -\frac {2 K[1] K[2] F'(-(K[1]-K[2]) (K[1]+K[2]))}{F(-(K[1]-K[2]) (K[1]+K[2]))^2} \, dK[1]-\frac {K[2]}{F(-(x-K[2]) (K[2]+x))}\right ) \, dK[2]+\int _1^x \left (\frac {K[1]}{F(-(K[1]-y(x)) (K[1]+y(x)))}+1\right ) \, dK[1]=c_1,y(x)\right ]\]
✓ Maple : cpu = 0.129 (sec), leaf count = 53
\[ \left \{ y \left ( x \right ) =\sqrt {{x}^{2}+{\it RootOf} \left ( -2\,x+\int ^{{\it \_Z}}\! \left ( F \left ( {\it \_a} \right ) \right ) ^{-1}{d{\it \_a}}+2\,{\it \_C1} \right ) },y \left ( x \right ) =-\sqrt {{x}^{2}+{\it RootOf} \left ( -2\,x+\int ^{{\it \_Z}}\! \left ( F \left ( {\it \_a} \right ) \right ) ^{-1}{d{\it \_a}}+2\,{\it \_C1} \right ) } \right \} \]