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