dsolve
find series solutions to ODE problems
Calling Sequence
Parameters
Options
Description
Examples
References
Compatibility
dsolve(ODE, y(x), 'series')
dsolve(ODE, y(x), 'series', x=pt)
dsolve({ODE, ICs}, y(x), 'series')
dsolve({sysODE, ICs}, {funcs}, 'series')
dsolve(ODE, y(x), 'type=series')
dsolve(ODE, y(x), 'type=series', x=pt)
dsolve({ODE, ICs}, y(x), 'type=series')
dsolve({sysODE, ICs}, {funcs}, 'type=series')
ODE
-
ordinary differential equation
y(x)
dependent variable (indeterminate function)
ICs
initial conditions for y(x) and/or its derivatives
sysODE
system of ODEs
{funcs}
set with indeterminate functions
pt
expansion point for series
'type=series'
to request a series solution
combined = true or false
If true, returns a single series. The default is false.
The dsolve command uses several methods when trying to find a series solution to an ODE or a system of ODEs. When initial conditions or an expansion point are given, the series is calculated at the given point; otherwise, the series is calculated at the origin.
The first method used is a Newton iteration based on a paper of Keith Geddes. See the References section in this help page.
The second method involves a direct substitution to generate a system of equations, which may be solvable (by solve) to give a series.
The third method is the method of Frobenius for nth order linear DEs. See the References section in this help page.
If the aforementioned methods fail, the function invokes LinearFunctionalSystems[SeriesSolution].
The series solutions may be expressed as the sum of several series (based on the solution method), rather than as a single series. If a single series is desired, the optional argument combined should be provided to assure that the output is provided in a single series data structure.
ode≔diff⁡y⁡t,t,t+diff⁡y⁡t,t2=0
ode≔ⅆ2ⅆt2y⁡t+ⅆⅆty⁡t2=0
When the initial conditions are not given, the answer is expressed in terms of the indeterminate function and its derivatives evaluated at the origin.
ans≔dsolve⁡ode,y⁡t,type=series
ans≔y⁡t=y⁡0+D⁡y⁡0⁢t−12⁢D⁡y⁡02⁢t2+13⁢D⁡y⁡03⁢t3−14⁢D⁡y⁡04⁢t4+15⁢D⁡y⁡05⁢t5+O⁡t6
If initial conditions are given, the series is calculated at that the given point:
ans≔dsolve⁡ode,y⁡a=Y_a,D⁡y⁡a=DY_a,y⁡t,type=series
ans≔y⁡t=Y_a+DY_a⁢t−a−12⁢DY_a2⁢t−a2+13⁢DY_a3⁢t−a3−14⁢DY_a4⁢t−a4+15⁢DY_a5⁢t−a5+O⁡t−a6
Alternatively, an expansion point can be provided, which is most useful when initial conditions cannot be given:
ans≔dsolve⁡1−t2⁢diff⁡y⁡t,t,t−2⁢t⁢y⁡t−y⁡t,y⁡t,series,t=1
ans≔y⁡t=c__1⁢t−1⁢1−34⁢t−1+748⁢t−12+1128⁢t−13−15715360⁢t−14+3371921600⁢t−15+O⁡t−16+c__2⁢ln⁡t−1⁢−32⁢t−1+98⁢t−12−732⁢t−13−3256⁢t−14+15710240⁢t−15+O⁡t−16+1−2916⁢t−12+2132⁢t−13−1313072⁢t−14−2219102400⁢t−15+O⁡t−16
The above solution is expressed as a sum of several series. The following calling sequence provides the combined form:
ans≔dsolve⁡1−t2⁢diff⁡y⁡t,t,t−2⁢t⁢y⁡t−y⁡t,y⁡t,series,combined,t=1
ans≔y⁡t=c__2+c__1−3⁢c__2⁢ln⁡t−12⁢t−1+−3⁢c__14+c__2⁢9⁢ln⁡t−18−2916⁢t−12+7⁢c__148+c__2⁢−7⁢ln⁡t−132+2132⁢t−13+c__1128+c__2⁢−3⁢ln⁡t−1256−1313072⁢t−14+−157⁢c__115360+c__2⁢157⁢ln⁡t−110240−2219102400⁢t−15+O⁡t−16
The order of the series expansion (default = 6) can be changed using (an environment variable - see Order). For example,
Order≔3
An example with a system of ODEs.
sys≔diff⁡x⁡t,t=y⁡t,diff⁡y⁡t,t=−x⁡t
sys≔ⅆⅆtx⁡t=y⁡t,ⅆⅆty⁡t=−x⁡t
ans≔dsolve⁡sysunionx⁡0=A,y⁡0=B,x⁡t,y⁡t,type=series
ans≔x⁡t=A+B⁢t−12⁢A⁢t2+O⁡t3,y⁡t=B−A⁢t−12⁢B⁢t2+O⁡t3
An example solved by LinearFunctionalSystems[SeriesSolution].
sys≔diff⁡y1⁡x,x−y1⁡x+x⁢y2⁡x=x3,x⁢diff⁡y2⁡x,x−2⁢y2⁡x
sys≔ⅆⅆxy1⁡x−y1⁡x+x⁢y2⁡x=x3,x⁢ⅆⅆxy2⁡x−2⁢y2⁡x
vars≔y1⁡x,y2⁡x
dsolve⁡op⁡sysuniony1⁡0=13,vars,series
y1⁡x=13+13⁢x+132⁢x2+O⁡x3,y2⁡x=D2⁡y2⁡02⁢x2+O⁡x3
Forsyth, A.R. Theory of Differential Equations. Cambridge: University Press, 1906. pp. 78-90
Geddes, Keith. "Convergence Behaviour of the Newton Iteration for First Order Differential Equations". Proceedings of EUROSAM '79. pp.189-199.
Ince, E.L. Ordinary Differential Equations. Dover Publications, 1956. pp. 398-406.
The dsolve command was updated in Maple 2023.
The combined option was introduced in Maple 2023.
For more information on Maple 2023 changes, see Updates in Maple 2023.
See Also
DEtools
dsolve,ICs
dsolve,inttrans
dsolve,Lie
dsolve,numeric
dsolve,system
LinearFunctionalSystems[SeriesSolution]
plots[odeplot]
series
Slode
Download Help Document