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))