series expansion once again...

Jens Vollinga vollinga at thep.physik.uni-mainz.de
Mon Apr 5 15:25:11 CEST 2004


Hi.

On Mon, Apr 05, 2004 at 05:11:39PM +0400, Sheplyakov Alexei wrote:
> Hello!
> My naive expectations are:
> 	1) first and second lines should be the same
> 	2) third line should be Order(epsilon)
> 
> But the output is:
> 
> tgamma(2-epsilon)*tgamma(3-epsilon)^(-1)*epsilon^2+tgamma(2-epsilon)*tgamma(3-epsilon)^(-1)*epsilon^(-1) = Order(epsilon^(-1))
> tgamma(2-epsilon)*tgamma(3-epsilon)^(-1)*epsilon^(-1) = 1/2*epsilon^(-1)+1/4+Order(epsilon)
> tgamma(2-epsilon)*tgamma(3-epsilon)^(-1)*epsilon^2 = Order(epsilon^(-1))
> 
> 	
> Is this a bug or a feature?

It's a bug, but it has already been fixed in CVS.

Bye,
Jens




More information about the GiNaC-list mailing list