DEtools - Maple Programming Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Differential Equations : DEtools : Solving Methods : DEtools/matrixDE

DEtools

  

matrixDE

  

find solutions of a linear system of ODEs in matrix form

 

Calling Sequence

Parameters

Description

Examples

Calling Sequence

matrixDE(A, B, t, method=matrixexp)

matrixDE(A, B, t, solution=solntype)

Parameters

A, B

-

coefficients of a system X' t=AtXt+Bt ; if B not specified, then assumed to be a zero vector

t

-

independent variable of the system

method=matrixexp

-

(optional) matrix exponentials

solution=solntype

-

(optional) where solution=polynomial or solution=rational

Description

• 

The matrixDE command solves a system of linear ODEs of the form X't=AtXt+Bt. If B is not specified then it is assumed to be the zero vector.

• 

An option of the form method = matrixexp can be specified to use matrix exponentials (in the case of constant coefficients).

• 

An option of the form solution = polynomial or solution = rational can be specified to search for polynomial or rational solution. In this case, the function invokes LinearFunctionalSystems[PolynomialSolution] or LinearFunctionalSystems[RationalSolution].

  

The command returns a pair St,Pt with St, which is an n by n Matrix, and Pt, which is an n by 1 Vector. A particular solution of the system can be then written in the form Ft=StC0+Pt where C0 is n by 1 and F0=C0+P0. If B is zero then P will also be zero.

• 

If a system is expressed in terms of equations, dsolve can be used instead.

Examples

withDEtools:

Nonconstant homogeneous system

AMatrix2,2,1,t2,t,1

A1t2t1

(1)

solmatrixDEA,t

solⅇtt32BesselI35,2t525ⅇtt32BesselK35,2t525ⅇtBesselI25,2t525tⅇtBesselK25,2t525t,00

(2)

Matrix of arbitrary coefficients

CMatrix2,1:

Verification of solution

Fevalmsol1 &* C+sol2:rhevalmA &* F:

simplifynormaltF1,1rh1,1,symbolic

0

(3)

simplifynormaltF2,1rh2,1,symbolic

0

(4)

Nonhomogeneous system of two variables with constant coefficients

AMatrix2,2,1,1,0,1;BMatrix2,1,tk,tl

A1101

Btktl

(5)

solmatrixDEA,B,t

solⅇtⅇtt0ⅇt,ⅇt2tl2WhittakerMl2,l2+12,tkl+tl2+1WhittakerMl2,l2+12,tkⅇt2+ⅇt2tk2WhittakerMk2,k2+12,tlⅇt2tl2WhittakerMl2,l2+12,tkⅇt2tl2WhittakerMl2,l2+12,tl+tl+1kl+tl2+1WhittakerMl2,l2+12,tⅇt2+ⅇt2tk2WhittakerMk2,k2+12,tⅇt2tl2WhittakerMl2,l2+12,t+tl+1k+tl+1l+tl+1k+1l+1ⅇt2tl2WhittakerMl2+1,l2+12,tl2+tl2+1ⅇt2WhittakerMl2+1,l2+12,tl+ⅇt2tk2WhittakerMk2+1,k2+12,tl2ⅇt2tl2WhittakerMl2+1,l2+12,tl+tl+1l2tl+1lt+tl2+1WhittakerMl2,l2+12,tⅇt2+tl2+1ⅇt2WhittakerMl2+1,l2+12,t+ⅇt2tk2WhittakerMk2+1,k2+12,tⅇt2tl2WhittakerMl2+1,l2+12,ttklt+2tl+1ltl+1ttkt+tl+1l+1t

(6)

Verification of solution

Fevalmsol1 &* C+sol2:rhevalmA &* F+B:

simplifynormaltF1,1rh1,1,symbolic

0

(7)

simplifynormaltF2,1rh2,1,symbolic

0

(8)

Nonconstant homogeneous system with unknown coefficients

AMatrix2,2,1,0,1,ft

A101ft

(9)

solmatrixDEA,t

solftDESolⅆ2ⅆt2_Ytⅆⅆt_Ytftⅆⅆt_Yt_Ytⅆⅆtft+ft_Yt,_Yt+ⅆⅆtDESolⅆ2ⅆt2_Ytⅆⅆt_Ytftⅆⅆt_Yt_Ytⅆⅆtft+ft_Yt,_YtDESolⅆ2ⅆt2_Ytⅆⅆt_Ytftⅆⅆt_Yt_Ytⅆⅆtft+ft_Yt,_Yt,00

(10)

General nonhomogeneous system of two variables with constant coefficients

AMatrix2,2,a,b,c,d

Aabcd

(11)

BMatrix2,1,ft,gt

Bftgt

(12)

solmatrixDEA,B,t

solⅇad+a22ad+4bc+d2t2ⅇa+d+a22ad+4bc+d2t2ⅇad+a22ad+4bc+d2t2d+a+a22ad+4bc+d22bⅇa+d+a22ad+4bc+d2t2da+a22ad+4bc+d22b,ftd+ⅆⅆtft+bgtⅇa+d+a22ad+4bc+d2t2ⅆtⅇta22ad+4bc+d2ftd+ⅆⅆtft+bgtⅇad+a22ad+4bc+d2t2ⅆtⅇad+a22ad+4bc+d2t2a22ad+4bc+d22ftⅇta22ad+4bc+d2ⅇad+a22ad+4bc+d2t2ⅇa+d+a22ad+4bc+d2t2d+2ⅇta22ad+4bc+d2ⅇad+a22ad+4bc+d2t2gtⅇa+d+a22ad+4bc+d2t2b+2ⅆⅆtftⅇta22ad+4bc+d2ⅇad+a22ad+4bc+d2t2ⅇa+d+a22ad+4bc+d2t2+2ftⅇad+a22ad+4bc+d2t2ⅇad+a22ad+4bc+d2t2d+ⅇta22ad+4bc+d2ftd+ⅆⅆtft+bgtⅇa+d+a22ad+4bc+d2t2ⅆtⅇad+a22ad+4bc+d2t2a22ad+4bc+d2ⅇta22ad+4bc+d2ftd+ⅆⅆtft+bgtⅇa+d+a22ad+4bc+d2t2ⅆtⅇad+a22ad+4bc+d2t2a+ⅇta22ad+4bc+d2ftd+ⅆⅆtft+bgtⅇa+d+a22ad+4bc+d2t2ⅆtⅇad+a22ad+4bc+d2t2d2ⅇad+a22ad+4bc+d2t2gtⅇad+a22ad+4bc+d2t2b2ⅆⅆtftⅇad+a22ad+4bc+d2t2ⅇad+a22ad+4bc+d2t2+ftd+ⅆⅆtft+bgtⅇad+a22ad+4bc+d2t2ⅆtⅇad+a22ad+4bc+d2t2a22ad+4bc+d2+ftd+ⅆⅆtft+bgtⅇad+a22ad+4bc+d2t2ⅆtⅇad+a22ad+4bc+d2t2aftd+ⅆⅆtft+bgtⅇad+a22ad+4bc+d2t2ⅆtⅇad+a22ad+4bc+d2t2d2fta22ad+4bc+d22ba22ad+4bc+d2

(13)

solmatrixDEA,t,method=matrixexp

sola22ad+4bc+d2ⅇa+d+a22ad+4bc+d2t2+aⅇa+d+a22ad+4bc+d2t2dⅇa+d+a22ad+4bc+d2t2+ⅇad+a22ad+4bc+d2t2a22ad+4bc+d2aⅇad+a22ad+4bc+d2t2+dⅇad+a22ad+4bc+d2t22a22ad+4bc+d2bⅇa+d+a22ad+