Numerical Multiple Integration

1. Apr 2, 2016

Korbid

I'm trying to solve numerically this multiple integral. But i don't know how to calculate it with Mathamtica or Sage software.

$$\int{e^{-(\vec{v}^2_1+\vec{v}^2_2)}e^{-E(\tau)}}d\vec{r}_1d\vec{r}_2d\vec{v}_1d\vec{v}_2$$
$$E(\tau)=\frac{k}{\tau^2}e^{-\tau/\tau_0}$$
$$\tau(\vec{r}_{12};\vec{v}_{12})=\frac{b-\sqrt{b^2-ac}}{a}$$
$$a=||\vec{v}_{12}||$$
$$b=\vec{r}_{12}\cdot\vec{v}_{12}$$
$$c=||\vec{r}_{12}|| - (2R)^2$$

Thank you!

2. Apr 2, 2016

SteamKing

Staff Emeritus
What you have presented is insufficient. The bounds of the integration are not apparently specified, nor are the variables identified in any way.

Much more context is needed to understand what you want to do.

3. Apr 2, 2016

Korbid

I'm trying to calculate the second virial coefficient for a potential E that depends on positions and velocities. The N particles are moving inside LxL square.

$$\int^L_0\int^L_0\int^{\infty}_0\int^{\infty}_0{e^{-(\vec{v}^2_1+\vec{v}^2_2)}e^{-E(\tau)}}d\vec{r}_1d\vec{r}_2d\vec{v}_1d\vec{v}_2$$
$$E(\tau)=\frac{k}{\tau^2}e^{-\tau/\tau_0}$$
$$\tau(\vec{r}_{12};\vec{v}_{12})=\frac{b-\sqrt{b^2-ac}}{a}$$
$$a=||\vec{v}_{12}||$$
$$b=\vec{r}_{12}\cdot\vec{v}_{12}$$
$$c=||\vec{r}_{12}|| - (2R)^2$$

R, tau_0 and k are constants.

4. Apr 2, 2016

kreil

I don't have much experience using Mathematica or Sage. But most numerical integration algorithms use iteration to solve such problems anyway, so as long as you know how to calculate a single integral with one of those programs you can use the same technique to calculate 4 nested ones. Just evaluate the integrals from the inside out, and feed the result of the first to the second, from the second to the third, and so on. If the constants give you trouble, just set them all to 1 to start with. (There might be a function supplied that is designed for doing an arbitrary number of integrations as well, so look for that.)

Also, you said you want to numerically integrate these equations, but the information you've provided is only really sufficient for a symbolic integration. For example, how would the software evaluate b without numeric expressions for r12 and v12?

5. Apr 2, 2016

Korbid

I'm sorry, i forgot it.
r12 is the relative position and v12 is the relative velocity
$$r_{12}=r_1-r_2$$
it's the same for v12