series - Maple Help
For the best experience, we recommend viewing online help using Google Chrome or Microsoft Edge.

Online Help

All Products    Maple    MapleSim


pdsolve/series

computes formal power series solutions for differential equation systems

 

Calling Sequence

Parameters

Options

Description

Examples

Calling Sequence

pdsolve(DE, F, series, opts)

Parameters

DE

-

a system of differential equations (ODEs, PDEs and/or algebraic constraints), possibly including initial values for the unknowns, their derivatives and/or an expansion point

F

-

optional, may be required, a function, or a set or list of them, indicating the dependent variables of the problem

disregard = ...

-

optional, to disregard whether the initial values cancel the equations and/or inequations of DE, the right-hand-side can be: nothing (default), inequations, equations, all

order = n

-

optional, a nonnegative integer specifying the order of the series expansion, default value is Order

Options

• 

The opts arguments may contain one or more of the following options.

• 

disregard = <nothing (default), inequations, equations, all>. This option permits to restrict the test for the initialvalues (see above).

– 

equations: the equations are disregarded; pdsolve only checks that the inequations are satisfied by initialvalues.

– 

inequations: the inequations are disregarded; pdsolve only checks that the equations are satisfied by initialvalues.

– 

all: equations and inequations are disregarded.

• 

order = n. This option specifies the truncation order for the series expansion; if not given, the value of n used is that of Order.

Description

• 

The function call pdsolve(DE, series) returns an integral (exponents are nonnegative integers), formal (the convergence issue is not addressed), power series solution of the differential ideal defined by DE, truncated at order = n.

• 

The unknowns may or not be specified. When the unknowns are specified (as a function, or a set or list of them), only these ones specified are considered unknowns and any other function is considered arbitrary. The solution is then expected to be valid for arbitrary values of these arbitrary functions. On the other hand, when the unknowns are not specified, all functions - differentiated or not - are considered unknowns.

• 

The initial values for the unknowns, and/or their derivatives, and/or an expansion point, can be provided in DE by enclosing all the arguments into a set or a list. The set of initial values must have only one expansion point and must be linear in the unknowns and derivatives evaluated at the expansion point. Typically, these initial values are specified as equations of the form u0=value_of_u0, where the u0 are dependent variables and/or their derivatives evaluated at the expansion point - say u0&comma;t,ux&comma;0 or u0&comma;0 when the unknown is ux&comma;t - and value_of_u0 represent their values, possibly being expressions involving independent variables (in which case keep the independent variable as is, in the left-hand-side u0).

• 

If initial values are not specified, the expansion point is the origin, all the independent variables equal to 0, and all the functions and their derivatives different from zero. It is also possible to specify only an expansion point different than zero directly in DE - say including equations of the form x=x0,t=t0.

• 

The values attributed to the unknowns and their derivatives at the expansion point, also cannot be chosen freely. They must satisfy the system F=0,S0, returned by the DifferentialAlgebra command PowerSeriesSolution when using the option conditions. These conditions, which are checked by default, can be disregarded using the option disregard = ...; see the Examples section.

– 

If the attributed values are solutions of the system F=0,S0, then, the returned power series are truncated solutions of the differential equation system DE.

– 

If these values do not satisfy F=0, then, the returned power series are not truncated solutions of the differential ideal.

– 

If they satisfy F=0, but do not satisfy S0, then, the returned power series may, or may not, be truncated solutions of the of DE. If they are, they do not need to be unique.

Examples

The heat equation in two dimensions

pde1diffux&comma;t&comma;t&comma;tdiffux&comma;t&comma;x=0

pde12t2ux&comma;txux&comma;t=0

(1)

Classically, one says that the solutions of the heat equation depend on two arbitrary functions ux&comma;0=fx and utx&comma;0=gx. Our choice for initial values is then

iv1ux&comma;0=fx,D2ux&comma;0=gx

iv1ux&comma;0=fx,D2ux&comma;0=gx

(2)

pdsolvepde1&comma;iv1&comma;series&comma;order=2

ux&comma;t=t2+2xDf02+Dg0xt+D2f0x22+g0t+f0

(3)

A general nonlinear PDE system

pde2diffux&comma;t&comma;x&comma;x2vx&comma;t+diffux&comma;t&comma;x&comma;xvx&comma;t+diffux&comma;t&comma;x=0,diffux&comma;t&comma;x&comma;t=0,diffux&comma;t&comma;t&comma;t21=0

pde22x2ux&comma;t2vx&comma;t+2x2ux&comma;tvx&comma;t+xux&comma;t=0,2txux&comma;t=0,2t2ux&comma;t21=0

(4)

The initial values are chosen to satisfy the conditions that can be computed using PowerSeriesSolution with the option condition

iv2vx&comma;0=1+fx,D2,2u0&comma;0=1,u0&comma;0=c0,D1u0&comma;0=0,D2u0&comma;0=c2,D1,1u0&comma;0=1

iv2vx&comma;0=1+fx,D2,2u0&comma;0=1,u0&comma;0=c0,D1u0&comma;0=0,D2u0&comma;0=c2,D1,1u0&comma;0=−1

(5)

pdsolvepde2&comma;iv2&comma;u&comma;vx&comma;t&comma;series&comma;order=3

ux&comma;t=c0+c2t12t2&comma;vx&comma;t=1+f0+Df0x+D2f0x22+D3f0x36

(6)

Note than unlike PowerSeriesSolution, pdsolve can handle any mathematical function (or compositions of them) in the coefficients of the differential equations provided that it admits a differential equation representation, for example

pde3diffux&comma;t&comma;x&comma;x=expsinxdiffux&comma;t&comma;t

pde32x2ux&comma;t=&ExponentialE;sinxtux&comma;t

(7)

pdsolvepde3&comma;t=t0&comma;x=x0&comma;series&comma;order=2

ux&comma;t=xx02&ExponentialE;sinx0+2t2t0D2ux0&comma;t02+D1,2ux0&comma;t0xx0tt0+D2,2ux0&comma;t0tt022+D1ux0&comma;t0xx0+ux0&comma;t0

(8)

pde4diffux&comma;t&comma;x&comma;x=BesselJν&comma;sinxdiffux&comma;t&comma;t

pde42x2ux&comma;t=BesselJν&comma;sinxtux&comma;t

(9)

pdsolvepde4&comma;x=x0&comma;series&comma;order=2

ux&comma;t=xx02BesselJν&comma;sinx0+2tD2ux0&comma;02+D1,2ux0&comma;0xx0t+D2,2ux0&comma;0t22+D1ux0&comma;0xx0+ux0&comma;0

(10)

See Also

pdsolve

pdsolve/boundaryconditions

DifferentialAlgebra

PowerSeriesSolution