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