dsolve(diff(diff(y(x),x),x) = -b^2/(a^2+x^2)^2*y(x),y(x))