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


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