New in Wolfram Mathematica 6: New Generation Numerical Integration  previous | next 
Use Adaptive Monte Carlo in n Dimensions
Mathematica 6 provides state-of-the-art adaptive Monte Carlo integration in any number of dimensions.
In[1]:=

Click for copyable input
ListPlot[Reap[

   NIntegrate[

    Exp[-((x - 2 \[Pi])^2 + (y - \[Pi])^2)], {x, -15, 15}, {y, -15, 

     15}, Method -> "AdaptiveMonteCarlo", 

    EvaluationMonitor :> Sow[{x, y}]]][[2, 1]], 

 PlotRange -> {{-15, 15}, {-15, 15}}, AspectRatio -> 1]
Out[1]=