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