MathGroup Archive 2011

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

Search the Archive

Re: Numerical/analytical integration around a pole

  • To: mathgroup at smc.vnet.net
  • Subject: [mg118911] Re: Numerical/analytical integration around a pole
  • From: Peter Pein <petsie at dordos.net>
  • Date: Sun, 15 May 2011 07:05:43 -0400 (EDT)
  • References: <201105110828.EAA01400@smc.vnet.net> <iqg5t7$d0i$1@smc.vnet.net>

Am 12.05.2011 10:31, schrieb A. Lapraitis:
> Thanks to Gary Wardall and Murray Eisenberg for quick responses.
> 
> I forgot to say I was having the problem on 7.0.1.0 on Linux x86
> (64-bit). Since people are not experiencing the strange behaviour on
> v8, it is fair to conclude that the bug got fixed.
> 

And while using Version 7, try to Rationalize[..,0] the integrand in
Integrate:

In[2]:= Through[{Composition[N,Integrate],NIntegrate}
                  [1/(Rationalize[y[x],0]+I/1000),{x,1/2,7/10}]]
Out[2]= {0.323456-3.37701 I,0.323456-3.37701 I}

but this lasts a little while.

Peter


  • Prev by Date: Combinations and Counting
  • Next by Date: Re: Protein Sequence Alignment efficiency
  • Previous by thread: Re: Combinations and Counting
  • Next by thread: Find Previous doesn't work???