[GiNaC-list] series of x^n around x=0

Vladimir V. Kisil kisilv at maths.leeds.ac.uk
Thu Dec 19 21:32:34 CET 2019


>>>>> On Thu, 19 Dec 2019 09:59:19 +0000, Feng Feng <f.feng at outlook.com> said:
  > Hi all,
  > I got a question while using series on x^n with n a non-integer.
  > here is an example:
  > ------
  > symbol x("x");
  > ex expr = pow(x, ex(3)/2);  // 3/2 may be a general non-integer
  > cout << expr.series(x,2) << endl;
  > ------
  > The output is: Order(x^2)
  > so is there a way to handle this?

  Note, that series expansion of x^(3/2) at x==0 is not well-defined in
  mathematical sense. GiNaC signals this by throwing pole error at
  expr.series(x,3).

  Thus it is not clear what do you mean by "is there a way to handle
  this?"
-- 
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


More information about the GiNaC-list mailing list