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