specializes ODE solutions to match given initial or boundary conditions
a set or list of initial or boundary conditions for an ODE problem. This value can include the ODE itself.
a solution of the form y(x) = ... where the lhs is the unknown of the ODE problem
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](..).
ode1 ≔ x⁢ⅆ2ⅆx2⁢y⁡x+4⁢y⁡x2=0
A general solution for ode 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.
ic1 ≔ y⁡1=−12,D⁡y⁡1=12
However, suppose that you find a solution for ode that has been computed by another means, for example
sol1 ≔ y⁡x=_C12⁢x
You can pass the initial conditions and the solution directly to the DEtools[IVPsol] command.
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ⅆx2⁢y⁡x−2⁢y⁡x3=0
gensol2 ≔ dsolve⁡ode2
Consider the following initial conditions for ode:
ic2 ≔ y⁡0=1,D⁡y⁡0=1
When solving the complete ODE-IVP, instead of directly using the general solution gensol, 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.
If you use the following command to view the results of trying to specialize the general solution gensol, 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.
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.1 ≔ dsolve⁡ode2,Lie,explicit
As illustrated, the solution returned by dsolve without extra arguments (that is, dsolve([ode, ic])) is the simplest one. This solution is obtained by first computing a particular solution from the symmetries (see DEtools,particularsol).
partsol2 ≔ DEtoolsparticularsol⁡ode22
This particular solution, which depends on only one arbitrary constant _C1, can also be specialized.
Download Help Document
What kind of issue would you like to report? (Optional)