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