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