Integration by change to polar coordinates

Hey,

I have the following integral:
$$I = \int_0^\infty dx_1\ldots \int_0^\infty dx_{3N} \int_0^\infty dy_1\ldots \int_0^\infty dy_{3N}\Theta\left(1-\sum_{j=1}^{3N}\left(|x_j|^2+|y_j|^2\right)\right)$$
Now I want to change to polar coordinates by the following substitution:

$$x_j &=& r_j \cos \varphi_j; y_j &=& r_j \sin\varphi_j$$
Then the integral becomes

$$\left[\int_0^{2\pi}d\varphi\right]^{3N}\int_0^\infty dr_1\ldots \int_0^\infty dr_{3N}r_1r_2\ldots r_{3N} \Theta\left(1-\sum_{j=1}^{3N}|r_j|^2\right)$$
Now this factors $x_j$ in front of the theta function are bothering me, because if they weren't there, the integral would be simply the volume of sphere with radius 1 in 3N dimensions, and this is a known formula.
But, they are there so this is a problem...
Or not, I don't know, is there another maybe easier way to solve this integral?
Is my substitution perhaps not wel chosen?

(there are some typo's in the formula, but not important, seems like the preview doesn't renew the latex code)

mathman
I don't see any xj only rj, which makes perfect sense.

I mean the $r_j$(instead of $x_j$) coming from the jacobian due to the coordinate transformation.
I wish they weren't there because the integral would be a volume of a 3N-sphere.

mathman
For better or worse, switching from rectangular to polar coordinates will always look like:

dxdy = rdrdφ. Intuitively you can look at the r factor being necessary to preserve the dimension of the differential as length squared.

Yes, of course.

Is there maybe anothor way to factorize the original integral in x and y so that I can solve it an a closed form?

The change to polar coordinates seemed logical to me because I would get a product of two indepent integrals.

mathman

$\Theta$ is Heaviside's stepfunction, equals 1 when it's argument is larger than zero, and 0 otherwise.
That's why it would be the volume of a sphere with radius 1 when being factorized without the r_j in front

mathman
Have you tried u_j=r_j2? Then let v=Σu_j and play with the limits.

The substitution $r_j^2=u_j$
for the r-integral only,
to
$$\frac{1}{2^{3N}}\int_0^\infty du_1\ldots \int_0^\infty du_{3N}\Theta\left(1-\sum_{j=1}^{3N}|u_j|\right)$$

Playing with the limits... I don't see exactly what you mean.

mathman
I was suggesting a change of variables where w=Σuj replaces one of the u's. Now the problem is adjusting the integration limits to accomodate the change of variable.

I still don't get it.
I think you mean this the following?

$$\frac{1}{2^{3N}}\int_0^\infty du_1\ldots \int_0^\infty du_{3N}\Theta\left(1-\sum_{j=1}^{3N}|u_j|\right)$$
with the substitution $w=\sum_j u_j$ you get $dw = \sum du_j$. How exactly do you enter this differential in the integral then...this become nested integrals then?

mathman
No! I meant substitute w for (for example) u1, so dw =du1, while all the other u's remain the same. As a result, the problem then becomes setting up the appropriate integration limits.

One further suggestion (possible alternative) - try setting up the integration limits so that you can dispense with the Heaviside function. Then u1 has limits (0,1), u2 has limits (0,1-u1), u3 has limits (0,1-u1-u2), etc. Then you won't need any change of variables. It looks like a somewhat messy arithmetic problem.

mathman
I finally found the solution. The surface ∑uk=1 is a simplex in 3N dimensions. Your integral is the volume between 0 and the simplex (with all sides limited to uk ≥ 0). This happens to be a solved problem in simplex analysis (V = 1/(3N)!). It is relatively easy to prove by induction.

http://en.wikipedia.org/wiki/Simplex

Thank you for it!
I didn't even know it was a known "standard"-integral.
Even never heard of simplex...

mathman
A simplex is an n dimensional surface with ∑uk=1 and all uk non-negative.

Just curious, but what kind of course is this? I'm finding these advanced calculus topics interesting, but I don't know what they're called. Is it analysis?

Thanks for the link mathman. What I meant by course was where do people learn about this kind of math? i.e. what's it called? I've taken up to Calc III and diff eq so far, but I haven't seen anything like this.

HallsofIvy