DOCSLIB.ORG

Variation of Parameters Method for Nonlinear Diffusion Equations

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Variation of Parameters Method for Nonlinear Diffusion Equations International Journal of Modern Applied Physics, 2013, 3(1): 48-56 International Journal of Modern Applied Physics ISSN: 2168-1139 Journal homepage:www.ModernScientificPress.com/Journals/ijmep.aspx Florida, USA Article Variation of Parameters Method for Nonlinear Diffusion Equations Rahmatullah and Syed Tauseef Mohyud-Din* Department of Mathematics, HITEC University, Taxila Cantt, Pakistan *Author to whom correspondence should be addressed; E-Mail: [email protected] Article history: Received 30 July 2013, Received in revised form 30 August 2003, Accepted 10 December 2013, Published 27 December 2013. Abstract: In this paper, we apply Variation of Parameters Method (VPM) to solve nonlinear diffusion equations which arise very frequently in various models related to different branches of physics. The suggested algorithm proved to be very efficient and finds the solution without any discretization, linearization, perturbation or restrictive assumptions. Numerical results reveal the complete reliability of the proposed VPM. Keywords Variation of Parameters Method; nonlinear problems; diffusion equation; exact solution. PACS: 02.30 Jr, 02.00.00. 1. Introduction Diffusion [1-3, 10-15] is an important phenomenon in physical models and one of the important equations translating such models is the following nonlinear partial differential equation: 푢푡=(퐷(푢)푢푥)푥, (1) subject to the initial condition 푢(푥, 0) = 푓(푥), (2) Copyright © 2013 by Modern Scientific Press Company, Florida, USA Int. J. Modern App. Physics. 2013, 3(1): 48-56 49 where D(u) denotes the diffusion term which plays a key role in many applications of diffusion processes and appears in several functional forms including power laws and exponential forms. Moreover, slow diffusion process is represented [1-3, 10-15]: 퐷(푢) = 푢푛 푛 > 0, (3) and fast diffusion processes is described as: 퐷(푢) = 푢푛, 푛 < 0. (4) Due to its physical importance, several techniques including decomposition [15, 16], Variational Iteration [11], Homotopy Perturbation [11, 12] have been applied by various researchers to find appropriate solutions. Inspired and motivated by the ongoing research in this area, we apply Variation of Parameters Method (VPM) [4-10] to find solutions of such equations. It is worth mentioning that Ma et. al. [4-10] used Variation of Parameters (VPM) for solving involved non-homogeneous partial differential equations and obtained solution formulas helpful in constructing the existing solutions coupled with a number of other new solutions including rational solutions, solitons, positions, negatons, breathers, complextions and interaction solutions of the KdV equations. It is observed that proposed scheme is highly efficient, very reliable and accurate. Moreover, suggested algorithm is also compatible with the physical nature of such problems and hence can be extended to find appropriate solutions of wide range of problems with diversified physical nature. 2. Variation of Parameters Method (VPM) [4-10] Consider the following second-order partial differential equation ytt ftxyzyyyy(,,,,,,, x y z xx , y yy , y zz ), (5) where t such that t is time, and f is linear or non-linear function of y,,,,,, yx y y y z y xx y yy y zz . The homogeneous solution of (1) is y t,,,, x y z A Bt where A and B are functions of x,, y z and t . Using Variation of parameters method we have following system of equations AB 0, tt Copyright © 2013 by Modern Scientific Press Company, Florida, USA Int. J. Modern App. Physics. 2013, 3(1): 48-56 50 B f , t and hence tt Axyzt ,,,,,,,,,,,, Dxyz sfds Bxyzt Cxyz fds 00 therefore, t yxyzt,,, yxyz ,,,0 tyxyz ,,,0 tsfsxyzyyyy (,,,,, ,,,), y yds t x y z xx yy zz 0 which can be solved iteratively as [4-10] t yk1 xyzt,,, yxyz ,,,0 tyxyz ,,,0 tsfsxyzy (,,,,,,, k y k y k y k , y k , y k ), ds t x y z xx yy zz 0 k 0,1,2, . 3. Numerical Examples In this section, we apply Variation of Parameters Method (VPM) to solve partial differential equations describing various types of diffusion process. Example 3. 1 Consider the fast diffusion process as follows −1 푢푡=(푢 푢푥)푥, (6) subject to the initial condition 2푐 푢(푥, 0) = , (7) (푎+푥)2 where ‘a’ and c≠0 are arbitrary constants. Applying Variation of Parameters Method (VPM), we get 푡 푢 = 퐴 (푥, 푡) + (푢 −1 푢 ) 푑푠. (8) 푛+1 1 ∫0 푛 푛푥 푥 2푐 Here 퐴 (푥, 푡) = 푢 (푥, 푡) = . 1 0 (푎+푥)2 Consequently, following approximants are obtained 2푐 2푡 푢 (푥, 푡) = + , (9) 1 (푎+푥)2 (푎+푥)2 2푐 2푡 푢 (푥, 푡) = + , 2 (푎+푥)2 (푎+푥)2 . Copyright © 2013 by Modern Scientific Press Company, Florida, USA Int. J. Modern App. Physics. 2013, 3(1): 48-56 51 . .. The closed form solution is given as 2(푐+푡) 푢(푥, 푡) = . (10) (푎+푥)2 Fig 1. Graphic depiction of solution of (5, 6) Example 3.2 Consider a slow diffusion process 푢푡 = (푢푢푥)푥 , (11) subject to the initial condition 1 2 푢(푥, 0) = 푥 , 푥 > 0, (12) 푐 where 푐 > 0 is an arbitrary constant. Applying Variation of Parameters Method (VPM), we get 푡 푢 = 퐴 (푥, 푡) + (푢 푢 ) 푑푠, (13) 푛+1 1 ∫0 푛 푛푥 푥 1 2 Here 퐴 (푥, 푡) = 푢 (푥, 푡) = 푥 1 0 푐 푡 푢 (푥, 푡) = 푢 + (푢 푢 ) 푑푠 1 0 ∫0 0 0푥 푥 2 1 6푡 푢 (푥, 푡) = 푥 ( + ) , (14) 1 푐 푐2 Consequently, following approximants are obtained 2 3 2 1 6푡 36푡 72푡 푢 (푥, 푡) = 푥 ( + + + ) , (15) 2 푐 푐2 푐3 푐4 2 3 4 5 2 1 6푡 36푡 216푡 864푡 2592푡 푢 (푥, 푡) = 푥 ( + + + + + + ⋯ ), (16) 3 푐 푐2 푐3 푐4 푐5 푐6 and so on. In the same way rest of the components can be obtained. Therefore solution of equation (1) in the closed form is Copyright © 2013 by Modern Scientific Press Company, Florida, USA Int. J. Modern App. Physics. 2013, 3(1): 48-56 52 푥2 푢(푥, 푡) = . (17) 푐−6푡 Fig. 2: Graphic depiction of solution of (11, 12) Example 3.3 Consider another slow diffusion process in the form of 2 푢푡=(푢 푢푥)푥 , (18) subject the Initial condition 푥+ℎ 푢(푥, 0) = . (19) 2√푐 Applying Variation of Parameters Method (VPM), we get 푡 2 푢 = 퐴 (푥, 푡) + (푢 푢 ) 푑푠, (20) 푛+1 1 ∫0 푛 푛푥 푥 푥+ℎ Here 퐴 (푥, 푡) = 푢 (푥, 푡) = . 1 0 2√푐 Consequently, following approximants are obtained 푡 2 푢 (푥, 푡) = 푢 + (푢 푢 ) 푑푠 (21) 1 0 ∫0 0 0푥 푥 푥+ℎ 푡 푢 (푥, 푡) = (1 + ) , (22) 1 2√푐 2푐 푥+ℎ 푡 3푡2 푡3 푡4 푢 (푥, 푡) = (1 + + + + ) , (23) 2 2√푐 2푐 8푐2 8푐3 64푐4 푥+ℎ 푡 3푡2 5푡3 13푡4 9푡5 푢 (푥. 푡) = (1 + + + + + + ⋯ ), (24) 3 2√푐 2푐 8푐2 16푐3 64푐4 18푐5 Copyright © 2013 by Modern Scientific Press Company, Florida, USA Int. J. Modern App. Physics. 2013, 3(1): 48-56 53 and so on. In the same way rest of the components can be obtained. Therefore solution of equation (18) in the closed form is 1 푥+ℎ 푢(푥, 푡) = . (25) 2 √푐−푡 Fig. 3: Graphic depiction of solution of (18,19) Example 3.4 Consider the next example where 1 푢 = ( 푢 ) , (26) 푡 1+푢2 푥 푥 subject to the initial condition 푢(푥, 0) = tan 푥. (27) Applying Variation of Parameters Method (VPM), we get 푡 1 푢푛+1 = 퐴1(푥, 푡) + ∫ ( 2 푢푛푥)푥 푑푠, (28) 0 1+푢푛 Here 퐴1(푥, 푡) = 푢0(푥, 푡) = tan 푥, 푡 1 푢1(푥, 푡) = 푢0 + ∫ ( 2 푢0푥)푥 푑푠 0 1+푢0 푢1(푥, 푡) = tan 푥, (29) In the same way, we find Copyright © 2013 by Modern Scientific Press Company, Florida, USA Int. J. Modern App. Physics. 2013, 3(1): 48-56 54 푢2(푥, 푡) = tan 푥, (30) and so on. Therefore solution of equation (26) in the closed form is 푢(푥, 푡) = tan 푥. (31) Fig 4: Graphic depiction of solution of (26, 27) Example 3.5 Finally, we consider another problem of the following form 1 푢 = ( 푢 ) , (32) 푡 푢2−1 푥 푥 subject to the initial condition 푢(푥, 0) = − coth 푥. (33) Proceeding as before, we get 푡 1 푢푛+1 = 퐴1(푥, 푡) + ∫ ( 2 푢푛푥)푥 푑푠, (34) 0 푢푛 −1 Here 퐴1(푥, 푡) = 푢0(푥, 푡) = − coth 푥 푡 1 푢1(푥, 푡) = 푢0 + ∫ ( 2 푢0푥)푥 푑푠, 0 푢0 −1 푢1(푥, 푡) = − coth 푥 , (35) 푢2(푥, 푡) = − coth 푥, (36) and so on. Hence the closed form solution is given as Copyright © 2013 by Modern Scientific Press Company, Florida, USA Int. J. Modern App. Physics. 2013, 3(1): 48-56 55 푢(푥, 푡) = − coth 푥. (37) Graph of above solution is shown in figure 5. Fig 5: Graphic depiction of solution of (32, 33) 4. Conclusion Variation of Parameters Method (VPM) is applied to solve nonlinear diffusion equations. The suggested algorithm which is applied without any discretization, linearization, perturbation or restrictive assumptions is very effective, reliable and accurate. Moreover, proposed scheme (VPM) can be implemented on other physical problems also. References [1] Q. Changzheng, Exact solutions to nonlinear diffusion equations obtained by a generalized conditional symmetry method, IMA J. Appl. Math. 62 (1999): 283-302. [2] L. Dresner, Similarity Solutions of nonlinear partial differential equations, Pitman, New York, 1983. [3] J. H. He, An elementary introduction of recently developed asymptotic methods and nanomechanics in textile engineering, Int. J. Mod. Phys. B, 22 (21) (2008): 3487-4578. [4] W. X.
Load more

Recommended publications
  • Modifications of the Method of Variation of Parameters
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector International Joumal Available online at www.sciencedirect.com computers & oc,-.c~ ~°..-cT- mathematics with applications Computers and Mathematics with Applications 51 (2006) 451-466 www.elsevier.com/locate/camwa Modifications of the Method of Variation of Parameters ROBERTO BARRIO* GME, Depto. Matem~tica Aplicada Facultad de Ciencias Universidad de Zaragoza E-50009 Zaragoza, Spain rbarrio@unizar, es SERGIO SERRANO GME, Depto. MatemAtica Aplicada Centro Politdcnico Superior Universidad de Zaragoza Maria de Luna 3, E-50015 Zaragoza, Spain sserrano©unizar, es (Received February PO05; revised and accepted October 2005) Abstract--ln spite of being a classical method for solving differential equations, the method of variation of parameters continues having a great interest in theoretical and practical applications, as in astrodynamics. In this paper we analyse this method providing some modifications and generalised theoretical results. Finally, we present an application to the determination of the ephemeris of an artificial satellite, showing the benefits of the method of variation of parameters for this kind of problems. ~) 2006 Elsevier Ltd. All rights reserved. Keywords--Variation of parameters, Lagrange equations, Gauss equations, Satellite orbits, Dif- ferential equations. 1. INTRODUCTION The method of variation of parameters was developed by Leonhard Euler in the middle of XVIII century to describe the mutual perturbations of Jupiter and Saturn. However, his results were not absolutely correct because he did not consider the orbital elements varying simultaneously. In 1766, Lagrange improved the method developed by Euler, although he kept considering some of the orbital elements as constants, which caused some of his equations to be incorrect.
    [Show full text]
  • MTHE / MATH 237 Differential Equations for Engineering Science
    i Queen’s University Mathematics and Engineering and Mathematics and Statistics MTHE / MATH 237 Differential Equations for Engineering Science Supplemental Course Notes Serdar Y¨uksel November 26, 2015 ii This document is a collection of supplemental lecture notes used for MTHE / MATH 237: Differential Equations for Engineering Science. Serdar Y¨uksel Contents 1 Introduction to Differential Equations ................................................... .......... 1 1.1 Introduction.................................... ............................................ 1 1.2 Classification of DifferentialEquations ............ ............................................. 1 1.2.1 OrdinaryDifferentialEquations ................. ........................................ 2 1.2.2 PartialDifferentialEquations .................. ......................................... 2 1.2.3 HomogeneousDifferentialEquations .............. ...................................... 2 1.2.4 N-thorderDifferentialEquations................ ........................................ 2 1.2.5 LinearDifferentialEquations ................... ........................................ 2 1.3 SolutionsofDifferentialequations ................ ............................................. 3 1.4 DirectionFields................................. ............................................ 3 1.5 Fundamental Questions on First-Order Differential Equations ...................................... 4 2 First-Order Ordinary Differential Equations ..................................................
    [Show full text]
  • Math 6730 : Asymptotic and Perturbation Methods Hyunjoong
    Math 6730 : Asymptotic and Perturbation Methods Hyunjoong Kim and Chee-Han Tan Last modified : January 13, 2018 2 Contents Preface 5 1 Introduction to Asymptotic Approximation7 1.1 Asymptotic Expansion................................8 1.1.1 Order symbols.................................9 1.1.2 Accuracy vs convergence........................... 10 1.1.3 Manipulating asymptotic expansions.................... 10 1.2 Algebraic and Transcendental Equations...................... 11 1.2.1 Singular quadratic equation......................... 11 1.2.2 Exponential equation............................. 14 1.2.3 Trigonometric equation............................ 15 1.3 Differential Equations: Regular Perturbation Theory............... 16 1.3.1 Projectile motion............................... 16 1.3.2 Nonlinear potential problem......................... 17 1.3.3 Fredholm alternative............................. 19 1.4 Problems........................................ 20 2 Matched Asymptotic Expansions 31 2.1 Introductory example................................. 31 2.1.1 Outer solution by regular perturbation................... 31 2.1.2 Boundary layer................................ 32 2.1.3 Matching................................... 33 2.1.4 Composite expression............................. 33 2.2 Extensions: multiple boundary layers, etc...................... 34 2.2.1 Multiple boundary layers........................... 34 2.2.2 Interior layers................................. 35 2.3 Partial differential equations............................. 36 2.4 Strongly
    [Show full text]
  • Modified Variation of Parameters Method for Differential Equations
    World Applied Sciences Journal 6 (10): 1372-1376, 2009 ISSN 1818-4952 © IDOSI Publications, 2009 Modified Variation of Parameters Method for Differential Equations 1Syed Tauseef Mohyud-Din, 2Muhammad Aslam Noor and 2Khalida Inayat Noor 1HITEC University, Taxila Cantt, Pakistan 2Department of Mathematics, COMSATS Institute of Information Technology, Islamabad, Pakistan Abstract: In this paper, we apply the Modified Variation of Parameters Method (MVPM) for solving nonlinear differential equations which are associated with oscillators. The proposed modification is made by the elegant coupling of traditional Variation of Parameters Method (VPM) and He’s polynomials. The suggested algorithm is more efficient and easier to handle as compare to decomposition method. Numerical results show the efficiency of the proposed algorithm. Key words: Variation of parameters method • He’s polynomials • nonlinear oscillator INTRODUCTION the inbuilt deficiencies of various existing techniques. Numerical results show the complete reliability of the The nonlinear oscillators appear in various proposed technique. physical phenomena related to physics, applied and engineering sciences [1-6] and the references therein. VARIATION OF Several techniques including variational iteration, PARAMETERS METHOD (VPM) homotopy perturbation and expansion of parameters have been applied for solving such problems [1-6]. Consider the following second-order partial He [2-9] developed the homotopy perturbation differential equation: method for solving various physical problems. This reliable technique has been applied to a wide range of ytt = f(t,x,y,z,y,y,y,yxyzxx ,yyy ,y)zz (1) diversified physical problems [1-23] and the references therein. Recently, Ghorbani et al. [10, 11] introduced He’s polynomials by splitting the nonlinear term into a where t such that (-∞<t<∞) is time and f is linear or non series of polynomials.
    [Show full text]
  • Step Functions, Delta Functions, and the Variation of Parameters Formula
    STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA STEPHEN SCHECTER 1. The unit step function and piecewise continuous functions The Heaviside unit step function u(t) is given by 0 if t< 0, u(t)= (1 if t> 0. The function u(t) is not defined at t = 0. Often we will not worry about the value of a function at a point where it is discontinuous, since often it doesn’t matter. u(t) u(t−a) 1−u(t−b) u(t−a)−u(t−b) 1 1 1 1 a b a b Figure 1.1. Heaviside unit step function. The function u(t) turns on at t = 0. The function u(t − a) is just u(t) shifted, or dragged, so that it turns on at t = a. The function 1 − u(t − b) turns off at t = b. The function u(t − a) − u(t − b), with a < b, turns on a t = a and turns off at t = b. We can use the unit step function to crop and shift functions. Multiplying f(t) by u(t−a) crops f(t) so that it turns on at t = a: 0 if t < a, u(t − a)f(t)= (f(t) if t > a. Multiplying f(t) by u(t − a) − u(t − b), with a < b, crops f(t) so that it turns on at t = a and turns off at t = b: 0 if t < a, (u(t − a) − u(t − b))f(t)= f(t) if a<t<b, 0 if t > b.
    [Show full text]
  • Variation of Parameters
    Overview An Example Double Check Further Discussion Variation of Parameters Bernd Schroder¨ logo1 Bernd Schroder¨ Louisiana Tech University, College of Engineering and Science Variation of Parameters y = yp + yh 2. Variation of Parameters is a way to obtain a particular solution of the inhomogeneous equation. 3. The particular solution can be obtained as follows. 3.1 Assume that the parameters in the solution of the homogeneous equation are functions. (Hence the name.) 3.2 Substitute the expression into the inhomogeneous equation and solve for the parameters. Overview An Example Double Check Further Discussion Variation of Parameters 1. The general solution of an inhomogeneous linear differential equation is the sum of a particular solution of the inhomogeneous equation and the general solution of the corresponding homogeneous equation. logo1 Bernd Schroder¨ Louisiana Tech University, College of Engineering and Science Variation of Parameters 2. Variation of Parameters is a way to obtain a particular solution of the inhomogeneous equation. 3. The particular solution can be obtained as follows. 3.1 Assume that the parameters in the solution of the homogeneous equation are functions. (Hence the name.) 3.2 Substitute the expression into the inhomogeneous equation and solve for the parameters. Overview An Example Double Check Further Discussion Variation of Parameters 1. The general solution of an inhomogeneous linear differential equation is the sum of a particular solution of the inhomogeneous equation and the general solution of the corresponding homogeneous equation. y = yp + yh logo1 Bernd Schroder¨ Louisiana Tech University, College of Engineering and Science Variation of Parameters 3. The particular solution can be obtained as follows.
    [Show full text]
  • Asymptotic Analysis and Singular Perturbation Theory
    Asymptotic Analysis and Singular Perturbation Theory John K. Hunter Department of Mathematics University of California at Davis February, 2004 Copyright c 2004 John K. Hunter ii Contents Chapter 1 Introduction 1 1.1 Perturbation theory . 1 1.1.1 Asymptotic solutions . 1 1.1.2 Regular and singular perturbation problems . 2 1.2 Algebraic equations . 3 1.3 Eigenvalue problems . 7 1.3.1 Quantum mechanics . 9 1.4 Nondimensionalization . 12 Chapter 2 Asymptotic Expansions 19 2.1 Order notation . 19 2.2 Asymptotic expansions . 20 2.2.1 Asymptotic power series . 21 2.2.2 Asymptotic versus convergent series . 23 2.2.3 Generalized asymptotic expansions . 25 2.2.4 Nonuniform asymptotic expansions . 27 2.3 Stokes phenomenon . 27 Chapter 3 Asymptotic Expansion of Integrals 29 3.1 Euler's integral . 29 3.2 Perturbed Gaussian integrals . 32 3.3 The method of stationary phase . 35 3.4 Airy functions and degenerate stationary phase points . 37 3.4.1 Dispersive wave propagation . 40 3.5 Laplace's Method . 43 3.5.1 Multiple integrals . 45 3.6 The method of steepest descents . 46 Chapter 4 The Method of Matched Asymptotic Expansions: ODEs 49 4.1 Enzyme kinetics . 49 iii 4.1.1 Outer solution . 51 4.1.2 Inner solution . 52 4.1.3 Matching . 53 4.2 General initial layer problems . 54 4.3 Boundary layer problems . 55 4.3.1 Exact solution . 55 4.3.2 Outer expansion . 56 4.3.3 Inner expansion . 57 4.3.4 Matching . 58 4.3.5 Uniform solution . 58 4.3.6 Why is the boundary layer at x = 0? .
    [Show full text]
  • Who Solved the Bernoulli Differential Equation and How Did They Do It? Adam E
    Who Solved the Bernoulli Differential Equation and How Did They Do It? Adam E. Parker Adam Parker ([email protected]) is an associate professor at Wittenberg University in Springfield, Ohio. He was an undergraduate at the University of Michigan and received his Ph.D. in algebraic geometry from the University of Texas at Austin. He teaches a wide range of classes and often tries to incorporate primary sources in his teaching. This paper grew out of just such an attempt. Everyone loves a mystery; mathematicians are no exception. Since we seek out puzzles and problems daily, and spend so much time proving things beyond any reasonable doubt, we probably enjoy a whodunit more than the next person. Here’s a mystery to ponder: Who first solved the Bernoulli differential equation dy C P.x/y D Q.x/yn? dx The name indicates it was a Bernoulli, but which? Aren’t there 20 Bernoulli mathe- maticians? (Twenty is probably an exaggeration but we could reasonably count nine!) Or, as is so often the case in mathematics, perhaps the name has nothing to do with the solver. The culprit could be anyone! Like every good mystery, the clues contradict each other. Here are the prime suspects. Was it Gottfried Leibniz—the German mathematician, philosopher, and developer of the calculus? According to Ince [12, p. 22] “The method of solution was discovered by Leibniz, Acta Erud. 1696, p.145.” Or was it Jacob (James, Jacques) Bernoulli—the Swiss mathematician best known for his work in probability theory? Whiteside [21, p. 97] in his notes to Newton’s papers, states, “The ‘generalized de Beaune’ equation dy=dx D py C qyn was given its complete solution in 1695 by Jakob Bernoulli.” Or was it Johann (Jean, John) Bernoulli—Jacob’s acerbic and brilliant younger brother? Varignon [11, p.
    [Show full text]
  • Nonhomogeneous Equations and Variation of Parameters
    Nonhomogeneous Equations and Variation of Parameters June 17, 2016 1 Nonhomogeneous Equations 1.1 Review of First Order Equations If we look at a first order homogeneous constant coefficient ordinary differential equation by0 + cy = 0: then the corresponding auxiliary equation ar + c = 0 has a root r1 = −c=a and we have a solution r1t −ct=a yh(t) = ce = c1e If the equation is nonhomogeneous by0 + cy = f: Then, we introduce the integrating factor ect=b d (ect=by) = ect=bf dt Z ct=a ct=a e y(t) = c1 + e f(t)dt Z −ct=b −ct=b ct=b y(t) = c1e + e e f(t)dt y(t) = yh(t) + yp(t) The solution is a sum of 0 • yh(t), the solution to the homogeneous equation. (byh + cyh = 0). It has the constant that will be determined by the initial condition. • yp(t), a solution that involves f. Then, 0 0 0 b(yh + yp) + c(yh + yp) = (byh + cyh) + (byp + cyp) = 0 + f = f: 1 We next take an similar, but less formal approach to second order equations, writing, y = yh + yp where yh is a general solution to 00 0 ayh + byh + cyh = 0: and yp is a particular solution to 00 0 ayp + byp + cyp = f: 1.2 Examples We gain intuition in the nature of particular solution through some illustrative exampels Example 1. For y00 + y0 + 4y = 2t; we try a particular solution yp(t) = At + B. Then yp(t) = At +B 0 yp(t) = A 00 yp (t) = 0 4yp(t) = 4At +4B 0 yp(t) = A 00 yp (t) = 0 Thus, 4At + (A + 4B) = 2t 4A = 2;A + 4B = 0 1 1 A = B = − 2 8 and 1 1 y (t) = t − : p 3 8 we can turn this suggestion into a strategy for the case that f is a polynomial.
    [Show full text]
  • Variation of Parameters
    156 ● CHAPTER 4 HIGHER-ORDER DIFFERENTIAL EQUATIONS In Problems 27–34 find linearly independent functions that 55. yЉϩ25y ϭ 20 sin 5x 56. yЉϩy ϭ 4 cos x Ϫ sin x are annihilated by the given differential operator. 57. yЉϩyЈϩy ϭ x sin x 58. yЉϩ4y ϭ cos2x 27. D5 28. D2 ϩ 4D yٞϩ8yЉϭϪ6x2 ϩ 9x ϩ 2 .59 29. D Ϫ D ϩ 30. D2 Ϫ D Ϫ ( 6)(2 3) 9 36 Ϫ yٞϪyЉϩyЈϪy ϭ xex Ϫ e x ϩ 7 .60 31. D2 ϩ 5 32. D2 Ϫ 6D ϩ 10 yٞϪ3yЉϩ3yЈϪy ϭ ex Ϫ x ϩ 16 .61 33. D3 Ϫ 10D2 ϩ 25D 34. D2(D Ϫ 5)(D Ϫ 7) Ϫ 2yٞϪ3yЉϪ3yЈϩ2y ϭ (ex ϩ e x)2 .62 In Problems 35–64 solve the given differential equation by Ϫ ٞϩ Љϭ x ϩ (4) undetermined coefficients 63. y 2y y e 1 64. y(4) Ϫ 4yЉϭ5x2 Ϫ e2x 35. yЉϪ9y ϭ 54 36. 2yЉϪ7yЈϩ5y ϭϪ29 .yЉϩyЈϭ3 38. yٞϩ2yЉϩyЈϭ10 In Problems 65–72 solve the given initial-value problem .37 Љϩ Јϩ ϭ ϩ 39. y 4y 4y 2x 6 65. yЉϪ64y ϭ 16, y(0) ϭ 1, yЈ(0) ϭ 0 40. yЉϩ3yЈϭ4x Ϫ 5 66. yЉϩyЈϭx, y(0) ϭ 1, yЈ(0) ϭ 0 yٞϩyЉϭ8x2 42. yЉϪ2yЈϩy ϭ x3 ϩ 4x .41 67. yЉϪ5yЈϭx Ϫ 2, y(0) ϭ 0, yЈ(0) ϭ 2 43. yЉϪyЈϪ12y ϭ e4x 44. yЉϩ2yЈϩ2y ϭ 5e6x 68. yЉϩ5yЈϪ6y ϭ 10e2x, y(0) ϭ 1, yЈ(0) ϭ 1 45. yЉϪ2yЈϪ3y ϭ 4ex Ϫ 9 Ϫ 69.
    [Show full text]
  • Theory of Ordinary Differential Equations
    theoryofodes July 4, 2007 13:20 Page i ¨ © Theory of Ordinary Differential Equations ¨ ¨ © © ¨ © theoryofodes July 4, 2007 13:20 Page ii ¨ © ¨ ¨ © © ¨ © theoryofodes July 4, 2007 13:20 Page i ¨ © Theory of Ordinary Differential Equations Christopher P. Grant Brigham Young University ¨ ¨ © © ¨ © theoryofodes July 4, 2007 13:20 Page ii ¨ © ¨ ¨ © © ¨ © theoryofodes July 4, 2007 13:20 Page i ¨ © Contents Contents i 1 Fundamental Theory 1 1.1 ODEsandDynamicalSystems . 1 1.2 ExistenceofSolutions . ... ... ... ... .. ... ... 6 1.3 UniquenessofSolutions . 10 1.4 Picard-LindelöfTheorem. 14 ¨ ¨ © 1.5 IntervalsofExistence. 16 © 1.6 DependenceonParameters . 20 2 Linear Systems 27 2.1 ConstantCoefficientLinearEquations . 27 2.2 UnderstandingtheMatrixExponential . 30 2.3 GeneralizedEigenspaceDecomposition . 33 2.4 OperatorsonGeneralizedEigenspaces . 37 2.5 RealCanonicalForm .. ... ... ... ... .. ... ... 41 2.6 SolvingLinearSystems. 43 2.7 QualitativeBehaviorofLinearSystems . 50 2.8 ExponentialDecay. ... ... ... ... ... .. ... ... 54 2.9 NonautonomousLinearSystems . 56 2.10NearlyAutonomousLinearSystems. 61 2.11PeriodicLinearSystems . 65 3 Topological Dynamics 71 3.1 InvariantSetsandLimitSets . 71 3.2 RegularandSingularPoints. 75 i ¨ © theoryofodes July 4, 2007 13:20 Page ii ¨ © Contents 3.3 DefinitionsofStability . 80 3.4 PrincipleofLinearizedStability . 85 3.5 Lyapunov’sDirectMethod. 90 3.6 LaSalle’sInvariancePrinciple . 94 4 Conjugacies 101 4.1 Hartman-GrobmanTheorem: Part1 . 101 4.2 Hartman-GrobmanTheorem: Part2 . 103 4.3 Hartman-GrobmanTheorem: Part3
    [Show full text]
  • Mathematical Methods of Physics I Instructor: Predrag Cvitanovi´C Fall Semester 2012
    Georgia Tech PHYS 6124 Mathematical Methods of Physics I Instructor: Predrag Cvitanovi´c Fall semester 2012 Homework Set #7 due October 30 2012 == show all your work for maximum credit, == put labels, title, legends on any graphs == acknowledge study group member, if collective effort [All problems in this set are from Goldbart] Problem 1) Motion of a classical particle Consider a classical particle of unit mass moving along the x-axis. Suppose that the motion is free, except that at time t = t, with 0 < t < T, the particle receives an impulse of unit strength. a) Write Newton’s equation describing the motion of the particle. Suppose that at time t = 0 the particle is located at position x1 and that at time t = T it is located at position x2. b) Sketch the position, velocity and acceleration of the particle as a function of time for 0 < t < T. c) Compute the Green function for the motion of the particle, i.e., solve d2G(t, t)/dt2 = d(t − t). d) Consider the applied force f (t) (with 0 ≤ t ≤ T) to be a sequence of impulses. Hence establish the trajectory of the particle in terms of an integral over the applied force. e) Suppose that the force takes the form f (t) = t2/2. Find the motion of the particle. f)( optional) Rather than solve for the Green function directly, as you did in part (c), construct the Green function using the eigenfunction expansion technique. Show, by Fourier analysis, that the two schemes for comput- ing the Green function give equivalent results.
    [Show full text]