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