obtain the Taylor series for an ODE or PDE system
rtaylor(solved, vars, options)
system in solved form
(optional) solving variables of the system
(optional) sequence of options to specify the ranking for the solved form, initial data, and the order of the Taylor series
The rtaylor command is part of the Rif subset of commands in the DEtools package.
The rtaylor command uses an output rifsimp form to obtain local Taylor series expansions for all dependent variables in the ODE or PDE system simultaneously. The Taylor series output is a list containing equations of the form depvar⁡indepvars=Taylor series.
The ranking related options that are accepted by rtaylor include the specification of the vars as a ranking, and the ranking and indep options described in rifsimp[ranking], rifsimp[options], and rifsimp[adv_options].
Note: Specification of vars different from those used to obtain the result from rifsimp can give incomplete results.
Note: Specification of vars as a set of dependent variables, or as a single dependent variable, makes the use of ranking options unpredictable (i.e. if using the vars as part of the ranking, it must be ordered, so must be a list).
Initial Data Options
The table resulting from a call to initialdata can be given as an option to rtaylor, in which case the Taylor series is given in terms of the functions present in the initial data.
In addition, the arbitrary functions and constants on the right-hand sides of the specified initial data can be given specific values, and the expansion can be computed for these values.
Taylor Series Order Options
The order=n option specifies the order that the Taylor series should be computed to, and must be a non-negative integer. The default value is 2.
The point=[ivar1=val1,...] option specifies an expansion point for the series. When this option is used, every independent variable must be given a value.
In general, any Constraint or DiffConstraint relations (see rifsimp[nonlinear]) in the rif form cannot be used in an automatic way, so they are ignored. These relations must be accounted for manually after the Taylor series calculation. Special care must be taken when DiffConstraint relations are present, because all derivatives of these relations must be manually accounted for. This is not the case for Constraint relations, as they are purely algebraic.
The requirement that the input solved form be in rif form can be relaxed mildly, but rtaylor still requires that the equations are in a valid solved form that matches the input ranking (given in the options), and have no integrability conditions remaining. Only when these conditions hold is the resulting Taylor series an accurate representation of the local solution.
A simple ODE
A PDE system with a single dependent variable
A PDE system with two dependent variables
An example using initial data
An example using specified initial data and an expansion point
Download Help Document
What kind of issue would you like to report? (Optional)