[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