DSolve[Derivative[2][y][x] == -(((B + A*x^2)*y[x])/(x*(-a1 + x^2)*(-a2 + x^2)*(-a3 + x^2))) - ((-((-a1 + x^2)*(-a2 + x^2)*(-a3 + x^2)) + x^2*((-a1 + x^2)*(-a2 + x^2) + (-a1 + x^2)*(-a3 + x^2) + (-a2 + x^2)*(-a3 + x^2)))*Derivative[1][y][x])/(x*(-a1 + x^2)*(-a2 + x^2)*(-a3 + x^2)),y[x],x]