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