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