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