# Order of integration changes output of indefinite multiple integral in Mathematica 7

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]) *)

-

## 1 Answer

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.

-