Computational Science people:
I originally posted this question at Math Stack Exchange and someone commented that I might get "much better" answers here:
I am a novice at numerical methods and Matlab. I am attempting to evaluate the following sum of two triple integrals (it can obviously be written more simply, but you still cannot evaluate it symbolically (?)). I am having trouble getting the $\LaTeX$ to work here, so I reluctantly broke it up into pieces here: I want to find the sum of
$$\frac{2}{((1/0.3) - 1)^2}\left(\int_1^{1/0.3}\int_1^{r_1}\int_0^{r_1-r_0}F_1(r_0,r_1,t)\exp(-\frac{(0.3)^2 t^2}{4})\,dt\,dr_0\,dr_1 \right),$$
and
$$\frac{2}{((1/0.3) - 1)^2}\left(\int_1^{1/0.3}\int_1^{r_1}\int_{r_1-r_0}^{r_1+r_0} F_2(r_0,r_1,t)\exp(-\frac{(0.3)^2 t^2}{4})\,dt\,dr_0\,dr_1 \right),$$
where
$$F_1(r_0,r_1,t)=\frac{t^2 r_0^3*(0.3)^3}{2r_1^3\sqrt{\pi}}$$
and
$$F_2(r_0,r_1,t)=\frac{(0.3)^3\pi^{3/2}(r_0+r_1-t)^4 (t^2+2t(r_0+r_1)-3(r_1-r_0)^2)^2}{288(\frac{4}{3}\pi r_0^3)(\frac{4}{3}\pi r_1^3)}.$$
EDIT (March 2 2013): Someone responded that they got Mathematica to do the integrals symbolically. I just attempted to do this (with simplified versions of the integrals) and Mathematica could only do the outer two of the first one, and stalled on the second one. I would appreciate some help. Here is what I did.:
I attempted to evaluate
$$\int_1^2 \int_1^{r_2} \int_0^{r_2-r_1} \frac{r_1^3 t^2 \exp(-t^2)}{r_2^3}\,dt\,dr_1\,dr_2$$ via
Integrate[r1^3/r2^3*t^2*Exp (-t^2), {t, 0, r2 - r1}, {r1, 1, r2}, {r2, 1, 2}]
and Mathematica returns (I had trouble with the $\LaTeX$ here because the result is long. I broke it into two equations. if anyone knows a good way to display this please tell me):
$$\int_1^2 \frac{1}{64r2^2} e^{-1-r2^2}(2e^{2r2}(25+r2(19+2r2(1+r2)))-$$
$$e^{1+r2^2}(32r2(2+r2^2)) +\sqrt{\pi}(11+4r2^2(9+r2^2))\operatorname{Erf}[1-r2])\,dr2.$$
Then I tried to evaluate
$$\int_1^2\int_1^{r_2}\int_{r_2-r_1}^{r_2+r_1} \ldots \qquad \qquad \qquad $$
$$\ldots\frac{\exp(-t^2)(r_1+r_2-t)^4(t^2+2t(r_1+r_2)-3(r_2-r_1)^2)^2}{r_1^3 r_2^3}\,dt\,dr_\,dr_2$$
using
Integrate[(r1 + r2 - t)^4*(t^2 + 2*t*(r1 + r2) - 3*(r2 - r1)^2)^2* Exp[-t^2]/r1^3/r2^3, {r2, 1, 2}, {r1, 1, r2}, {t, r2-r1, r2 + r1}]
just now, and Mathematica has not returned an answer after about half an hour (but I am having computer network problems right now, and they may be to blame).
[END OF MARCH 2 EDIT]
I used Matlab's "triplequad" command, with no extra options. I handled the variable limits of integration by means of heaviside functions, because I didn't know any other way to do it. Matlab gave me $0.007164820144202$.
I know Matlab is good software, but I have heard that numerical triple integrals are hard to do accurately, and mathematicians are supposed to be skeptical, so I want some way to verify the accuracy of this answer. The integrals give the expected value of a certain experiment (if anyone wants, I can edit this question to describe the experiment): I implemented the experiment in Matlab using appropriately randomly generated numbers, a million times, and averaged the results. I repeated this process four times. Here are the results (I apologize if I have used the word "trial" improperly):
Trial 1: $0.007133292603256$
Trial 2: $0.007120455071989$
Trial 3: $0.007062595022049$
Trial 4: $0.007154940168452$
Trial 5: $0.007215000289130$
Although each trial used a million samples, the simulation values only agree in the first significant digit. They are not close enough to each to each other for me to determine whether the numerical triple integral is accurate.
So can anyone tell me whether I can trust the result of "triplequad" here, and under what circumstances one can trust it in general?
One suggestion I got at Math Stack Exchange was to try other software like Mathematica, Octave, Maple, and SciPy. Is this good advice? Do people actually do numerical work in Mathematica and Maple? Octave is kind of a Matlab clone, so can I assume it uses the same integration algorithms? I haven't even heard of SciPy before and would appreciate any opinions about it.
UPDATE: Someone from Math Stack Exchange did it in Maple and got $0.007163085468$. That is agreement to three significant figures. That is a good sign.
Also, I would appreciate suggestions on how to enter long, multi-line expression in $\LaTeX$ in Stack Exchange. Can you use the "aligned" environment here? I tried, and I couldn't get it to work.