computes formal power series solutions for differential equation ideals - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Differential Equations : Differential Algebra : DifferentialAlgebra/PowerSeriesSolution

DifferentialAlgebra[PowerSeriesSolution] - computes formal power series solutions for differential equation ideals

Calling Sequence

PowerSeriesSolution(regchain, initialvalues, opts)

PowerSeriesSolution(conditions, regchain, initialvalues, opts)

Parameters

regchain

-

a regular differential chain or a list of them representing an ideal

initialvalues

-

optional, a list or a set of equations specifying initial values for the independent and dependent variables and their derivatives

disregard = ...

-

optional, to disregard whether the initial values cancel the equations and/or inequations of regchain, the right-hand-side can be: nothing (default), inequations, equations, all

order = n

-

optional, a nonnegative integer specifying the order of the series expansion, default value is Order

notation

-

optional, to request the output to be in any of jet, tjet, diff or Diff notations

Description

• 

The function call PowerSeriesSolution(regchain, initialvalues, order = n) returns an integral (exponents are nonnegative integers), formal (the convergence issue is not addressed), power series solution of the differential ideal defined by regchain, truncated at order = n. The initialvalues argument provides both the expansion point of the series and the (generalized) initial values. The result is a list of equations u=series (one for each dependent variable) in the diff or the Diff notation (default is diff).

Conditions

• 

The function call PowerSeriesSolution(conditions, regchain) returns a list containing two lists, one with equations and another one with inequations, that the initial values must satisfy so that an integral formal power series can be computed. The function call PowerSeriesSolution(conditions, regchain, initialvalues) returns the same system (equations and inequations) but evaluated at the expansion points specified in the given initial values.

• 

The list of equations returned is the set of the regular differential chain polynomials which have leading degrees greater than or equal to 2. See Equations. The list of inequations is the set of initials and separants of the regular differential chain polynomials. See Inequations.

• 

The initialvalues argument is a list or a set of equations of the form vi=ai, where the vi are independent variables (specify the expansion point) or dependent variables and/or their derivatives and the ai are general expressions. Independent variables which do not occur among the vi are taken to be 0.

• 

The right-hand-sides of the vi=ai, such that the vi are derivatives, provide the arbitrary constants and/or functions the series depends on. They provide also values for the constants which are algebraic over the other ones. The algorithm proceeds as follows.

– 

First a generic sum is computed, where each derivative is replaced by its normal form.

– 

Second, some substitutions are performed according to the following rules. Consider a derivative u, occurring in some term of the generic sum.

• 

If u is not the derivative of any vi, then u is left unevaluated.

• 

If u is the derivative of exactly one vi, say, u=θvi for some derivation operator θ, then, θai is substituted to vi. The differentiation θai is performed using diff or Diff, depending on the notation (see below in the Options section).

• 

If u is the derivative of many different vi, then, an equation such that vi has maximal order is selected. The substitution is then performed as explained just above.

– 

Last, the resulting expression is evaluated, at the expansion point.

• 

The values attributed to the derivatives at the expansion point, cannot be chosen freely. They must satisfy the system F=0,S0, returned by PowerSeriesSolution when using the option conditions. These conditions, which are checked by default, can be disregarded using the option disregard = ...; see the Examples section.

– 

If the attributed values are solutions of the system F=0,S0, then, the returned power series are truncated solutions of the differential ideal.

– 

If these values do not satisfy F=0, then, the returned power series are not truncated solutions of the differential ideal.

– 

If they satisfy F=0, but do not satisfy S0, then, the returned power series may, or may not, be truncated solutions of the differential ideal. If they are, they do not need to be unique.

• 

This command is part of the DifferentialAlgebra package. It can be called using the form PowerSeriesSolution(...) after executing the command with(DifferentialAlgebra). It can also be directly called using the form DifferentialAlgebra[PowerSeriesSolution](...).

Examples

withDifferentialAlgebra:

A linear ODE

R:=DifferentialRingderivations=t,blocks=u

R:=differential_ring

(1)

The ideal for the ODE for ⅇt consists of a single element, pick it up

exponential:=RosenfeldGroebnerutu,R1

exponential:=regular_differential_chain

(2)

There are no constraints on the initial values for linear systems. Here is the generic series; when the order is not specified, it is set to the value of Order

PowerSeriesSolutionexponential

ut=u0+u0t+12u0t2+16u0t3+124u0t4+1120u0t5+1720u0t6

(3)

The generic series, at t=t0 up to order 3, the initial values can be entered as a list [...] or a set {...}

PowerSeriesSolutionexponential,ut0=c,order=3

ut=c+ctt0+12ctt02+16ctt03

(4)

The initial values can also be expressed directly in jet notation, as in t=t0,u=c, where the expansion point needs to be specified because it is not the origin. When the order is not specified, the value of Order is used; to illustrate we enter the initial values in jet notation

PowerSeriesSolutionexponential,t=t0,u=c

ut=c+ctt0+12ctt02+16ctt03+124ctt04+1120ctt05+1720ctt06

(5)

A nonlinear ODE

The ideal for the `'`2u23u=0 with u0 contains only one case; pick it up

ideal:=RosenfeldGroebnerut24u,u0,R1

ideal:=regular_differential_chain

(6)

There are conditions to be satisfied by the initial values of a formal power series for this problem

PowerSeriesSolutionconditions,ideal

Du024u0=0,Du00

(7)

The initial values u0=c2,ut0=2c satisfy these conditions; to verify that you can enter these initial values as u0=c2,Du0=2c, or recalling that when the expansion point is not specified it is taken equal to 0, you can also enter them in jet notation and without specifying t=0, as in

PowerSeriesSolutionconditions,ideal,u=c2,ut=2c

0=0,c0

(8)

The series, which turns out to be a polynomial, can thus be computed; to illustrate, instead of jet notation we use here the equivalent, standard Maple notation for initial values

PowerSeriesSolutionideal,u0=c2,Du0=2c

ut=c2+2ct+t2

(9)

The conditions computed imply on the inequation u0. Let's verify for instance that initial value u=0 leads to a contradiction

PowerSeriesSolutionconditions,ideal,u=0,ut=0

0=0,00

(10)

One can still compute a formal series solution around t=0 and with u=0, disregarding that the inequations required are not satisfied by using the option disregard = inequations, but the condition ut0, passed to RosenfeldGroebner when computing the ideal, is not satisfied (around t=0) by the solution returned

PowerSeriesSolutionideal,u=0,ut=0,disregard=inequations

ut=t2

(11)

In the example under consideration, if you do not specify initial values, their default values are used, that is: each function and their derivatives are all taken at the origin and are assumed to be all different from zero. These default values do not annihilate all the equations of each element (regular chain) of the ideal, and thus the series cannot be computed

PowerSeriesSolutionideal,disregard=inequations

Error, (in DifferentialAlgebra:-PowerSeriesSolution) the initial conditions do not annihilate all the equations of the ideal

One may still wish to compute a generic formula for the solution, keeping in mind that the initial values must annihilate the equations of the ideal and ignoring for now what these initial values are (say, to be adjusted afterwards). For this purpose use the option disregard = equations (ignore the conditions on the equations)

PowerSeriesSolutionideal,disregard=equations

ut=u0+Du0t+t2

(12)

The heat equation

The ranking is orderly. Thus, the formal power series solutions of this equation, involves all the derivatives of the form uxn and ut,xn where n stands for any nonnegative integer. Classically, one says that the solutions of the heat equation depend on two arbitrary functions u0,x=fx and ut0,x=gx. Such a formal power series can be generated by FormalPowerSeries.

R:=DifferentialRingderivations=t,x,blocks=u

R:=differential_ring

(13)

heat:=RosenfeldGroebnerut,tux,R1

heat:=regular_differential_chain

(14)

PowerSeriesSolutionheat,u=fx,ut=gx,order=2

ut,x=f0+Df0x+g0t+12D2f0x2+Dg0tx+12Df0t2

(15)

A general PDE

The regular differential chains computed by the RosenfeldGroebner function permit to compute integral formal power series solutions of the input system. One considers, here, a polynomial PDE system.

R:=DifferentialRingderivations=x,y,blocks=u,v

R:=differential_ring

(16)

syst:=ux,x2v+ux,xv+ux,ux,y,uy,y21,ux0

syst:=vux,x2+vux,x+ux,ux,y,uy,y21,ux0

(17)

The radical of the differential ideal generated by syst is an intersection of two differential ideals presented by regular differential chains. The solution set of syst is thus the union of the solution sets of the two chains.

ideal:=RosenfeldGroebnersyst,R

ideal:=regular_differential_chain,regular_differential_chain

(18)

Inspecting the left-hand sides of the equations of the first component, one sees that the solutions depend on one arbitrary function vx,0=fx, three arbitrary constants u0,0=c0, ux0,0=c1, uy0,0=c2, and, two constants which are algebraic over the other ones: uy,y0,0 and ux,x0,0.

Equationsideal1,solved

ux,x2=vux,x+uxv,ux,y=0,uy,y2=1,vy=0

(19)

Let us choose some values. Some of them are symbolic, the other ones are numeric. First one verifies that they satisfy the above system. Then, the series is computed.

iv:=v=1+xfx,uy,y=1,u=c0,ux=0,uy=c2,ux,x=1

iv:=v=1+xfx,uy,y=1,u=c0,ux=0,uy=c2,ux,x=1

(20)

PowerSeriesSolutionconditions,ideal1,iv

0=0,0=0,

(21)

PowerSeriesSolutionideal1,iv,order=3

ux,y=c0+c2y+12y212x216x3,vx,y=1+f0x+Df0x2+12D2f0x3

(22)
• 

Inspecting the left-hand sides of the equations of the second component, one sees that the solutions depend on three arbitrary constants u0,0=c0, uy0,0=c1, v0,0=c2, and, one constant uy,y0,0 which is equal to +/-1.

Equationsideal2,solved

uy,y2=1,ux=14v,vx=2,vy=0

(23)

sols:=PowerSeriesSolutionideal2,u=c0,uy=c1,v=c2,uy,y=1,order=10

sols:=ux,y=c0+c1y+14c2x+12y214x2,vx,y=c22x

(24)

These solutions are polynomials. One can plug them into the input equations, and verify the result.

sys:=Tools:-FromJetsyst,u,vx,y

sys:=2x2ux,y2vx,y+2x2ux,yvx,y+xux,y,2yxux,y,2y2ux,y21,xux,y0

(25)

evalsys,sols

0,0,0,14c212x0

(26)

See Also

DifferentialAlgebra, NormalForm, DEtools[initialdata],DEtools[rtaylor]


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