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