From: Rainer Schöpf <rainer.schoepf@gm...>  20120929 12:11:18

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 