DEtools - Maple Programming Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Differential Equations : DEtools : Linear DE Manipulation : DEtools/indicialeq

DEtools

  

indicialeq

  

compute the indicial polynomial of a homogeneous linear ODE

 

Calling Sequence

Parameters

Description

Examples

References

Calling Sequence

indicialeq(des, ivar, alpha, dvar)

Parameters

des

-

homogeneous linear ordinary differential equation or its list form

ivar

-

independent variable

alpha

-

point

dvar

-

(optional) dependent variable

Description

• 

This routine is used to determine either the indicial equation at the point ivar=α for a given homogeneous linear ordinary differential equation, or an equivalent list form (either as generated by DEtools[convertAlg] or as a coefficient list). If the ODE is regular singular (also called Fuchsian) at ivar=α, then the degree of the indicial equation equals the order of the ODE; otherwise it is smaller than the order of the ODE. The roots of the indicial equation are called the exponents of the ODE at ivar=α.

• 

des may be in the standard differential equation form, or in one of two list forms: 1) a list as generated by DEtools[convertAlg], or 2) a coefficient list (as in the first element of a list generated by DEtools[convertAlg]).

• 

dvar must be specified when des is input as a differential equation.

• 

The definition of indicial equation is as follows. Let n be a new variable, substitute dvar=ivarαn into the equation, take the series expansion at ivar=α, and take the coefficient of the first non-zero term. The result will be a polynomial in n, which by definition is the indicial equation (indicialeq replaces n by ivar so that no new variable is needed, and also divides by the leading coefficient). This definition is illustrated by an example below.

• 

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

Examples

withDEtools:

Coefficient list form:

coefs21x2x+1,0,100x2x12:

indicialeqcoefs,x,0

x2x+21100=0

(1)

indicialeqcoefs,x,1

x2x+21100=0

(2)

Differential equation form:

DE7211x173+58+263yx+72xx1720x19+3x4ⅆⅆxyx+x7x4x1ⅆ2ⅆx2yx:

indicialeqDE,x,47,yx

x25931764x=0

(3)

indicialeqDE,x,1,yx

x234x=0

(4)

indicialeqDE,x,0,yx

x249x=0

(5)

This could also have been computed directly from the definition of the indicial equation, as follows:

PDEyx=xn|DEyx=xn:

PseriesnormalPxn,x=0:

Psubsn=x,tcoeffconvertP,polynom,x:

PlcoeffP,x

x249x

(6)

ConvertAlg form:

Y2x2+5x3ⅆ2ⅆx2yx+5xx2ⅆⅆxyx+1+xyx=0:

YconvertAlgY,yx

Y1+x,x2+5x,5x3+2x2,0

(7)

indicialeqY,x,25,yx

x23710x=0

(8)

indicialeqY,x,0,yx

x2+32x+12=0

(9)

References

  

[1] Ince, E.L. Ordinary Differential Equations, Chapters XVI and XVII. New York: Dover Publications, 1956.

See Also

DEtools

DEtools/gen_exp

DEtools[convertAlg]