On Sat, 29 Sep 2012 at 18:30 +0700, Andrey G. Grozin wrote:
> In the current reduce (svn 1759), if I say
>
> load_package specfn; int(exp(x^2),x,0,infinity);
>
> this produces
>
> sqrt(pi)*erf(infinity)
> 
> 2
>
> If I don't load specfn, I get
>
> sqrt(pi)
> 
> 2
>
> So, loading specfn reduces the quality :( Can this be avoided?
Yes. The problem are these two rules in the specfn package:
int(1/e^(~tt^2),~tt,0,~z) => erf(z)/2*sqrt(pi),
int(1/e^(~tt^2),~tt,~z,infinity) => erfc(z)/2*sqrt(pi),
which are used even if z is infinity. The obvious solution is to replace them by
int(1/e^(~tt^2),~tt,0,~z) => erf(z)/2*sqrt(pi) when z freeof infinity,
int(1/e^(~tt^2),~tt,~z,infinity) => erfc(z)/2*sqrt(pi) when z freeof infinity
I'll make the change later today.
Rainer
