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