Convert to Formal Power Series - Maple Help
For the best experience, we recommend viewing online help using Google Chrome or Microsoft Edge.

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Power Series : Convert to Formal Power Series

convert/FormalPowerSeries

convert to formal Laurent-Puiseux series

 

Calling Sequence

Parameters

Returns

Options

Description

Examples

References

Compatibility

Calling Sequence

convert(f, FormalPowerSeries, eq, n, opts)

convert(f, FormalPowerSeries, eq, a(n), opts)

convert(f, FPS, eq, n, opts)

convert(f, FPS, eq, a(n), opts)

Parameters

f

-

algebraic expression

eq

-

equation (e.g. z=z0) or name (e.g. z); optional if f contains only one variable

n

-

(optional) name of the summation variable in the output

a(n)

-

(optional) name for the nth series coefficient in the output, which will be a recursive formula for these coefficients instead of a formal power series

opts

-

sequence of optional equations of the form keyword=value; possible keywords are method, differentialorder, recurrence, or output

Returns

• 

When a hypergeometric type representation is returned, it can either be a single inert Sum, or a linear combination of such Sums, possibly plus some initial terms that are linear combinations of rational powers of zz0 or 1z.

• 

When such a representation cannot be obtained, or when a recursive formula is requested, the output is a sequence of two expressions. The first return value is similar to the representation described above, except that the coefficients of one of the Sums are indeterminate, such as, e.g., n=0anzn. The second return value is an REsol data structure giving a recurrence relation for an, including initial conditions.

Options

• 

differentialorder: positive integer N or list of two positive integers N1,N2; upper bound(s) for the order of the differential equation(s) (DE) internally computed for f. Default: N1=10 for holonomic DEs and N2=4 for quadratic DEs. If differentialorder=N is given, then the upper bound for holonomic DEs is N1=max10,N, and the upper bound N2 for quadratic DEs is 4 if N4, 9 if N9, and N otherwise.

• 

method: one of hypergeometric, holonomic, quadratic, or default; specifies the underlying computational method that will be used. method=hypergeometric looks for a closed-form formula of the representation whereas holonomic and quadratic look for a recursive formula from a holonomic or a quadratic differential equation satisfied by f. The default method applies the three methods one after the other and stops once a result is found.

• 

recurrence: either true or false (default). If recurrence=true then a recursive formula with the coefficient a(n) is sought. This reduces the methods applied to only the two options holonomic and/or quadratic. If a(n) is specified as the 4th argument, then the default is recurrence=true. If both options method and recurrence are given, the former takes precedence.

• 

output: expanded or combined (default); controls the presentation of hypergeometric type power series (method=hypergeometric). output=combined gathers series of same power in one summation, whereas output=expanded puts linearly (rationally) independent coefficients in different summations.

Description

• 

This command computes the power series representation of linear combinations of meromorphic functions and some special functions.

• 

convert(f, FormalPowerSeries, z=z0), or, equivalently, convert(f, FPS, z=z0), tries to find a formal power series expansion for f with respect to the variable z at the point of expansion z0.

• 

If eq is just a name z, then the expansion point z=0 is assumed.

• 

The output series is a hypergeometric type power series, or a recursive formula of the representation sought when a hypergeometric type representation cannot be computed. The user may also specify the desired representation using the method option.

• 

See the references for a precise definition of what is called hypergeometric type power series. It generalizes the family of Laurent-Puiseux series having hypergeometric term coefficients.

• 

The command first looks for a holonomic differential equation satisfied by f (see DEtools[FindODE]); hence Maple must know the derivatives of f.

• 

If the summation variable n is not specified, then one of of the variable names n, n0, n1, etc. is chosen.          

• 

For a complete list of known functions, see inifcns.

Examples

convertcosz,FormalPowerSeries,z=π2

n=0−1nzπ22n+12n+1!

(1)

convertexpx+hypergeoma,b,c,x2,FormalPowerSeries,x,k

k=0xkk!+k=0pochhammerb,kpochhammera,kx2kk!pochhammerc,k

(2)

convertln1+sqrtz+z+z32,FormalPowerSeries

n=0−1nzn2+12n+1+n=0−1nzn+1n+1

(3)

convertzcosz32+arcsinz132,FormalPowerSeries

n=024nn!2z2n3+232n+2!+n=0−1nz3n+12n!

(4)

convert1z2pz3q,FormalPowerSeries,z

n=0p2n2+2q13n3p52+−1np2n22qpn2+12−1np1n2q2+2pq23n3p1n2q2zn2p3q2p32+q+n=0pq1nz3np2+q23p+q43+n=0q23nz3n+1p2+q23p+q43

(5)

Illustration for the output option

convertsinz+cosz5,FormalPowerSeries

n=025n+59n10−1nz2n42n!+n=0525n39n2−1nz2n+142n+1!

(6)

convertsinz+cosz5,FormalPowerSeries,output=expanded

n=0−1n25nz2n42n!+n=05−1n9nz2n42n!+n=05−1nz2n22n!+n=05−1n25nz2n+142n+1!+n=015−1n9nz2n+142n+1!+n=05−1nz2n+122n+1!

(7)

Illustration for the method option. All the outputs below can be obtained by the default option. However, one can directly specify the method that yields the result in order to save time. ⅇz3+z and arcsinz3  are holonomic but not of hypergeometric type.

convertexpz+z3,FormalPowerSeries,method=holonomic

n=0Anzn,RESoln+3An+33AnAn+2=0,An,A0=1,A1=1,A2=12,INFO

(8)

convertarcsinz3,FormalPowerSeries,method=holonomic

n=0Anzn+1,RESoln4+4n3+6n2+4n+1An+2n418n362n298n60An+2+n4+14n3+71n2+154n+120An+4=0,An,A0=0,A1=0,A2=1,A3=0,INFO

(9)

zⅇz1 and 1+tanz2 are not holonomic. Bernoulli numbers:

convertzexpz1,FormalPowerSeries,Bn,method=quadratic

n=0Bnzn,RESolBn+3+Bn+2+_k=1n+2B_kBn+3_kn+4,Bn,B0=1,B1=12,B2=112,INFO

(10)

convert1+tanz2,FormalPowerSeries,method=quadratic

n=0Anzn,RESolAn+2+8An+_k=1n16A_kAn_kn+1n+2,An,A0=1,A1=0,INFO

(11)

Different results from the three methods

convertarcsinz2,FormalPowerSeries,method=hypergeometric

n=024nn!2z2n+22n+2!

(12)

convertarcsinz2,FormalPowerSeries,method=holonomic

n=0Anzn+1,RESoln22n1An+n2+5n+6An+2=0,An,A0=0,A1=1,INFO

(13)

convertarcsinz2,FormalPowerSeries,method=quadratic

n=0Anzn,RESolAn+44n+2An+2+_k=2n_k+1A_k+1n+3_kAn+3_k+_k=2n+2_k+1A_k+1n+5_kAn+5_k4n+3,An,A0=0,A1=0,A2=1,A3=0,INFO

(14)

convertexparcsinz,FormalPowerSeries,method=hypergeometric

n=0pochhammerI2,npochhammerI2,n4nz2n2n!+n=04npochhammer12I2,npochhammer12+I2,nz2n+12n+1!

(15)

convertexparcsinz,FormalPowerSeries,method=holonomic

n=0Anzn,RESoln21An+n2+3n+2An+2=0,An,A0=1,A1=1,INFO

(16)

convertexparcsinz,FormalPowerSeries,method=quadratic

n=0Anzn,RESolAn+3+2n+1An+1+2An+2+_k=1n1_k+1A_k+1n+1_kAn+1_k+_k=1n+1_k+1A_k+1n+3_kAn+3_k+_k=1n+1A_kAn+2_k2n6,An,A0=1,A1=1,A2=12,INFO

(17)

The quadratic method detects non-trivial identities

f2arctanhsin2z1+cos2z:gln1+tanz1tanz:

convertfg,FormalPowerSeries,method=quadratic

0

(18)

Illustration for the recurrence option. One of holonomic and quadratic is used for the method.

convertsinz+cosz,FormalPowerSeries,an,recurrence

n=0anzn,RESolan+n2+3n+2an+2=0,an,a0=1,a1=1,INFO

(19)

convertln1+z1,FormalPowerSeries,an,recurrence

n=0anzn1,RESolan+3+n+1an+2+_k=1n+2a_kan+3_kn+4,an,a0=1,a1=12,a2=112,INFO

(20)

convertln1+sinz,FormalPowerSeries,an,recurrence

n=0anzn,RESolan+4+n+2an+2+n+2n+3an+3+_k=1n_k+1_k+2a_k+2n+2_kan+2_kn+2n+3n+4,an,a0=0,a1=1,a2=12,a3=16,INFO

(21)

convertarctanzexpz,FormalPowerSeries,an,recurrence

n=0anzn,RESolan+2n4an+1+n2+5n+7an+2+2n6an+3+n2+7n+12an+4=0,an,a0=0,a1=1,a2=1,a3=16,INFO

(22)

More examples.

convertexpzlnz,FormalPowerSeries

n=0lnzznn!

(23)

convertexpzsinz,FormalPowerSeries

n=04n−1nz4n+14n+1!+n=024n−1nz4n+24n+2!+n=024n−1nz4n+34n+3!

(24)

convert1p+qz+rz3,FormalPowerSeries,z

n=0_α=RootOfr_Z3+q_Z+p_α1n3_α2r+qzn

(25)

convertsqrtπexpx1erfsqrtx,FPS,x=

n=0−1n2n!4n1xn+12n!

(26)

convertrandpolyz,degree=2log1+z,FormalPowerSeries

7z+167z22+n=02−1n42n2+163n+118zn+3n+1n+2n+3

(27)

convert1+z2expz+arcsechsqrtz,FormalPowerSeries

1+ln2lnz2+n=0n+1!n2+n+1!n22n22n+1!+n+1!zn+1n+1!2

(28)

References

  

Bertrand Teguia Tabuguia and Wolfram Koepf, Symbolic conversion of holonomic functions to hypergeometric type power series. Computer Algebra issue of the Journal of Programming and Computer Software, February 2022.

  

Bertrand Teguia Tabuguia, Power Series Representations of Hypergeometric Type and Non-Holonomic Functions in Computer Algebra. Ph.D. thesis, University of Kassel, Germany, May 2020.

  

Bertrand Teguia Tabuguia and Wolfram Koepf, Power series representations of hypergeometric type functions. In Corless R., Gerhard J., Kotsireas I. (eds): Maple in Mathematics Education and Research. MC 2020. Communications in Computer and Information Science, Springer, July 2021.

  

Bertrand Teguia Tabaguia, Wolfram Koepf, On the representation of non-holonomic power series, https://arxiv.org/abs/2109.09574, September 2021.

  

Koepf, Wolfram, Power series in computer algebra. Journal of Symbolic Computation, Vol. 13, (1992): 581-603.

Compatibility

• 

The output option was added in Maple 2022.

• 

The method=holonomic, method=quadratic, and method=default options were added in in Maple 2022.

• 

The differentialorder=N1,N2 option was added in Maple 2022.

• 

method=exponential and method=rational were workaround methods for when method=hypergeometric could not find a representation. However, this is no longer the case; the new method=hypergeometric always finds a better and a nicer power series representation. method=exponential or method=rational are still accepted and are equivalent to method=hypergeometric.

• 

The dir option is no longer supported. It is ignored when specified. As a workaround, one can make corresponding assumptions on the main variable.

• 

The makereal option is not necessary because the algorithm prioritizes real (rational) coefficients. It is ignored when specified.

• 

The convert/FormalPowerSeries command was updated in Maple 2022.

See Also

DEtools[FindODE]

gfun[holexprtodiffeq]

LREtools[hypergeomsols]

LREtools[mhypergeomsols]

RESol

series

Slode[mhypergeom_series_sol]