dsolve(diff(y(x),x) = -(-1/x-_F1(y(x)^2-2*x))/(y(x)^2)^(1/2)*x,y(x))