ODEs and PDEs - Maple Help

Online Help

All Products    Maple    MapleSim


Ordinary and Partial Differential Equations

Maple is the world leader in finding exact solutions to ordinary and partial differential equations. Maple 2021 extends that lead even further with new algorithms and techniques for solving more ODEs and PDEs.

For Maple 2021, there are significant improvements in dsolve for the exact solution of 2ndorder linear ODEs using hypergeometric functions. The algorithms implemented are at the frontier of the understanding of this problem, and handle classes of extended equations with apparent singularities as well as the most common linear equations with 4 and 5 singularities.

For ODEs and PDEs, the LieAlgebrasOfVectorFields package in Maple 2021 has a new command MapDE , for analyzing the possible linearization of polynomially nonlinear equations, and determining the corresponding mapping when it exists, allowing in that way for the computation of more approximate and exact solutions.

The new Student[ODEs] package covers the material in a standard first course in ODEs and provides step-by-step tools for solving ODEs as well as interactive visualization.  For details, see Student Packages.

 

Advanced type of Hypergeometric solutions for 2nd order linear ODEs

Linearization of polynomially nonlinear ODE and PDE

Advanced type of Hypergeometric solutions for 2nd order linear ODEs

When they exist, 2F1, 1F1 and 0F1 hypergeometric solutions for 2ndorder linear ODEs are now computable in a rather general case. An equation with 3 regular singularities and any number of apparent singularities can now be solved in terms of 2F1 equations using the algorithms presented in (6) and (7). Likewise, an equation with 2 or 1 singularities, one of which is irregular, can now be solved in general using 1F1 and 0F1 hypergeometric functions using the algorithms presented in (1), (2) and (3). Common equations with up to 5 regularities (plus any number of apparent singularities) can now be solved using mappings and the tables from references (4) and (5). All this significantly extends the previous capabilities in solving linear equations, and also of higher order and nonlinear ODE and PDE that through a reduction or symmetry process require the solving of intermediate 2ndorder linear ODEs.

 

The new algorithms are now automatically used by dsolve and also accessible through a new DEtools command, hypergeometricsols.

 

Examples

Despite the powerful algorithms for hypergeometric solutions of previous Maple releases, none of the following examples could be solved. Now they can. This equation admits 2F1 solutions where the argument is polynomial in x of degree 3

ode__1ⅆ2ⅆx2yx=1296a56x2+4a3+109x+4a949xx1b+216x288432x231+bx1x2x1ⅆⅆxyx9a56a16x232b41+bx1x2x1 yx:

dsolveode__1

yx=_C1F12112a2,a2+512;12;bx1x2+_C2x1xF12712a2,1112a2;32;bx1x2

(1.1)

Verify this solution

odetest,ode__1

0

(1.2)

This other equation admits solutions in terms of modified Bessel functions of the first and second kind

ode__2ⅆ2ⅆx2yx=8x482x92x15ⅆⅆxyx+4x696x5+931x44452x3+10539x212474x+97204x2x322x92x15 yx:

dsolveode__2

yx=_C1x4x336x2+141x225I1x2+4x3+24x263x+27I2x2x3+_C2x4x3+36x2141x+225K1x2+4x3+24x263x+27K2x2x3

(1.3)

Moreover, note that each independent solution involves a linear combination of functions. Indeed, while Bessel (or more generally, 0F1) solutions are related to equations with 1 singularity of irregular kind, this example has other regular singularities:

DEtoolssingularitiesode__2

regular=0,3,92,152,irregular=

(1.4)

thus requiring the use of a linear combination of pFq functions to construct a solution. A similar situation happens with the next example where, due to the presence of apparent singularities, a linear combination of - this time 1F1 - hypergeometric functions (KummerM and KummerU) is required to solve the problem

ode__3ⅆ2ⅆx2yx=16x4+32x236x+816x416x3+32x24xⅆⅆxyx12x3+13x2+x9416x416x3+32x24 x yx:

dsolveode__3

yx=_C112x122M54,3,x+20x+5M94,3,x+_C216x122U54,3,x+20x5U94,3,x

(1.5)

In addition to computing new solutions out of reach in previous releases, when the linear ODE is of Heun type, so it has 4 regular singularities or one of its confluent cases, and it happens to be one of the special Heun function cases that can be expressible using hypergeometric functions, both kinds of solutions are now computable. This equation is of the Heun triconfluent type

ode__4  ⅆ2ⅆx2yx = 4 x44 x2+10 x2yx:

The algorithms in dsolve perceive that and solve it accordingly

dsolveode__4

yx=_C1ⅇ23x3+xHT36234,152,613,x6233+_C2ⅇ23x3xHT36234,152,613,x6233

(1.6)

By using the option of indicating the method, a solution in terms of linear combinations of Airy functions (of the 0F1)  is now also computable

dsolveode__4,hypergeometricsols

yx=_C12x+1Aix21+2x2+x1Aix21+_C22x+1Bix21+2x2+x1Bix21

(1.7)

To express this solution in terms of the more general 0F1 form you can use

convert,hypergeom

yx=_C12x+13316Γ232F10;13;x2+139π+x212313F10;53;x2+1396Γ23+2x2+x1313F10;23;x2+1393Γ23x21316Γ23F10;43;x2+1392π+_C22x+13323Γ232F10;13;x2139π+356x12x+12F10;53;x2139x1x+16x26Γ23+2x2+x1356F10;23;x2+1393Γ23+323x21Γ23F10;43;x2+1392π

(1.8)

A similar situation, this time with regards to this equation of the Heun biconfluent class

ode__5ⅆ2ⅆx2yx=x2+3x+43xⅆⅆxyx+x+19x yx:

dsolveode__5

yx=_C1HB13,6,−3,1469,6x6ⅇxx+66+_C2HB13,6,−3,1469,6x6ⅇxx+66x13

(1.9)

The solution is expressible in terms of linear combinations of modified Bessel functions of the 1st and 2nd kinds

dsolveode__5,hypergeometricsols

yx=_C1ⅇxx+612x2+6x+16I134+x4+xx12+4+x4+xxI434+x4+xx124+xx16+_C2ⅇxx+612x2+6x+16K134+x4+xx124+x4+xxK434+x4+xx124+xx16

(1.10)

The following example is of the same kind but more general: it belongs to the Heun general class, but again it is one of those special cases where pFq function solutions exist, in this case 2F1 with rational coefficients

ode__6ⅆ2ⅆx2yx=345x2+14x+14x27x2+14x+3ⅆⅆxyx+3x+12b3b+364x227x2+14x+3 yx:

dsolveode__6

yx=_C1HG4I2+74I2+7,45927b+4I2b143