arxiv: v1 [math.gm] 16 Feb 2018

Transcription

1 Differential Operator Method of Finding A Particular Solution to An Ordinary Nonhomogeneous Linear Differential Equation with Constant Coefficients Wenfeng Chen arxiv: v math.gm 6 Feb 208 Department of Mathematics and Physics, College of Arts and Sciences, SUNY Polytechnic Institute, Utica, NY 3502, USA We systematically introduce the idea of applying differential operator method to find a particular solution of an ordinary nonhomogeneous linear differential equation with constant coefficients when the nonhomogeneous term is a polynomial function, exponential function, sine function, cosine function or any possible product of these functions. In particular, different from the differential operator method introduced in literature, we propose and highlight utilizing the definition of the inverse of differential operator to determine a particular solution. We suggest that this method should be introduced in textbooks and widely used for determining a particular solution of an ordinary nonhomogeneous linear differential equation with constant coefficients in parallel to the method of undetermined coefficients.

2 2 I. INTRODUCTION When we determine a particular solution of an ordinary nonhomogeneous linear differential equation of constant coefficients with nonhomogeneous terms being a polynomial function, an exponential function, a sine function, a cosine function or any possible products of these functions, usually only the method of undetermined coefficients is introduced in undergraduate textbooks 3. The basic idea of this method is first assuming the general form of a particular solution with the coefficients undetermined, based on the nonhomogeneous term, and then substituting the assumed solution to determine the coefficients. The calculation process of this method is very lengthy and time-consuming. In particular, in the case that the assumed particular solution duplicates a solution to the associated homogeneous equation, one must multiply a certain power function to erase the duplication, and students usually get confused and frustrated in applying this method to find a particular solution of the equation. Nevertheless, differential operator method provide a convenient and effective method of finding a particular solution of an ordinary nonhomogeneous linear differential equation of constant coefficients with the nonhomogeneous terms being a polynomial function, an exponential function, a sine function, a cosine function or any possible products of these functions. The efficiency of this method in determining a particular solution is based on the following facts: exponential function is an eigenfunction of D d/dx; sine and cosine functions are eigenfunctions of second-order differential operator D 2 ; There exists an exponential shift theorem when a polynomial of differential operator acts on a product of exponential function with another continuous function. Especially, if we apply the inverse of differential operator acting on a differentiable function, which is an integration over the function, this method can greatly reduce the complication and arduousness of the calculation. However, the introduction to this method only scatters partially in some online lecture notes of differential equations 4. Hence, here we give a systematic introduction to the method including both rigorous mathematical principles and detailed calculation techniques. In Sect.II, we introduce the definitions of differential operator and its higher order versions as well their actions on differentiable functions. We emphasize the fact that exponential functions are eigenfunctions of differential operator and both sine and cosine functions are eigenfunctions of differential operators of even orders. Further, we define the polynomial of differential operator and list its properties when it acts on differentiable functions. We highlight two remarkable properties including eigenvalue substitution rule and exponential shift rule, which will play important roles in solving differential equations.

3 3 In addition, we introduce the kernel of the polynomial operator and point out it is the solution to the associated linear homogeneous differential equations. In Sect. III, based on the Fundamental Theorem of Calculus, we define the inverses of differential operator and its higher order versions. Then according to the definition of a function of differential operator given by the Taylor series we define the inverse of a polynomial of differential operator and show some of its properties when it acts on differentiable functions. Sect. IV contains the main content of the article. We use a large number of examples to show how the differential operator method works efficiently in determining a particular solution for a nonhomogeneous linear differential equation with constant coefficients when the nonhomogeneous term is a polynomial function, an exponential function, a sine, a cosine function or any possible product of these functions. Sect V is a a brief summary on the application of differential operator method in solving nonhomogeneous linear ordinary differential equations with constant coefficients. II. DIFFERENTIAL OPERATOR AND ACTION ON A DIFFERENTIABLE FUNCTION A. Definition of differential operator A differential operator D acting a differentiable function y fx) on R takes the form D d dx ) The action of D and its higher order versions D n on a at least n-times differentiable function y fx) is just to take the derivatives of the function: For example, let Dfx) dfx) dx, 2) D 2 fx) DDfx) d dfx) d2 fx) dx dx dx 2,. D n fx) DD n fx) d dx d n fx) dx n dn fx) dx n 3) y x 3 +x+3e 2x 5sin3x) 4) then Dy d x 3 +x+3e 2x 5sin3x) 3x 2 ++6e 2x +5cos3x) 5) dx

4 4 B. Action of differential operator on elementary functions and eigenvalue of differential operator The actions of differential operator on elementary functions including exponential function, sine and cosine function, and polynomial functions are listed as follows:. Exponential function De λx d dx eλx λe λx 6) It shows e λx is an eigenfunction of D with eigenvalue λ in the representation space composed of differentiable function on R. Straightforwardly, for the higher order version D n, when n 2 is a positive integer, there exists D n e λx dn dx neλx d ) d d dx dx dx eλx λ n e λx 7) 2. Action of D 2 on sine and cosine functions D 2 sinβx) β 2 sinβx), D 2 cosβx) β 2 cosβx) 8) This shows both sinβx) and cosβx) are eigenfunctions of D 2 with eigenvalue β 2. Further, consider the higher order version D 2n D 2) n, we have D 2n sinβx) D 2) n sinβx) D 2 D 2 D ) 2 sinβx β 2) n sinβx), D 2n cosβx) D 2) n cosβx) D 2 D 2 D ) 2 cosβx β 2) n cosβx) 9) 3. Action of D on power function Dx k d dx xk kx k, D n x k dn dx nxk kk ) k n+)x k n 0)

5 5 where k is a positive integer. Especially, D n x k 0 for k < n ) For example, D 5 x 3 d5 dx 5x3 0 2) C. Polynomial of Differential Operator D and Property Definition : Let P n x) a n x n +a n x n + +a x+a 0 3) be a polynomial of real variable x of degree n, where a 0, a,, a n are real constants. Then P n D) a n D n +a n D n +a D + +a 0 4) is called a polynomial of differential operator D of degree n. In the representation space composed of at least n times differentiable functions fx) on R, ) d d n P n D) P n a n dx dx n +a d n n dx n + +a D +a 0 5) The properties of a polynomial of differential operator D are listed as follows:. Linearity Theorem : Let P n D) be a polynomial of differential operator D of degree n defined in 4). Then there exists P n D)afx)+bgx) ap n D)fx)+bP n D)gx) 6) where a and b are constants, and both fx) and gx) are differentiable functions of at least n-times. 2. Sum rule and product rule Theorem 2: Let P n D) and Q m D) be polynomials of differential operator D of degree n and m, respectively. Then there exist ). Sum rule P n D)+Q m D)fx) Q m D)+P n D)fx) P n D)fx)+Q m D)fx) 7)

6 6 where fx) is a differentiable function of at least maxn,m) times. 2). Product rule P n D)Q m D)gx) Q m D)P n D)gx) P n D)Q m D)gx) Q m D)P n D)gx) 8) where fx) is a differentiable function of at least maxn,m) times. Theorems and 2 can be proved straightforwardly with the explicit representation forms of P n D) and Q m D) on the space of differentiable functions. 3. Eigenvalue substitution rule Theorem 3: Let P n D) be a polynomial of differential operator D of degree n. Then there exist ). P n D)e λx P n λ)e λx e λx P n λ) 9) 2). P n D 2 )sinβx P n β 2 )sinβx sinβxp n β 2 ), P n D 2 )cosβx P n β 2 )cosβx cosβxp n β 2 ) 20) where λ and β are real constants. Proof: ). First we have d P n D)e λx P n dx ) e λx d n a n d n dx n +a n dx n + +a D +a 0 ) a n λ n +a n λ n + +a λ+a 0 e λx P n λ)e λx e λx P n λ) 2) ) e λx 2). First applying Eq.8), we obtain P n D 2 )sinβx a n D 2 ) n +a n D 2 ) n + +a D 2 +a 0 sinβx a n β 2 ) n +a n β 2 ) n + +a β 2 )+a 0 sinβx P n β 2 )sinβx sinβxp n β 2 ) 22) Similarly there is the result for P n D 2 )cosβx.

7 7 4. Exponential shift rule Theorem 4: Let P n D) be a polynomial of differential operator D of degree n. Then there exists D n e λx fx) e λx D +λ) n fx), 23) P n D) e λx fx) e λx P n D +λ)fx) 24) where fx) is a differentiable function of at least n times on R. Proof: ). We use the method of mathematical induction to prove the result 23). First, the result arises for n : D e λx fx) d e λx fx) e λ+ λx d ) fx) e λx D +λ)fx) 25) dx dx Second, assume the result holds for n k: D k e λx fx) e λx D +λ) k fx) 26) Then D k+ e λx fx) D ) D k e λx fx) e λx D +λ) D D +λ) k fx) e λx D +λ) k fx) e λx D +λ) k+ fx) 27) Therefore, we have the general result 23). 2). The result 24) follows straightforwardly from 23): P n D) e λx fx) a n D n +a n D n + +a D +a 0 e λx fx) e λx a n D +λ) n +a n D +λ) n + +a D +λ)+a 0 fx) e λx P n D +λ)fx) 28) D. Kernel of Differential Operator According to the general definition on the kernel of an operator, the kernel of differential operator D k, k,2,, in the space of differentiable real functions is defined as follows:

8 8 Definition 3: kerd k { } fx) D k fx) 0 29) Due to Eq.), kerd k is the set of polynomial functions of degree at most k, kerd k span {,x,,x k } { } c k x k +c k 2 x k 2 c0 + +c x+c 0,c,,c k are constants 30) Further, we can write down the kernel of the polynomial P n D) of differential operator D, which is the solution space of a linear n-th order ordinary homogeneous differential equation with constant coefficients, P n D)yx) a n y n) x)+a n y n ) x)+ +a y x)+a 0 yx) 0 3) The kernel kerp n D) takes different forms depending on the form of P n D). Applying the exponential shift rule 23) and eigenvalue substitution rule, we can easily find the following results: Theorem 5: ). If P n D) is product of n distinct linear factors for simplicity, we take a n for now and later), n P n D) D r )D r 2 ) D r n ) D r i ), i r r 2 r n 32) then kerp n D) span{e rx,e r2x,,e rnx } { n } c i e r ix c,c 2,,c n are constants i 33) 2). If P n D) is product of k distinct repeated linear factors, then k P n D) D r ) n D r 2 ) n2 D r k ) n k D r i ) n i, i r r 2 r k, k n i n 34) i { } kerp n D) span xe rx,x 2 e rx,,x n e rx ; ;xe rkx,x 2 e rkx,,x nk e r kx 35)

9 9 3). If P n D) is product of k distinct quadratic factors, and in this case there must be n 2k, then P n D) D α ) 2 +β 2 D α 2 ) 2 +β2 2 D α k ) 2 +βk 2 k i D α i ) 2 +βi 2, α i α j, or β i β j for i j 36) kerp n D) span{e α x cosβ x),e α x sinβ x); ;e α kx cosβ k x),e α kx sinβ k x)} 37) 4). If P n D) is product of p repeated distinct quadratic factors, then P n D) D α ) 2 +β 2 p i n D α 2 ) 2 +β2 2 n2 D α p ) 2 +βp 2 np D α i ) 2 +βi 2 ni, α i α j, or β i β j for i j 38) kerp n D) span{e αx cosβ x),e αx sinβ x),xe αx cosβ x),xe αx sinβ x),,x n e αx cosβ x),x n e αx sinβ x); e α2x cosβ 2 x),e α2x sinβ 2 x),xe α2x cosβ 2 x),xe α2x sinβ 2 x),,x n2 e α2x cosβ 2 x),x n2 e α2x sinβ 2 x); ;e αpx cosβ p x),e αpx sinβ p x),xe αpx cosβ p x),xe αpx sinβ p x), },x np e αpx cosβ p x),x np e αpx sinβ p x) 39) If P n D) is a mixture of the above four cases, so are the corresponding terms from each of the kernels. For example, let then P 2 D) D 2)D 5) 3 D +3) 2 +4 D 7) ) { kerp 2 D) span e 2x ;e 5x,xe 5x,x 2 e 5x ;e 3x cos2x),e 3x sin2x);e 7x cos4x),e 7x sin4x), } xe 7x cos4x),xe 7x sin4x),x 2 e 7x cos4x),x 2 e 7x sin4x),x 3 e 7x cos4x),x 3 e 7x sin4x) a 0 e 2x + b 0 +b x+b 2 x 2) e 5x +c cos2x+c 2 sin2x)e 3x + d 0 +d x+d 2 x 2 +d 3 x 3) e 7x cos4x+ where a 0, b 0, b, b 2, c, c 2, d 0,, d 3 and e 0,, e 3 are real numbers. e 0 +e x+e 2 x 2 +e 3 x 3) e 7x sin4x 4)

10 0 III. INVERSE OF DIFFERENTIA OPERATOR AND ACTION ON A CONTINUOUS FUNCTION A. Definition of the Inverse of Differential Operator The inverse of differential operator D can be defined according to the Fundamental Theorem of Calculus. Because there exists d x x fx )dx D fx )dx fx) 42) dx x 0 x 0 wherefx) is a continuous function on afiniteinterval a,b andx 0 is anarbitrary constant on a,b, therefore, the action of the inverse D of the differential operator D on a continuous function fx) can be defined in terms of the following definite integral: D fx) x x 0 fx )dx 43) Further, we can successively define the action of higher order inverse differential operator on a continuous function on a finite interval: D 2 fx) D 3 fx) D ) 2 x x fx) D D fx) dx fx 2 )dx 2, x 0 x 0 44) D ) 3 x x x2 fx) dx dx 2 fx 3 )dx 3, x 0 x 0 x 0 45). D n fx) D ) n x fx) dx x 0 x x 0 dx 2 xn x 0 fx n )dx n 46) The justification of these definitions can be demonstrated by the following examples:. D cos3x) x x 0 cos3x )dx 3 sin3x) sin3x 0) 3 sin3x)+c sin3x)+kerd 47) 3 x x x 2. D 2 ) x dx 2 x 2 dx 2 x 2 x 2 0 x 0 x x3 2 x2 0x+ 3 x3 0 x 0 6 x3 +c x+c 2 6 x3 +kerd 2 48) B. Function of Differential Operator As a generalization of the polynomial P n D), a function of the differential operator D can be defined in a way similar to the definition of a square matrix 5.

11 Definition 2: Let fx) be a uniformly convergent function on R. Then a function of fd) of the differential D is defined in terms of the Taylor series expansion of fx) about x 0: fd) n0 n! fn) 0)D n n0 n! d n fx) dx n D n x0 n0 For examples, two typical functions of D are listed as follows: ). n D 2). e D n0 n0 n! d n fx) dx n d n x0 dx n 49) D n +D +D 2 + +D n + 50) n! Dn +D + 2! D2 + + n! Dn + 5) C. Inverse of Polynomial of Differential Operator Using the Taylor series representation 49) on a function of differential operator D, the product rule 8) of polynomials of differential operator D, and the definitions 43), 44), 45), 46) of the action of the inverse operator D k on a continuous function, we can define the inverse of a polynomial operator P n D) as follows: Definition 4: Let P n D)fx) a n D n +a n D n + +a D +a 0 ) fx) gx) 52) where both fx) and gx) are differentiable functions. ). If a 0 0, then fx) P n D) gx) P n D) gx) a n D n +a n D n gx) + +a D +a 0 ) p an D n + + a ) p D gx) 53) a 0 a 0 a 0 p0 2). If all a 0 a a k 0 k n) and a k 0, then fx) P n D) gx) P n D) gx) a n D n +a n D n + +a k D k gx) a n D n +a n D n + +a k D k gx) a n D n k +a n D n k D k gx) 54) + +a k

12 2 The definition of P n D) and the product rule 8) yields the following property Corollary : Let P n D) and Q m D) be two polynomials of D of degree n and m, respectively. Then P n D)Q m D) fx) Q m D)P n D) fx) fx) Q m D) Q m D) P n D) P n D) where fx) is a differentiable function on a certain interval. fx) 55) Corollary 2: Let P n D) and Q m D) be two polynomials of D of degree n and m, respectively. Then P n D) fx) P n D)Q m D) Q md)fx) Q m D) P n D)Q m D) fx) 56) where fx) is a differentiable function on a certain interval. IV. DIFFERENTIAL OPERATOR METHOD OF SOLVING NONHOMOGENEOUS LINEAR ORDINARY DIFFERENTIAL EQUATION WITH CONSTANT COEFFICIENTS The general form of an n-th order nonhomogeneous linear ordinary differential equation with constant coefficients takes the following form: d n y a n dx n +a d n y n dx n + +a dy dx +a 0y gx) 57) where a 0, a,, a n are constants. Expressed in terms of differential operator, the equation reads a n D n +a n D n + +a D +a 0 ) y P n D)y gx) 58) Then formally a particular solution of the equation is Yx) gx) 59) P n D) Due to the exponential shift rule 24) and the eigenvalue substitution rule 9), 20) as well as the feature ), this approach, like the method of undetermined coefficients, works only when gx) is a polynomial function P k x) of degree k, an exponential function e αx, sinβx), cosβx) or any possible product of these functions. In the following we discuss each case of gx)

13 3 A. The case gx) is an exponential function We first give the following results for gx) Ae αx in order to show how the differential operator approach works in this case, where A is a constant. Theorem 6: A particular solution to a differential equation P n D)yx) Ae αx 60) reads as follows:. If P n α) 0, then Yx) P n D) Aeαx ) P n α) Aeαx ) 6) 2. If P n α) 0, i.e., P n D) D α) k P n k D), k n then Yx) P n D) Aeαx ) ) A P n k α) k! xk +kerd k e αx 62) Proof: According to the rule of eigenvalue substitution 9), we have P n D)Ae αx AP n α)e αx 63) ). Since P n α) 0, we divide the equation 63) by P n α), then Ae P n α) P nd)ae αx αx ) ) P n D) Ae αx 64) P n α) In comparison with the differential equation 60), it shows Yx) Ae αx /P n α). 2). Since P n α) 0, the polynomial P n D) must take the following form for a certain k n, P n D) P n k D)D α) k, P n k α) 0, 65) and the equation becomes P n D)yx) D α) k P n k D)yx) Ae αx 66)

14 4 Because D α) k P n k D) k! xk +kerd )e k αx ) P n k α)d α) k e αx k! xk +kerd k P n k α)e αx D k k! xk +kerd k ) P n k α)e αx 67) therefore, yx) given in 62) is the solution to the differential equation. Theorem 6 shows a formal method of finding a particular solution of differential equation when the nonhomogeneous term is an exponential function, Yx) P n D) Aeαx ) Aeαx P n k α) D k Axk e αx k!p n k α) P n k D)D α) k ) Aeαx ) 68) The following two examples shows how this method works. ). A particular solution to the differential equation 3y 2y +6y 5e 3x 69) is Yx) 3D 2 5e 3x) 2D +8 5e 3x e3x 70) 2). A particular solution to the differential equation D )D +5)D 2) 3 yx) 3e 2x 7) is Yx) D )D +5)D 2) 3 3e 2x) 3e2x 7 3e 2x 2 )2+5) D 3 3! x3 4 x3 e 2x 72) A straightforward result comes from Theorem 6, which is called the Exponential Input Theorem. Corollary 3 Exponential Input Theorem):

15 5 Let P n D) be a polynomial of D of degree n. Then the nonhomogeneous linear differential equation 60) has a particular solution of the following form: Ae αx P n α), if P nα) 0 Axe αx P n α), if P nα) 0 but P nα) 0 Ax 2 e αx Yx) P n α), if P nα) P n α) 0 but P n α) 0. Ax k e αx P n k) α), if P nα) P nα) P n k ) α) 0 but P n k) α) 0 73) Proof: We still start from Eq.63). ). For the case P n α) 0, the proof is the same as that for Theorem 6. 2). If P n α) 0, while P nα) 0 we first take the derivative on the the equation 63) with respect to α, then d dα P nd)ae αx ) d dα AP nα)e αx, P n D)Axe αx ) A P n α)eαx +P n α)xe αx ) AP n α)eαx 74) We divide the above equation by P nα) and obtain Axe P n α) P nd)axe αx αx ) P n D) P n α) Ae αx 75) Therefore, a particular solution of the differential equation 60) is yx) Axe αx /P nα). 3). If P n α) P nα) 0, while P nα) 0, we take the derivative the equation 63) twice with respect to α, d 2 dα 2 P nd)ae αx ) d2 dα 2 AP nα)e αx, P n D) Ax 2 e αx) A P nα)e αx +2P nα)xe αx +P nα)x 2 e αx AP nα)e αx 76) Then dividing the above equation by P nα), we have Ax2e αx P n D) P n α) Ae αx 77) Hence in this case the solution to the differential equation is yx) Ax 2 e αx /P n α) 60).

16 6 4). For the general case, P n α) P nα) P n k ) α) 0, while P n k) α) 0, we take the derivative on the equation 63) k times with respect to α, d k dα k P nd)ae αx ) dk dα k AP nα)e αx 78) With an application of the product rule for higher order derivative, it yields P n D) Ax k e αx) A d k k dx k fx)gx) p0 k p k k P n k p) α)x p e αx p0 p A P n k) α)e αx +kp n k ) α)xe αx + AP k) n α)eαx f k p) x)g p) x), 79) kk ) P n k 2) α)x 2 e αx + +P n α)x k e αx 2 Therefore, a particular solution of the equation in this case is yx) Ax k e αx /P k) n α). 80) One can also apply Corollary 3 to find a particular solution of a differential equation. For example, determine a particular solution of the equation Since D 2)D 4) 3 y 5e 4t 8) P D) D 4) 2 4D 0), P D) D 4)2D 36), P 3) D) 2D 7), P4) P 4) P 4) 0, P 3) 4) 0 82) so a particular solution of the equation is Yx) 5x3 e 4x P 3) 4) 5 36 x3 e 4x 83) B. The case gx) is a polynomial function P k x) In this case, due to the fact ) it is most convenient to use Definition 4 see Eqs.53) and 53)) to find a particular solution of a differential equation. We use the following two examples to illustrate the idea.

17 7 Example : Determine a particular solution of the equation y 5y +3y +2y x 2 +3x 2 84) The differential operator form of the equation is ) D 3 5D 2 +3D +2 y 2x 3 +4x 2 6x+5 85) and consequently, a particular solution of the equation is ) Yx) D 3 5D 2 2x 3 +4x 2 6x+5 +3D +2 ) 2+D 3 5D 2 2x 3 +4x 2 6x+5 +3D)/2 2 2 D3 5D 2 +3D)+ 4 5D2 +3D) 2 ) 8 3D)3 2x 3 +4x 2 6x D D2 9 ) 8 D3 2x 3 +4x 6x+5) 2 x x x ) Example 2: Determine a particular solution of the equation y 3y +2y x 3 2x 2 87) A particular solution of the above equation is Yx) D 3 3D 2 x 3 2x 2) x 3 2x 2) +2D D D)2 D) +D +D 2 +D 3) x 3 2x 2) D2 D) + D 2 2 D + 4 D2 + ) 8 D3 x 3 +x 2 +2x+2) 2 D x x x+ 7 ) 4 8 x x x x 88) C. The case gx) is either a sine or cosine function cases: For gx) Asinβx or gx) Acosβx, the solution is formally given by Eq.59) as previous There are two cases: Yx) Asinβx+Bsinβx) 89) P n D)

18 8 ). P n D)sinβx 0 and P n D)cosβx 0 In this case, we first employ the result of Corollary 2 and the simple algebraic operation a+b)a b) a 2 b 2 to make the denominator become a function of D 2, and then apply the eigenvalue substitution rule 20). We use the following example to illustrate this method. Example: Determine a particular solution of the equation 2y +y 5y +3y 3sin2x 90) A particular solution of the equation is Yx) 2D 3 +D 2 5D +3 3sin2x) 2D3 5D) D 2 +3) 2D 3 5D) 2 D 2 +3) 2 3sin2x) 2D 3 D 2 5D 3 D 2 2D 2 5) 2 D 2 +3) 2 3sin2x) 3 4)2 4) d3 4+3) 2 dx 3 d2 dx 2 5 d ) dx 3 sin2x 3 26cos2x sin2x) 9) 677 2). P n D)sinβx 0 and P n D)cosβx 0 In this case, due to following form, D 2 +β 2) sinβx 0 and D 2 +β 2) cosβx 0, P n D) must takes the P n D) P n 2k D) D 2 +β 2) k, k 92) To determine a particular solution in this case, we first introduce the following result. Theorem 7: There exist )p D 2 +β 2 2p cosβx ) 2p)!2β) 2px2p cosβx, )p D 2 +β 2 2p sinβx ) 2p)!2β) 2px2p sinβx, p,2, 93) D 2 +β 2 ) 2p+ cosβx ) p 2p+)!2β) 2p+x2p+ sinβx, D 2 +β 2 ) 2p+ sinβx ) p+ 2p+)!2β) 2p+x2p+ cosβx, p 0,,2, 94)

19 9 Proof: Applying the exponential shift rule, we have D 2 +β 2) k x k e iβx) e iβx D +iβ) 2 +β 2 k x k e iβx D +2iβ) k D k x k e iβx i k 2β) k k! 95) For k 2p, p,2,, using the Euler formula, we obtain D 2 +β 2) 2p ) x 2p cosβx+ix 2p sinβx 2p)!2β) 2p ) p cosβx+isinβx) 96) The real and imaginary parts yield D 2 +β 2) 2p ) x 2p cosβx 2p)!2β) 2p ) p cosβx, D 2 +β 2) 2p ) x 2p sinβx 2p)!2β) 2p ) p sinβx 97) and lead to the result 93). Similarly, for k 2p+, p 0,,, there arises D 2 +β 2) 2p+ ) x 2p+ cosβx+ix 2p+ sinβx 2p+)!2β) 2p ) p icosβx sinβx) 98) The real and imaginary parts give the result 94). We hence can determine a particular solution 89) for the case P n D) 0 as follows: Yx) P n 2k D)D 2 +β 2 ) k Acosβx+Bsinβx) D 2 +β 2 ) k P n 2k D) Acosβx+Bsinβx) 99) For example, we determine a particular solution of the differential equation D ) 2 D 2)D 2 +4) 2 yx) 4sin2x) 00) Then Yx) D ) 2 D 2)D 2 +4) 2 4sin2x) 4 D 2 +4) 2 D ) 2 D 2) sin2x 4 D +) 2 D +2) D 2 +4) 2 D 2 ) 2 D 2 4) sin2x d 3 d2 50 D 2 +4) 2 +4 dx3 dx 2 +5 d ) dx +2 sin2x 50 D 2 +4) 2 2cos2x 4sin2x) 800 x2 cos2x 7sin2x) 0)

20 20 An alternative method for solving the case gx) Acosβx) + Bsinβx is first using the Euler formula to express cos βx and sin βx in terms of complex exponential functions, cosβx 2 e iβx +e iβx), sinβx e iβx e iβx) 02) 2 and then determining a particular solution in the same way as gx) being a complex exponential function. For example, the particular solutions of the equations 90) and 00) can be worked out as follows: ). For the equation 90), we have Yx) 2D 3 +D 2 5D +3 3sin2x) 3 e 2ix 2D 3 +D 2 e 2ix) 5D +3 2i 3 ) e 2ix 2i 2D +2i) 3 +D +2i) 2 5D +2i)+3 ) e 2ix 2D 2i) 3 +D 2i) 2 5D 2i)+3 3 ) e 2ix 2i 26i e 2ix +26i 3 2i i)cos2x+isin2x) 26i)cos2x isin2x) 3 26cos2x sin2x) 03) 677 2). For the equation 00), we have Yx) D ) 2 D 2)D 2 +4) 2 4sin2x) 4 e 2ix D ) 2 D 2)D 2 +4) 2 e 2ix) 2i 2 ) e 2ix i D +2i ) 2 D +2i 2)D +2i) ) e 2ix D 2i ) 2 D 2i 2)D 2i) e 2ix ) i 2i ) 2 2i 2) D 2 D +4i) 2 e 2ix ) 2i ) 2 2i 2) D 2 D 4i) 2 e 2ix 6i 2i) 2 i) + e 2ix +2i) 2 D 2 +i) 800i 3+4i)+i)cos2x+isin2x)+ 3 4i) i)cos2x isin2x) 2 x2 800 x2 cos2x 7sin2x) 04)

21 2 D. The case gx) is a product of exponential function and polynomial function For gx) P m x)e αx, one should first apply the exponential shift rule to move out of the exponential function, and then proceed to determine a particular solution as the case gx) being a polynomial function: For examples, Yx) P n D) P mx)e αx ) e αx P n D +α) P mx) 05) ). A particular solution of the equation D 3) 2 D 2 2D +5)D +2)yx) x 2 3x+)e 2x 06) is Yx) x 2 D 3) 2 D ) 2 3x+)e 2x +4D +2) e 2x D ) 2 D +) 2 +4D +4) x2 3x+) e 2x 20 27D +2D 3 +4D 4 +D 5x2 3x+) e 2x 20 27D x2 3x+) e 2x + 27 ) d dx + d 2 20 dx 2 x 2 3x+) x 2 3 ) 9 x+ e 2x 07) ). A particular solution of the equation D 3)D 2) 2 D +)yx) 4x 2)e 2x 08) is Yx) 4x 2)e 2x D 2) 2 e 2x D 3)D +) D )D +3)D 24x 2) e 2x D 2 D 2 +2D 3) 4x 2) e2x D 2 2D 3 4x 2) e 2x D d 3dx ) 4x 2) 3 e2x D 2 4x+ 2 ) 3 9 x2 2x+)e 2x 09)

22 22 E. The case gx) is a product of polynomial function and sine or cosine functions In this case, gx) P m x)acosβx + Bsinβx). The most convenient way of determining a particular solution is to first express sine and cosine function in terms of complex exponential functions, and then use the method of dealing with the product of a polynomial function and exponential function to find a particular solution: For examples, Yx) P n D) P mx)acosβx+bsinβx) 2 A ib)eiβx P n D +iβ) P mx) + 2 A+iB)e iβx P n D iβ) P mx) 0) ). A particular solution to the differential equation y 4y x 2 3)sin2x ) is ) Yx) D 2 x 2 ) 3 sin2x 4 D 2 x 2 3 e 2ix e 2ix) ) x i e 2ix 2i D +2i) 2 4 e 2ix D 2i) 2 x 3) 2 4 e 2ix 2i D 2 +4iD 8 e 2ix D 2 x 3) 2 4iD 8 i e + 2ix 6 2 id ) e 8 D2 2ix 2 id ) 8 D2 x 2 3) ) 4x 2 3 sin2x+4xcos2x 2) 32 2). A particular solution to the equation y +4y 4x 2 cos2x 3) is Yx) D i 24 ) 4x 2 cos2x e 2ix 2x 2 D 2 e 2ix +e 2ix) +4 D +2i) 2 +4 x2 +e 2ix D 2i) 2 +4 x2 e 2ix DD +4i) x2 +e 2ix DD 4i) x2 e 2ix D 2x 2 +ix ) e 2ix D 2x 2 ix 4 4 6x 2 cos2x+x8x 2 3)sin2x ) 4)

23 23 F. The case gx) is a product of exponential function and sine or cosine functions For gx) e αx Acosβx+Bsinβx), as other cases involving sine or cosine functions, the best way of determining a particular solution for an equation is to first express sine and cosine function in terms of complex exponential functions. There are two possibilities. ). If P n α ± iβ) 0, then one can directly apply the eigenvalue substitution rule to get a particular solution: Yx) P n D) eαx Acosβx+Bsinβx) A ib 2P n α+iβ) eα+iβ)x + A+iB 2P n α iβ) eα iβ)x 5) For example, a particular solution of the equation y 2y +2y e 2x 2cosx 6sinx) 6) is Yx) D 2 e 2x 2cosx 6sinx) 2D +2 +3i)e 2+i)x D 2 + 3i)e 2 i)x 2D +2 e 2+i)x e 2 i)x +3i) 2+i) 2 22+i)+2 + 3i) 2 i) 2 22 i) e2x 7cosx sinx) 7) An alternative way is to first utilize the exponential rule to move out e αx and then use the method of dealing with sine or cosine function to determine a particular solution. For example, a particular solution of the above equation in this approach is obtained as follows: Yx) e 2x D +2) 2 2D +2)+2 2cosx 6sinx) e2x e 2x D D D 2 +2) 2 4D 2 2cosx 6sinx) 5 e2x D 2 +2D +2 2cosx 6sinx) ) d 2 dx 2 2 d dx +2 2cosx 6sinx) 2 5 e2x 7cosx sinx) 8)

24 24 2). If P n α ±iβ) 0, one should first apply the exponential shift rule and then the definition for the inverse of differential operator. For example, we determine a particular solution of the equation D )D 2 2D +5)y 4e 3x cos2x 9) Then Yx) D )D 2 6D +3) D 3) 2 +4D ) D 3) 2 +4 { i e3+2i)x ) 4e 3x cos2x e 3+2i)x +e 3 2i)x} + D 3) i e3 2i)x 2 i)e3+2i)x DD +4i) + 2 +i)e3 2i)x DD 4i) 8i xe3x i)cos2x+isin2x) 8i xe3x +i)cos2x+isin2x) 4 xe3x sin2x cos2x) 20) One can also work first take out e αx with the exponential shift rule and then apply the results 93) and 94) of Theorem 7. The particular solution of the above example in this approach reads Yx) D )D 3) 2 +4) 4e 3x D 2 +4 ) 4e 3x cos2x 4e 3x D +2)D 2 +4) cos2x D 2 D 2 4 cos2x ) e 3x D 2 +4 sin2x+cos2x) 4 xe3x sin2x cos2x) 2) G. The case gx) is a product of polynomial function, exponential function and sine or cosine functions Finally, we consider the most general case, gx) P m x)e αx Acosβx+Bsinβx). The method of getting a particular is the same as the previous case that gx) is a product of a polynomial function and an exponential function: For examples, Yx) P n D) P mx)e αx Acosβx+Bsinβx) 2 A ib) P n D +α+iβ) P mx)+ 2 A+iB) P n D +α iβ) P mx) 22)

25 25 ). A particular solution to the differential equation y 5y +6y e x cos2x+e 2x 3x+4)sinx 23) is Yx) D 2 5D +6 Y x) D 2 5D +6 ex cos2x) e x cos2x+e 2x 3x+4)sinx Y x)+y 2 x), 24) D 2)D 3) ex cos2x) e x D +)D +2) cos2x ex D )D 2) D 2 )D 2 4) cos2x ) ) d d 40 ex dx + dx +2 cos2x 20 ex cos2x+3sin2x), 25) Y 2 x) D 2 3x+4)e 2x cosx 5D +6 3x+4) e 2+i)x +e 2 i)x) 2 D 2)D 3) 2 e2+i)x D +i)d +i) 3x+4) + 2 e2 i)x D i)d i) 3x+4) 2i e2+i)x D +i 3x+4+3i) 2i e2 i)x D i 3x+4 3i) 2 e2x 3x+0)cosx+3x+)sinx 26) 2). A particular solution of the equation y +2y +2y e x ) 3+2sinx+4x 2 cosx 27) is Yx) D 2 e x ) 3+2sinx+4x 2 cosx Y x)+y 2 x)+y 3 x) 28) +2D +2 where Y x) 3e x ) D 2 +2D +2 Y 2 x) 2e x D 2 sinx ) +2D +2 2e x Y 3 x) D 2 +2D +2 3e x ) e x, D +) 2 + 2e x sinx ) D 2 + sinx xe x cosx, 29) ) 4x 2 e x ) cosx D +) 2 4x 2 e x cosx +

26 26 2x 2 D +) 2 e +i)x +e i)x) + 2e +i)x DD +2i) x2 +2e i)x ie +i)x D x 2 +ix 2 3 e x 3x 2 cosx+x2x 2 3)sinx DD 2i) x2 )+ie +i)x D x 2 ix 2 ) 30) V. SUMMARY We have introduced both the mathematical principle and computing technique of using differential operator method to find a particular solution of an ordinary nonhomogeneous linear differential equation with constant coefficients when the nonhomogeneous term is a polynomial function, exponential function, sine function, cosine function or any possible product of these functions. We have reviewed some rules in the differential operations including the eigenvalue substitution rule and the exponential shift rule. Furthermore, we have used a number of examples to illustrate the practical application of this approach. In particular, comparing with the introductions on this differential operator method scattered in some lecture notes, we propose and highlight the application of the inverse of differential operator in determining a particular solution, and overcome the difficulty when a polynomial of differential operator in a differential equation is singular. The differential operator method has great advantages over the well-known method of undetermined coefficients introduced in textbooks in determining a particular solution. Therefore, it is meaningful to make a systematic introduction and put an emphasis on this method. We hope that this method can be systematically introduced in textbooks and widely used for determining a particular solution of an ordinary nonhomogeneous linear differential equation with constant coefficients, in parallel to the method of undetermined coefficients. ACKNOWLEDGMENTS I would like to thank my Colleagues in the Group of Applied Mathematics at SUNY Polytechnic Institute for useful remarks. I am especially indebted to Dr. Andrea Dziubek and Dr. Edmond

27 27 Rusjan for various discussions and encouragements. D.G. Zill, A First Course in Differential Equations with Modelling Applications,0th edition, Brooks/Cole, W.E. Boyce and R.C. DiPrima, Elementary Differential Equations and Boundary Value Problems, 0th edition, Wiley, R. Nagle, E.B. Saff and A.D. Snider, Fundamentals of Differential Equations, 8th edition, Addison- Wesley, See for example, math.mit.edu/ jorloff/suppnotes/suppnotes03/o.pdf ; esg/notes/topic6.pdf 5 See for example, D. Poole, Linear Algeba, 3rd edition, Brooks/Cole, 20.