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