DSolve[y[x]*(a^2 - b^2*f[x]^2 + (a*Derivative[1][f][x])/f[x]) - (2*a + Derivative[1][f][x]/f[x])*Derivative[1][y][x] + Derivative[2][y][x] == 0,y[x],x]