[GiNaC-devel] pseries::evalf() and pseries::eval()

Vladimir V. Kisil kisilv at maths.leeds.ac.uk
Sun Sep 29 22:35:14 CEST 2019


	Dear All,

	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.

  Best wishes,
  Vladimir
-- 
Vladimir V. Kisil                 http://www.maths.leeds.ac.uk/~kisilv/
  Book:     Geometry of Mobius Transformations     http://goo.gl/EaG2Vu
  Software: Geometry of cycles          http://moebinv.sourceforge.net/
  Jupyter (Colab):         https://github.com/vvkisil/MoebInv-notebooks
  Jupyter (CodeOcean):       https://codeocean.com/capsule/7952650/tree

-------------- next part --------------
A non-text attachment was scrubbed...
Name: 0001-Fix-pseries-eval-and-pseries-evalf-methods.patch
Type: text/x-diff
Size: 1203 bytes
Desc: GiNaC patch
URL: <http://www.ginac.de/pipermail/ginac-devel/attachments/20190929/0888c6f1/attachment.bin>


More information about the GiNaC-devel mailing list