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