Numerical Multiple Integration

  • Thread starter Korbid
  • Start date
  • #1
17
0

Main Question or Discussion Point

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!
 

Answers and Replies

  • #2
SteamKing
Staff Emeritus
Science Advisor
Homework Helper
12,796
1,667
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!
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
17
0
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.
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
kreil
Insights Author
Gold Member
668
67
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
17
0
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?
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
 

Related Threads on Numerical Multiple Integration

  • Last Post
Replies
1
Views
674
  • Last Post
Replies
3
Views
2K
Replies
8
Views
2K
  • Last Post
Replies
3
Views
26K
Replies
1
Views
6K
Replies
16
Views
628
  • Last Post
Replies
3
Views
847
Replies
27
Views
1K
  • Last Post
Replies
2
Views
2K
Replies
4
Views
61K
Top