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