DSolve[Derivative[2][y][x] == -(((s + a*r*x^b)*y[x])/(x^2*(-1 + a*x^b))) - ((q + a*p*x^b)*Derivative[1][y][x])/(x*(-1 + a*x^b)),y[x],x]