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