computes polynomial solutions for (differential or not) system of equations. - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Differential Equations : PDEtools : Symmetries : PDEtools/PolynomialSolutions

PDEtools[PolynomialSolutions] - computes polynomial solutions for (differential or not) system of equations.

Calling Sequence

PolynomialSolutions(S, PDESYS, DepVars)

Parameters

PDESYS

-

an equation or a set or list of them; it can include PDEs, ODEs, non-differential equations and inequations

DepVars

-

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

degree = ...

-

optional - to indicate the upper bound degree of the solution, could be a non-negative integer, or a set of equations where the lhs is one of unknowns (DepVars) and the rhs is the upper bound degree for its solution

dependency = ...

-

optional - to indicate the dependency of the solution, could be a non-negative integer, a range of them, a name, or a set of names, or a set or list containing any of the previous objects

HINT = ...

-

optional - a set of equations with unknowns in the lhs and related hints polynomial in the independent variables in the rhs

ivars = ...

-

optional, a name or set of them indicating additional independent variables, typically not present in the dependency of the unknowns DepVars

numberofsolutions = ...

-

optional - can be a nonnegative integer or the keyword all, to specify the number of polynomial solutions desired; default is all having degree as upper bound

simplifier = ...

-

optional - indicates the simplifier to be used instead of the default simplify/size

Description

• 

Given PDESYS, a system of equations for DepVars, possibly including inequations and non-differential equations, PolynomialSolutions computes solutions that are polynomial in the independent variables of the problem. The independent variables are typically the variables DepVars depend on, plus any others optionally specified with ivars = .... PolynomialSolutions also works with anticommutative variables set using the Physics package using the approach explained in PerformOnAnticommutativeSystem.

• 

If DepVars is not given, PolynomialSolutions will consider all the differentiated unknown functions in PDESYS as unknown of the problems. Specifying DepVars however permits not only restricting the unknowns in different ways but also specifying unknowns of the problems which do not appear differentiated in PDESYS.

• 

PolynomialSolutions uses a heuristic algorithm to compute upper bound degrees. The algorithm used extends, for non-autonomous systems of equations, the ideas used in TWSolutions (step 2. of The solving process) for computing solutions to autonomous systems. This extended algorithm is accessible at user level as PDEtools:-Library:-UpperBounds (see PDEtools,Library). To override these upper bounds see the degree option.

• 

PolynomialSolutions is automatically invoked by Infinitesimals, InvariantSolutions and SimilaritySolutions when they receive the optional argument typeofsymmetry = polynomial in which case the dependency and degree can also be specified directly to those commands.

• 

To avoid having to remember the optional keywords, if you type the keyword misspelled, or just a portion of it, a matching against the correct keywords is performed, and when there is only one match, the input is automatically corrected.

Examples

withPDEtools:

declareξ,ηx,y

ξx,ywill now be displayed asξ

ηx,ywill now be displayed asη

(1)

sys1:=xξx,yxηx,y+yξx,yyηx,y+ηx,y2x2ξx,y+2y2ξx,y=1

sys1:=ξxηx+ξyηy+ηξx,x+ξy,y=1

(2)

sol1:=PolynomialSolutionssys1

sol1:=η=_C3y2+_C2y+_C1_C5x_C5,ξ=_C5x+_C4,η=_C3x2+_C2x+_C1_C52+2_C3_C6x_C2_C61y_C5+_C3_C62y2_C52,ξ=_C5y+_C6x+_C4

(3)

Solutions can be tested as usual, with pdetest

mappdetest,sol1,sys1

0,0

(4)

Compute polynomial solutions of degree = 1 depending on 1 variable

PolynomialSolutionssys1,degree=1,dependency=1

η=_C1_C3x_C3,ξ=_C3x+_C2,η=_C1_C3y_C3,ξ=_C3y+_C2

(5)

Compute polynomial solutions depending on 0 (there are none) 1 and 2 variable

PolynomialSolutionssys1,dependency=0..2

η=_C1_C3x_C3,ξ=_C3x+_C2,η=_C1_C3y_C3,ξ=_C3y+_C2,η=_C3y2+_C2y+_C1_C5x_C5,ξ=_C5x+_C4,η=_C3x2+_C2x+_C1_C52+2_C3_C6x_C2_C61y_C5+_C3_C62y2_C52,ξ=_C5y+_C6x+_C4

(6)

PolynomialSolutions is invoked by InvariantSolutions and other commands when they receive the optional argument typeofsymmetry = polynomial. In this example, solutions are computed where the degree of the polynomial infinitesimals is requested to be zero

InvariantSolutionssys1,typeofsymmetry=polynomial,degree=0

η=x+_C1_F1x,ξ=_F1x,η=y+_C1_F1y,ξ=_F1y

(7)

An example involving mathematical functions (ln):

sys2:=Ix32y2ξx,yy+3yξx,yy,Ix22y2ηx,yxy2+22yxξx,yxy23ηx,yx+3xyηx,yy+2yξx,yy2y2,x2I2yxηx,yx2yI2x2ξx,yx2yIxxξx,yy6Ix2xηx,y+6yξx,yy4lnx+3yξx,yy2+ξx,yyIy,4ξx,yy3lnx2ξx,yy+6ηx,yxlnxy2Ix2xηx,y+4xxξx,yy3lnx+2xxξx,yy2xyηx,yy3lnxxyηx,yy+2ξx,yy3+2x2ηx,yx3I+ηx,yx

sys2:=Ix3yξy,y+3ξyy,Ix2xy2ηy,y+2xy2ξx,y+3xyηy+2y2ξy3ηxy2,x2Iηx,yx2yIξx,xx2yIxξxy6Ix2ηx+6ξyy4lnx+3ξyy2+Iξyy,4ξy3lnx2ξy+6ηxlnxy2Ix2ηx+4xξxy3lnx+2xξxy2xηyy3lnxxηyy+2ξy3+Iηx,xx3+ηx

(8)

The internal routine computes the following upper bound degrees

PDEtools:-Library:-UpperBoundssys2,ξ,ηx,y

η=3,ξ=2

(9)

When the degrees are not specified, PolynomialSolutions automatically calls that internal routine to decide about the degrees, resulting in the following solution

sol2:=PolynomialSolutionssys2

sol2:=η=y3,ξ=x

(10)

pdetestsol2,sys2

0,0,0,0

(11)

See Also

DEtools/polysols, pdetest, PDEtools, PDEtools/Library, PerformOnAnticommutativeSystem, Physics, Setup, TWSolutions


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