DEtools - Maple Help

Online Help

All Products    Maple    MapleSim


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

DEtools

  

initialdata

  

find initial data for a solved-form ODE or PDE system

 

Calling Sequence

Parameters

Description

Examples

Calling Sequence

initialdata(system,vars)

initialdata(rifresult,vars)

Parameters

system

-

list or set of polynomially nonlinear PDEs or ODEs in solved form

vars

-

(optional) dependent variables to compute the initial data with respect to

rifresult

-

single case result as returned by rifsimp

Description

• 

The initialdata command is part of the Rif subset of commands in the DEtools package.

• 

The initialdata command determines the initial data required to obtain a formal power series solution for a PDE or ODE system at a point (see rtaylor). The input system must be in standard or involutive form, so all equations must be of the form derivative = expression, and there must be no additional relations between the dependent variables that are not accounted for by the equations of the system. It can be used on any PDE or ODE system in isolated form, as well as for a single case result output by rifsimp. In general, all integrability conditions (equality of mixed partials) must be satisfied to obtain an accurate answer, so it is best to use rifsimp to simplify the system before calculation of the initial data.

• 

The optional vars argument specifies the variables for which to compute the initial data. This option is useful for classification problems, where different dependent variables play different roles in the system. This option also provides a way to specify the variables of a problem and their dependencies (in the event that one of the dependent variables does not appear in the Solved list, for example).

• 

The output is given in table form with the Finite entry containing the required zero dimensional initial data (constants), and the Infinite entry containing the required higher dimensional data (functions).

• 

When the input to initialdata is a rifresult, then additional information is provided on output. As well as the infinite and finite initial data, the resulting table may contain a Pivots and/or Constraint entry (if these entries are present in the rifresult). These returned lists of equations are given in terms of the initial data.

• 

When Constraint equations are present in the result, these equations must be satisfied by the initial data for the underlying system to be well posed. These Constraints must be accounted for manually.

• 

Notes:

  

When not using the rifresult form of the initialdata call, any nonlinear constraints must be converted to initial data form and imposed manually.

  

The initial data required for a system is independent of the right-hand side of the solved form equations, so for the following examples, ignore the right-hand side.

• 

Background Information: Formal Power Series and the initialdata algorithm

  

Being able to obtain a formal power series solution of a PDE or ODE system at a point essentially means that the values of all derivatives of the solution at that point are known, or can be found through differential relations from the system.

  

As a simple example, consider the ODE in ux, u'' = u u' (where prime notation is used to denote differentiation). Here the initial data is u(x0) and u'(x0). With this data and the ODE, you can then determine u''(x0) = u(x0) u'(x0), and from the derivative of the ODE, you can obtain u'''(x0) = u'(x0)^2+u(x0) u''(x0), and so on.

  

Knowledge of all derivatives then allows you to write u(x) = u(x0)+(x-x0) u'(x0)+(x-x0)^2 u(x0) u'(x0)/2+`...` to any desired order.

  

As an example of the initialdata algorithm in two dimensions, consider the following system.

sys:=[diff(f(x,y),x,x,x)=0,diff(f(x,y),x,y)=0];

sys:=3x3fx,y=0,2yxfx,y=0

(1)

DEtools[initialdata](sys);

tableInfinite=fx0,y=_F1y,Finite=D1fx0,y0=_C1,D1,1fx0,y0=_C2

(2)
  

In this case, we have one arbitrary function and two arbitrary constants. Plots can be used to visualize how the initialdata algorithm calculates this data. Consider an integer-valued graph where the x-axis represents differentiation by x, and the y-axis represents differentiation by y. Looking at f[xxx] from the original system, we can place dots at that value and at all of its derivatives:

d1:=[3,0],[4,0],[3,1],[4,1],[3,2],[4,2],[3,3],[4,3],[3,4],[4,4]:

plot([d1], x=-0.1..4.1, y=-0.1..4.1, style=point, tickmarks=[[0,1,2,3,4],[0,1,2,3,4]]);

  

Note, for example, that the point (4,1) corresponds to f[xxxxy], which can be obtained from the differentiation of f[xxx] with respect to x and y. Adding all derivatives that are known from the f[xy] equation gives the following graph.

d2:=[1,1],[2,1],[3,1],[4,1],[1,2],[2,2],[3,2],[4,2],[1,3],[2,3],[3,3],[4,3],[1,4],[2,4],[3,4],[4,4]:

plot([d1,d2], x=-0.1..4.1, y=-0.1..4.1, style=point, tickmarks=[[0,1,2,3,4],[0,1,2,3,4]]);

  

The objective is to fill in all missing points. From the graph, we can see that all points of the form (0,i) are missing, which explains the arbitrary function of y in the answer. This function describes the data for all y derivatives of f. We still require the points (1,0) and (2,0), which correspond to the arbitrary constants present in the answer for the values of f[x] and f[xx], as given in the Finite list.

Examples

withDEtools:

sys12x2fx,y,z=0,yfx,y,z=0,zfx,y,z=0

sys1:=2x2fx,y,z=0,yfx,y,z=0,zfx,y,z=0

(3)

This system is fully specified if the following derivatives are known.

initialdatasys1

tableInfinite=,Finite=fx0,y0,z0=_C1,D1fx0,y0,z0=_C2

(4)

For this next system, the initial data contains a number of arbitrary functions of one variable.

sys22yxfx,y,z=0,2zxfx,y,z=0,2zyfx,y,z=0

sys2:=2yxfx,y,z=0,2zxfx,y,z=0,2zyfx,y,z=0

(5)

initialdatasys2

tableInfinite=fx,y0,z0=_F1x,D2fx0,y,z0=_F2y,D3fx0,y0,z=_F3z,Finite=

(6)

For the next system, the initial data contains a number of arbitrary functions of two variables.

sys33zyxfx,y,z=0

sys3:=3zyxfx,y,z=0

(7)

initialdatasys3

tableInfinite=fx,y,z0=_F1x,y,D3fx,y0,z=_F2x,z,D2,3fx0,y,z=_F3y,z,Finite=

(8)

The one-dimensional heat equation.

sys4tux,t=2x2ux,t

sys4:=tux,t=2x2ux,t

(9)

initialdatasys4

tableInfinite=ux,t0=_F1x,Finite=

(10)

This example is a system that contains mixed data.

sys53y2xfx,y,z=0,2zxfx,y,z=0,3z2yfx,y,z=0

sys5:=3y2xfx,y,z=0,2zxfx,y,z=0,3z2yfx,y,z=0

(11)

initialdatasys5

tableInfinite=fx,y0,z0=_F1x,D2fx,y0,z0=_F2x,D3fx0,y0,z=_F3z,D2,2fx0,y,z0=_F4y,D2,2,3fx0,y,z0=_F5y,Finite=D2,3fx0,y0,z0=_C1

(12)

An example using initialdata with a nonlinear rif result,

sys6ⅆ2ⅆx2ux2ux+ⅆⅆxux3

sys6:=ⅆ2ⅆx2ux2ux+ⅆⅆxux3

(13)

rif6rifsimpsys6

rif6:=tablePivots=ⅆ2ⅆx2ux0,Case=ux0,ⅆ2ⅆx2ux2,ⅆ2ⅆx2ux0,ⅆ3ⅆx3ux,Solved=ⅆ3ⅆx3ux=12ⅆ2ⅆx2uxⅆⅆxux3ⅆⅆxux2ux,Constraint=ⅆ2ⅆx2ux2ux+ⅆⅆxux3=0

(14)

initialdatarif6

tablePivots=_C30,Infinite=,Constraint=_C1_C32+_C23=0,Finite=ux0=_C1,Dux0=_C2,D2ux0=_C3

(15)

where it is understood that the selected initial data must obey the nonlinear constraint (so in truth, you have only two free parameters) and must not violate the pivots (so _C30).

See Also

caseplot

DEtools

DifferentialAlgebra[PowerSeriesSolution]

Rif

rifsimp

rtaylor

 


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