MathGroup Archive 2005

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

Search the Archive

Re: High precision NIntegrate problem. Please help!

  • To: mathgroup at smc.vnet.net
  • Subject: [mg60328] Re: High precision NIntegrate problem. Please help!
  • From: Paul Abbott <paul at physics.uwa.edu.au>
  • Date: Tue, 13 Sep 2005 06:07:07 -0400 (EDT)
  • Organization: The University of Western Australia
  • References: <dfueda$1tu$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

In article <dfueda$1tu$1 at smc.vnet.net>, a_noether_theorem at yahoo.com 
wrote:

> I'm trying to numerically evaluate a particularly nasty integral, and
> no matter what I do, I can't get Mathematica to give me an answer without
> complaining, and the answer seems to change as I change the precision
> options, and continues to change up to arbitrary levels of precision.
> I could really use some help.
> 
> I've tried every combination of options I could, but nothing works.
> 
> In case you care what the function is, it's the following.  Otherwise,
> skip what's between the lines.

Actually, what's between the lines is most important. Sometimes, 
numerical integration is the only way to proceed, but if symbolic 
integration is possible, then it's likely to be more useful.

And the correct approach to using NIntegrate often depends on the form 
of the integrand.

> -------------------------------------------------------------
> I 1st evaluate the following sum:
> 
> 1/( I omega_n - e0) 1/( I omega_n - e1) 1/( I omega_n - e2) 1/( I
> omega_n - e3)
> 
> from n=-Infinity to Infinity, where omega_n = (2 n + 1) Pi/beta
> 
> The sum can be done analytically by contour integration and Mathematica can
> also do the sum.  Anyhoo, it leads to a longish expression that I need
> to integrate.  Here,
> 
> e0 = t1 Cos[x] + t2 Cos[y] - mu
> e1 = t1 Cos[x+X1] + t2 Cos[y+Y1] - mu
> e2 = t1 Cos[x+X2] + t2 Cos[y+Y2] - mu
> e3 = t1 Cos[x+X3] + t2 Cos[y+Y3] - mu
> 
> and t1, t2, mu, beta, X1.., Y1... are just some numbers I choose.  This
> is a Feyman diagram calculation, if anyone cares.   I evaluate this
> with a few {Xi,Yi} sets.
> ----------------------------------------------------------------------

You talk about integration -- but you don't mention what you are 
integrating (x and y I presume), or the integration range.

Also, can you compute the integral of the summand symbolically, and then 
perform the doubly-infinite summation?

Many computations involving Feynman diagrams can be found in closed form 
and solved using tricks such as parametric differentiation, and these 
should be investigated prior to using NIntegrate.

Cheers,
Paul

_______________________________________________________________________
Paul Abbott                                      Phone:  61 8 6488 2734
School of Physics, M013                            Fax: +61 8 6488 1014
The University of Western Australia         (CRICOS Provider No 00126G)    
AUSTRALIA                               http://physics.uwa.edu.au/~paul


  • Prev by Date: Set of strings reducing problem (Again!)
  • Next by Date: Re: Why am I getting this error?
  • Previous by thread: High precision NIntegrate problem. Please help!
  • Next by thread: Re: High precision NIntegrate problem. Please help!