Dear Vladimir, On 29.09.19 22:35, Vladimir V. Kisil wrote:
GiNaC tutorial in section "2.2 What it can do for you" says that Ginsh shall produce:
> series(tgamma(x),x==0,3); x^(-1)-Euler+(1/12*Pi^2+1/2*Euler^2)*x+ (-1/3*zeta(3)-1/12*Pi^2*Euler-1/6*Euler^3)*x^2+Order(x^3) > evalf(%); x^(-1)-0.5772156649015328606+(0.9890559953279725555)*x -(0.90747907608088628905)*x^2+Order(x^3)
However, with the current GiNaC versions the second output will be identical to the first. I am attaching a small patch which shall restore the expected output.
Thanks for reporting this! It is a regression introduced in f8c2455fbb. I assume you just added the eval() method to it.rest in pseries::eval() because of symmetry? There shouldn't be any need to do so since 1.7.0. All my best, -richy. -- Richard B. Kreckel <https://in.terlu.de/~kreckel/>