NIntegrate with AdaptiveMonteCarlo gives different results
- To: mathgroup at smc.vnet.net
- Subject: [mg123350] NIntegrate with AdaptiveMonteCarlo gives different results
- From: "Mantovani Sarti Valentina" <smantovani at fe.infn.it>
- Date: Fri, 2 Dec 2011 07:22:50 -0500 (EST)
- Delivered-to: l-mathgroup@mail-archive0.wolfram.com
- Reply-to: smantovani at fe.infn.it
Hello everybody,
I am dealing with a three dimensional integral and I am using
NIntegrate. Searching the web I found out that the best method is the
AdaptiveMonteCarlo, which also speed up the calculation.
this is the integral, no particular singular points:
test[d_] := (0.5 NIntegrate[(compnx[x, y, z, d] +
compny[x, y, z, d] + compnz[x, y, z, d]), {z, -R - d/2,
R + d/2}, {x, -R, R}, {y, -R, R},
Method -> "MultiDimensionalRule"]*hc) // Timing
so I want to calculate it for different values of the parameter d, but
for example if I successively calculate it for d=4 i get two different
results which differs of about 30-50 MeV.
Does anyone know if it is a bug or how could I fix it?
Thanks in advance.
--
Mantovani Sarti Valentina, Ph.D. student
e-mail: smantovani at fe.infn.it
tel. (+39)0532-974.353
INFN - Sezione di Ferrara
c/o Department of Physics (Blocco C)
Via Saragat 1
44100 Ferrara
Italy