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