[GiNaC-list] series((x+x^2)^2,x,0) is broken

Richard B. Kreckel kreckel at in.terlu.de
Sat Jun 24 17:47:32 CEST 2023


On 6/23/23 12:28, Vitaly Magerya wrote:
>> For one thing, I suppose that the Order(x) function is missing a power
>> evaluation of the kind Order(x)^e -> Order(x^e). Does this sound right?
> 
> Thanks for looking into this.
> 
> The substitution above is not correct if e<0. For example:
> 
>      x = Order(1/A) => |x| <= Const * 1/A
> 
> but
> 
>      x = 1/Order(A) => |x| >= 1 / (Const * A)
> 
> so
> 
>      Order(1/A) =/= 1/Order(A)
> 
> It is correct for non-negative integer e though.

Wait.

Consider a series which terminates in Order(x^n). We don't make any 
statement about the magnitude or the sign of the constant in front of 
the Order(x^n) term. Likewise for Order(x)^n. Hence, I think the 
transformation Order(x)^n -> Order(x^n) is correct, even if n<0.

Comments or objections?

    -richard.


More information about the GiNaC-list mailing list