[Date Index] [Thread Index] [Author Index]

Re: Random points in NIntegrate

• To: mathgroup at smc.vnet.net
• Subject: [mg56447] Re: Random points in NIntegrate
• From: David Bailey <dave at Remove_Thisdbailey.co.uk>
• Date: Tue, 26 Apr 2005 01:32:46 -0400 (EDT)
• References: <d4hvel$1g7$1@smc.vnet.net>
• Sender: owner-wri-mathgroup at wolfram.com

Ismail T. wrote:
> Hi,
> 
> I need to do the following. I am using mathematica to take a
> multi-dimensional integral with NIntegrate and QuasiMonteCarlo method
> is default for such cases. There are two things that I couldn't figure
> out. First, I want to print the random points that the mathematica
> chooses to take each data points. This is particularly important to
> see whether mathematica is stuck only some parts of the region of
> integration where it assumes that the biggest cntribution comes from.
> 
> Secondly, how can I make mathematica to take sample (random) points
> only in a specific part of the region of integration without applying
> any particular cut to the limits of integration?
> 
> Thanks a lot in advance,
> 
> ismail
> 
The easiest way to print out the random points chosen by NIntegrate is 
to write the integrand as a call to a function - say f[x,y,z] - which 
you define to print out (or otherwise record) its arguments and then 
evaluate and return the integrand at that point.

I am not sure I understand your second question.

David Bailey
dbaileyconsultancy.co.uk