computes formal power series solutions for differential equation systems - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Differential Equations : pdsolve : pdsolve/series

pdsolve/series - computes formal power series solutions for differential equation systems

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

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,t,ux,0 or u0,0 when the unknown is ux,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

pde1:=2t2ux,txux,t=0

pde1:=2t2ux,txux,t=0

(1)

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

iv1:=ux,0=fx,D2ux,0=gx

iv1:=ux,0=fx,D2ux,0=gx

(2)

pdsolvepde1,iv1,series,order=2

ux,t=12t2+2xDf0+Dg0xt+12D2f0x2+g0t+f0

(3)

A general nonlinear PDE system

pde2:=2x2ux,t2vx,t+2x2ux,tvx,t+xux,t=0,2txux,t=0,2t2ux,t21=0

pde2:=2x2ux,t2vx,t+2x2ux,tvx,t+xux,t=0,2xtux,t=0,2t2ux,t21=0

(4)

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

iv2:=vx,0=1+fx,D2,2u0,0=1,u0,0=c0,D1u0,0=0,D2u0,0=c2,D1,1u0,0=1

iv2:=vx,0=1+fx,D2,2u0,0=1,u0,0=c0,D1u0,0=0,D2u0,0=c2,D1,1u0,0=1

(5)

pdsolvepde2,iv2,u,vx,t,series,order=3

ux,t=c0+c2t12t2,vx,t=1+f0+Df0x+12D2f0x2+16D3f0x3

(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

pde3:=2x2ux,t=ⅇsinxtux,t

pde3:=2x2ux,t=ⅇsinxtux,t

(7)

pdsolvepde3,x=x0,t=t0,series,order=2

ux,t=12xx02ⅇsinx0+2t2t0D2ux0,t0+D1,2ux0,t0xx0tt0+12D2,2ux0,t0tt02+122x2x0D1ux0,t0+ux0,t0

(8)

pde4:=2x2ux,t=BesselJν,sinxtux,t

pde4:=2x2ux,t=BesselJν,sinxtux,t

(9)

pdsolvepde4,x=x0,series,order=2

ux,t=12xx02BesselJν,sinx0+2tD2ux0,0+D1,2ux0,0xx0t+122x2x0D1ux0,0+12D2,2ux0,0t2+ux0,0

(10)

See Also

pdsolve, pdsolve/boundaryconditions, DifferentialAlgebra, PowerSeriesSolution


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