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