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