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