convert/Sum - convert to inert sum representation - a formal power series
convert/sum - convert to sum representation
|
Calling Sequence
|
|
convert( expr, Sum, opts, Sum_opts );
convert( expr, sum, opts, Sum_opts );
|
|
Parameters
|
|
expr
|
-
|
valid Maple expression, equation, set, list, or similar structure
|
opts
|
-
|
(optional) arguments as described in convert/to_special_function
|
Sum_opts
|
-
|
(optional) any of the following arguments intrinsic to a Sum:
|
expansionvariable = x
|
-
|
indicates the expansion variable - the rhs could be any name
|
dummy = n
|
-
|
indicates the dummy (summation) variable - the rhs could be any name
|
x = a
|
-
|
the lhs and rhs respectively indicate the expansion variable and expansion point
|
b(k)
|
-
|
a function name with an expansion variable as argument, representing the "kth coefficient" of the series expansion
|
method = value
|
-
|
indicates the method to be used; <value> can be any of: Local, Global, rational, exponential, hypergeometric
|
makereal = value
|
-
|
indicates that a series with real coefficients should be returned; <value> can be any of: true, false
|
differentialorder = n
|
-
|
indicates an upper bound for the order the intermediate differential equation used to represent expr; <value> can be any positive integer
|
dir = value
|
-
|
indicates the direction of the limit computation for initial values; <value> can be any of: default, left, right, real, or complex
|
recurrence = value
|
-
|
indicates that the output should be a recurrence for b(k) when a closed form for these coefficients, to represent expr as a Sum, cannot be computed
|
|
|
|
|
Description
|
|
•
|
All the optional arguments accepted by the Maple convert/to_special_function network are allowed, so that one can restrict the application of convert to Sum in different manners.
|
•
|
The convert(expr, sum) command does the same as convert(expr, Sum) followed by an attempt to compute all the Sums introduced, and here again to override the default behavior you can use any of the optional arguments described in convert/to_special_function.
|
•
|
By default, a combination of the Global and Local methods is used; in this approach, when possible, all mathematical functions found expr are converted, and in cases of nested mathematical functions, or linear combination of functions which admit a hypergeometric representation, a Global conversion is attempted first, taking advantage of the properties of such expressions. By invoking this conversion with the optional argument method = Global (or skip = Local), or with method = Local (or skip = Global), respectively only the Global or the Local methods are used.
|
|
NOTE: The use of the following optional arguments automatically make convert/Sum to skip the Local method and only work using the Global method: x = a, makereal = true, differentialorder = n, dir = value, method = value, indicating one of the methods rational, exponential, hypergeometric, and recurrence = value. Also, powers of x are not converted unless you explicitly request that using the optional argument include = powers (all powers) or include = radicals (only fractional powers).
|
•
|
The types of expressions that Global conversions can handle are: - expressions of hypergeometric type, where b(k+m)/b(k) is a rational function of k for some integer m; - expressions of exponential type, which satisfy a linear homogeneous differential equation with constant coefficients; - expressions of rational type, which are either rational or have a rational derivative; - linear combinations of hypergeometric functions are treated by the Petkovsek-van-Hoeij algorithm; see LREtools[hypergeomsols].
|
|
|
Optional arguments for Global conversions
|
|
•
|
differentialorder: a positive integer n (default: n=4); upper bound for the order of the differential equation searched for. This controls the depth of the search for a differential equation for expr. Higher values of n will increase the chance to find the solution, but increase the running time as well.
|
•
|
dir: one of default, left, right, real, or complex; direction of the limit computation for initial values. If a, the expansion point, is finite, then the default is dir = complex. If a is either infinity or -infinity, then the default is dir = real. See also limit. This optional argument can be used to request real or one-sided (e.g. asymptotic) series.
|
•
|
makereal: either true or false (default); makereal = true (or just makereal for short) indicates that a series with real coefficients should be returned.
|
•
|
method: one of default, hypergeometric, rational, or exponential. Specifies the method that will be used; the default method uses an internal selection strategy. See the Examples below for an illustration of the various methods.
|
•
|
recurrence: either true or false (default). If recurrence = true (or recurrence for short) is given and the coefficients entering the formal power series can be computed in closed form, then the output is a recurrence for b(k).
|
|
|
Examples
|
|
>
|
|
The 7 functions in the "Ei_related" class are:
| |
| (1) |
>
|
|
| (2) |
>
|
|
| (3) |
>
|
|
| (4) |
Note that some conversions are possible only under assumptions.
>
|
|
* Partial match of "sum" against topic "sum_form".
| |
| (5) |
Therefore, the following is not converted.
>
|
|
| (6) |
The conversion occurs only when you specify the proper assumptions.
>
|
|
| (7) |
>
|
|
| (8) |
By default, powers, or rational expressions, are not converted to Sum unless explicitly requested, e.g. using the include = powers (could also be include = ratpoly or include = radicals) optional argument; make the dummy summation be
>
|
|
| (9) |
>
|
|
| (10) |
>
|
|
| (11) |
The output above is computed using the Global algorithm; let's see the output representing first, Locally, each of the powers found in as MeijerG functions, then representing these functions as Sums
>
|
|
| (12) |
Expansions around t = 0 or t = 1; note the indication of the expansion variable when indicating the expansion point
>
|
|
| (13) |
>
|
|
| (14) |
>
|
|
| (15) |
>
|
|
| (16) |
For composite mathematical functions, first the Global approach is attempted
>
|
|
| (17) |
>
|
|
| (18) |
The following examples illustrate the use of the method 'rational', 'exponential', 'hypergeometric'.
>
|
|
| (19) |
>
|
|
| (20) |
>
|
|
| (21) |
>
|
|
| (22) |
>
|
|
| (23) |
>
|
|
| (24) |
User-defined functions are handled provided their derivative is known. You define the derivative of the function as follows (see diff for more information).
>
|
`diff/g` := proc(a,x) g(a)*diff(a,x) end proc;
|
| (25) |
>
|
|
| (26) |
Indefinite integrals are handled.
>
|
|
| (27) |
>
|
|
| (28) |
Linear combinations of hypergeometric functions are recognized; compare the default output with the output requesting makereal = true
>
|
|
| (29) |
>
|
|
| (30) |
>
|
|
| (31) |
In the next example, the output is expressed in terms of algebraic numbers.
>
|
|
| (32) |
>
|
|
| (33) |
Maple's special functions are handled.
>
|
|
| (34) |
>
|
|
| (35) |
In this case, convert/Sum returns without expanding
>
|
|
| (36) |
>
|
|
| (37) |
because a conversion exists only for restricted values of the parameters; these restrictions are seen via (see FunctionAdvisor):
>
|
|
* Partial match of "sum" against topic "sum_form".
| |
| (38) |
>
|
|
| (39) |
>
|
|
| (40) |
>
|
|
| (41) |
>
|
|
* Partial match of "sum" against topic "sum_form".
| |
| (42) |
>
|
|
| (43) |
Hidden polynomials are frequently detected.
>
|
|
| (44) |
>
|
|
| (45) |
Asymptotic power series can be computed.
>
|
|
| (46) |
>
|
|
| (47) |
>
|
|
| (48) |
>
|
|
| (49) |
Real and one-sided (asymptotic) series can be computed using the dir option.
>
|
|
| (50) |
>
|
|
| (51) |
>
|
|
| (52) |
>
|
|
| (53) |
>
|
|
| (54) |
An example where a the command fails to compute a closed form for the coefficients in the formal power series and so it returns, under request, the recurrence equation for these coefficients.
>
|
|
| (55) |
>
|
|
| (56) |
>
|
|
| (57) |
Note however that in this case a conversion exists for the appropriate conditions on the function parameters:
>
|
|
* Partial match of "sum" against topic "sum_form".
| |
| (58) |
>
|
|
| (59) |
|
|
See Also
|
|
assuming, convert, convert/FormalPowerSeries, convert/to_special_function, dsolve,formal_series, dsolve,formal_solution, FunctionAdvisor, gfun,holexprtodiffeq, map, Slode,FPseries, Sum, sum
|
|
References
|
|
|
[1] Wolfram Koepf, Power series in computer algebra. Journal of Symbolic Computation 13, 1992, 581-603
|
|
[2] Dominik Gruntz, Wolfram Koepf: Maple package on formal power series. Maple Technical Newsletter 2 (2), 1995, 22-28
|
|
|