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