\[ y'(x)=\frac {1}{-\left (x^2 \left (\frac {1}{y(x)}+1\right ) \text {$\_$F1}\left (x \left (\frac {1}{y(x)}+1\right )\right )\right )+x^2 \text {$\_$F1}\left (x \left (\frac {1}{y(x)}+1\right )\right )+x \left (\frac {1}{y(x)}+1\right )-x} \] ✓ Mathematica : cpu = 0.506831 (sec), leaf count = 365
DSolve[Derivative[1][y][x] == (-x + x*(1 + y[x]^(-1)) + x^2*_F1[x*(1 + y[x]^(-1))] - x^2*(1 + y[x]^(-1))*_F1[x*(1 + y[x]^(-1))])^(-1),y[x],x]
\[\text {Solve}\left [\int _1^{y(x)}\left (\frac {x \text {$\_$F1}\left (x \left (1+\frac {1}{K[2]}\right )\right )-1}{x \text {$\_$F1}\left (x \left (1+\frac {1}{K[2]}\right )\right ) K[2]-K[2]+x \text {$\_$F1}\left (x \left (1+\frac {1}{K[2]}\right )\right )}-\int _1^x\left (\frac {\text {$\_$F1}\left (K[1] \left (1+\frac {1}{K[2]}\right )\right )-\frac {K[1] \text {$\_$F1}'\left (K[1] \left (1+\frac {1}{K[2]}\right )\right )}{K[2]}-\frac {K[1] \text {$\_$F1}'\left (K[1] \left (1+\frac {1}{K[2]}\right )\right )}{K[2]^2}}{K[1] \left (K[2] \text {$\_$F1}\left (K[1] \left (1+\frac {1}{K[2]}\right )\right )+\text {$\_$F1}\left (K[1] \left (1+\frac {1}{K[2]}\right )\right )\right )-K[2]}-\frac {\left (K[2] \text {$\_$F1}\left (K[1] \left (1+\frac {1}{K[2]}\right )\right )+\text {$\_$F1}\left (K[1] \left (1+\frac {1}{K[2]}\right )\right )\right ) \left (K[1] \left (\text {$\_$F1}\left (K[1] \left (1+\frac {1}{K[2]}\right )\right )-\frac {K[1] \text {$\_$F1}'\left (K[1] \left (1+\frac {1}{K[2]}\right )\right )}{K[2]}-\frac {K[1] \text {$\_$F1}'\left (K[1] \left (1+\frac {1}{K[2]}\right )\right )}{K[2]^2}\right )-1\right )}{\left (K[1] \left (K[2] \text {$\_$F1}\left (K[1] \left (1+\frac {1}{K[2]}\right )\right )+\text {$\_$F1}\left (K[1] \left (1+\frac {1}{K[2]}\right )\right )\right )-K[2]\right ){}^2}\right )dK[1]\right )dK[2]+\int _1^x\left (\frac {y(x) \text {$\_$F1}\left (K[1] \left (1+\frac {1}{y(x)}\right )\right )+\text {$\_$F1}\left (K[1] \left (1+\frac {1}{y(x)}\right )\right )}{K[1] \left (y(x) \text {$\_$F1}\left (K[1] \left (1+\frac {1}{y(x)}\right )\right )+\text {$\_$F1}\left (K[1] \left (1+\frac {1}{y(x)}\right )\right )\right )-y(x)}-\frac {1}{K[1]}\right )dK[1]=c_1,y(x)\right ]\] ✓ Maple : cpu = 0.213 (sec), leaf count = 40
dsolve(diff(y(x),x) = 1/(-x+(1/y(x)+1)*x+_F1((1/y(x)+1)*x)*x^2-_F1((1/y(x)+1)*x)*x^2*(1/y(x)+1)),y(x))
\[y \left (x \right ) = {\mathrm e}^{\operatorname {RootOf}\left (-\textit {\_Z} -\left (\int _{}^{\frac {x \,{\mathrm e}^{\textit {\_Z}}}{{\mathrm e}^{\textit {\_Z}}-1}}\frac {1}{\textit {\_a} \left (\textit {\_F1} \left (\textit {\_a} \right ) \textit {\_a} -1\right )}d \textit {\_a} \right )+c_{1}\right )}-1\]