obtain the Taylor series for an ODE or PDE system - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Differential Equations : Rif : DEtools/rtaylor

DEtools[rtaylor] - obtain the Taylor series for an ODE or PDE system

Calling Sequence

rtaylor(solved, vars, options)

Parameters

solved

-

system in solved form

vars

-

(optional) solving variables of the system

options

-

(optional) sequence of options to specify the ranking for the solved form, initial data, and the order of the Taylor series

Description

• 

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 depvarindepvars=Taylor series.

Ranking Options

• 

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.

Other Considerations

• 

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.

Examples

withDEtools:

A simple ODE

rtaylorⅆ2ⅆx2fx=fx,order=4

fx=fx0+Dfx0xx012fx0xx0216Dfx0xx03+124fx0xx04

(1)

A PDE system with a single dependent variable

rtaylor2y2fx,y=xfx,yfx,y,2x2fx,y=2fx,y,order=3

fx,y=fx0,y0+D1fx0,y0xx0+D2fx0,y0yy0+fx0,y0xx02+D1,2fx0,y0xx0yy0+12D1fx0,y0fx0,y0yy02+13D1fx0,y0xx03+D2fx0,y0xx02yy0+122fx0,y02+D1fx0,y02xx0yy02+16D1,2fx0,y0fx0,y0+D2fx0,y0D1fx0,y0yy03

(2)

A PDE system with two dependent variables

rtaylor2x2fx,y=ygx,y,2y2fx,y=xgx,y,xgx,y=ygx,y

fx,y=fx0,y0+D1fx0,y0xx0+D2fx0,y0yy0+12D2gx0,y0xx02+D1,2fx0,y0xx0yy0+12D1gx0,y0yy02,gx,y=gx0,y0+D2gx0,y0xx0+D2gx0,y0yy0+12D2,2gx0,y0xx02+D2,2gx0,y0xx0yy0+12D2,2gx0,y0yy02

(3)

An example using initial data

sys:=2x2fx,y=0,2yxfx,y=0

sys:=2x2fx,y=0,2yxfx,y=0

(4)

id:=initialdatasys

id:=tableFinite=D1fx0,y0=_C1,Infinite=fx0,y=_F1y

(5)

rtaylorsys,id,order=3

fx,y=_F1y0+_C1xx0+D_F1y0yy0+12D2_F1y0yy02+16D3_F1y0yy03

(6)

An example using specified initial data and an expansion point

ids:=evalid_F1y=siny,_C1=1|evalid_F1y=siny,_C1=1

ids:=tableFinite=D1fx0,y0=1,Infinite=fx0,y=siny

(7)

rtaylorsys,ids,order=3,point=x=1,y=π

fx,y=x1y+π+16yπ3

(8)

See Also

caseplot, DEtools, DifferentialAlgebra[PowerSeriesSolution], Rif, rifsimp, rifsimp[nonlinear]


Download Help Document

Was this information helpful?



Please add your Comment (Optional)
E-mail Address (Optional)
What is ? This question helps us to combat spam