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