Mathematica Stack Exchange is a question and answer site for users of Mathematica. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I'm trying to integrate a form-factor used in the calculation of radiation between two rectangles in perpendicular planes. While the integral is usually done over fixed limits, I am trying to do the indefinite multiple integral. When I change the order of integration the results differ in one term.

Because the result gets used in a summation that is essentially evaluating the integral at the end points representing the rectangle limits, the extra terms do not matter i.e. they end up cancelling out. But I'm bothered that they are being generated to begin with and that they differ when I change the order of integration.

In order to make sure that the difference is not because of the possibility of imaginary numbers or zeros, cropping up, I set assumptions that limit the variables to all be finite positive numbers. I'm guessing the differences are because of the log's or arctan's that crop up during the integration, but I can't figure out why.

In the example below the two integrals differ in the term x^2 , 2yw-y^2

The effects are not from FullSimplify as that just collects terms - it doesn't combine log or arctan terms.

Can anyone explain the results I get? I'm running Mathematica 7.0 Thanks in advance! (sorry I can't post an image but I'm a new user so it's not allowed)

Fperp = FullSimplify[
  1/Pi Integrate[(x*z)/(x^2 + z^2 + (y - w)^2)^2, x, y, z, w, 
    Assumptions -> {x > 0 && w > y0 & z > 0 && y > 0 && 
       x <= Infinity && y < Infinity && z < Infinity && 
       w < Infinity}]]

(* ===>  (x^2 + 4 (w - y) Sqrt[x^2 + z^2]
   ArcTan[(w - y)/Sqrt[x^2 + z^2]] + (-x^2 + (w - y)^2 - z^2) Log[
   x^2 + (w - y)^2 + z^2])/(8 \[Pi]) *)

Fperp2 = FullSimplify[(1/Pi)*
     Integrate[(x*z)/(x^2 + z^2 + (y - w)^2)^2, w, z, y, x, 
     Assumptions -> {x > 0 && w > 0 && z > 0 && y > 0 && 
       x < Infinity && y < Infinity && z < Infinity && w < Infinity}]]

(* ===>  (-w (w - 2 y) + 4 (w - y) Sqrt[x^2 + z^2]
   ArcTan[(w - y)/Sqrt[x^2 + z^2]] + (-x^2 + (w - y)^2 - z^2) Log[
   x^2 + (w - y)^2 + z^2])/(8 \[Pi]) *)
share|improve this question

There is no reason why the two integrals in your problem should be the same, and they aren't. They aren't definite integrals, only anti-derivatives (indefinite integrals). That is, the only property that they must have is the following:

D[Fperp, x, y, z, w] == D[Fperp2, x, y, z, w]

(* ==> True *)

Indeed, they are antiderivatives of the same function. The reason you got antiderivatives that differ in form is that the integration algorithm generally introduces different substitutions and defines different integration constants depending on the form in which you present the integrand. That's really all that needs to be said in general - now you could still ask for the reason for the specific differences, but that would best be explained by doing the multiple integrals one at a time, and I don't think it would be easy to draw any more general statements from that.

share|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.