DEtools/IVPsol - Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : DEtools/IVPsol

DEtools

  

IVPsol

  

specializes ODE solutions to match given initial or boundary conditions

 

Calling Sequence

Parameters

Description

Examples

Calling Sequence

IVPsol(ics, sol)

Parameters

ics

-

a set or list of initial or boundary conditions for an ODE problem. This value can include the ODE itself.

sol

-

a solution of the form y(x) = ... where the lhs is the unknown of the ODE problem

Description

• 

The IVPsol command specializes the integration constants of the form _Cn (where n is an integer) by entering the given solution sol to match the given initial or boundary conditions ics.

• 

This command is part of the DEtools package, and so it can be used in the form IVPsol(..) after executing the command with(DEtools). However, it can always be accessed through the long form of the command by using DEtools[IVPsol](..).

Examples

withDEtools:

ode1xⅆ2ⅆx2yx+4yx2=0

ode1xⅆ2ⅆx2yx+4yx2=0

(1)

A general solution for ode[1] is not known. Therefore, in principle, it would be not possible to solve an initial value problem that contains it, for example, by using the following initial conditions.

ic1y1=12,Dy1=12

ic1y1=12,Dy1=12

(2)

However, suppose that you find a solution for ode[1] that has been computed by another means, for example

sol1yx=_C12x

sol1yx=_C12x

(3)

You can pass the initial conditions and the solution directly to the DEtools[IVPsol] command.

DEtoolsIVPsolic1,sol1

yx=12x

(4)

For the following example, the general solution of the ODE involves the inverse of an elliptic integral, represented by the elliptic function JacobiSN.

ode2ⅆ2ⅆx2yx2yx3=0

ode2ⅆ2ⅆx2yx2yx3=0

(5)

gensol2dsolveode2

gensol2yx=_C2JacobiSNIx+_C1_C2,I

(6)

Consider the following initial conditions for ode[2]:

ic2y0=1,Dy0=1

ic2y0=1,Dy0=1

(7)

When solving the complete ODE-IVP, instead of directly using the general solution gensol[2], dsolve first searches for symmetries from where a particular solution could be computed. The computed solution could be easier to specialize to satisfy the given conditions y(0) = 1, D(y)(0) = 1, resulting in the following ODE-IVP solution.

dsolveode2,ic2

yx=1x1

(8)

If you use the following command to view the results of trying to specialize the general solution gensol[2], instead of using a short-hand form for the process using symmetries, the following results are returned. This is the solution returned in Maple 12.

DEtoolsIVPsolic2,gensol2

yx=RootOfI_Z2JacobiCNInverseJacobiSN1_Z,I,IJacobiDNInverseJacobiSN1_Z,I,I1JacobiSNIx+InverseJacobiSN1RootOfI_Z2JacobiCNInverseJacobiSN1_Z,I,IJacobiDNInverseJacobiSN1_Z,I,I1,IRootOfI_Z2JacobiCNInverseJacobiSN1_Z,I,IJacobiDNInverseJacobiSN1_Z,I,I1RootOfI_Z2JacobiCNInverseJacobiSN1_Z,I,IJacobiDNInverseJacobiSN1_Z,I,I1,I

(9)

Yet another general solution, involving elliptic integrals themselves, is obtained using Lie symmetry methods to compute it; this is the solution returned in Maple 11 and earlier. It can also be specialized.

gensol2.1dsolveode2,Lie,explicit

gensol2.1yx=RootOf_Z1_f4+_C1ⅆ_f+x+_C2,yx=RootOf_Z1_f4+_C1ⅆ_f+x+_C2

(10)

DEtoolsIVPsolic2,gensol2.11

yx=RootOf_Z11_f4ⅆ_f+x

(11)

As illustrated, the solution returned by dsolve without extra arguments (that is, dsolve([ode[2], ic[2]])) is the simplest one. This solution is obtained by first computing a particular solution from the symmetries (see DEtools,particularsol).

partsol2DEtoolsparticularsolode22

partsol2yx=1x+_C1

(12)

This particular solution, which depends on only one arbitrary constant _C1, can also be specialized.

DEtoolsIVPsolic2,partsol2

yx=1x+1

(13)

See Also

DEtools

dsolve

dsolve, Lie