Ilya Zakharevich on Sun, 07 Feb 2010 03:32:17 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: PARI stable release 2.3.5 (for real now) 

On Fri, Feb 05, 2010 at 11:50:05AM +0100, Bill Allombert wrote:
> On Wed, Feb 03, 2010 at 12:14:57AM -0800, Ilya Zakharevich wrote:
> >   *) intnum() is still borken (same as in 3.0):
> > 
> >     ? intnum(Z=[0,-1/2],[[1],-3/2],1/sqrt(Z)*(1+Z)^(-3/2))
> >     %1 = 7.150393303262919781690330041 + 1018.215229301285883563073358 I
> > 
> >      (I presume this is integral from 0 to infinity...)

> You need to split the integral:

I know that this works around the bug.

> ? intnum(Z=[0,-1/2],1,1/sqrt(Z)*(1+Z)^(-3/2))+intnum(Z=1,[[1],-3/2],1/sqrt(Z)*(1+Z)^(-3/2))
> %36 = 2.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000
> 
> This belongs to a class of problem inherent to the method. 

Then it should split automatically.

Yours,
Ilya