First order differential equation Linear equation; Method of integrating factors

Similar documents
SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.

Thomas Whitham Sixth Form

The Matrix Exponential

The Matrix Exponential

Thomas Whitham Sixth Form

MSLC Math 151 WI09 Exam 2 Review Solutions

y cos x = cos xdx = sin x + c y = tan x + c sec x But, y = 1 when x = 0 giving c = 1. y = tan x + sec x (A1) (C4) OR y cos x = sin x + 1 [8]

Mathematics 1110H Calculus I: Limits, derivatives, and Integrals Trent University, Summer 2018 Solutions to the Actual Final Examination

u r du = ur+1 r + 1 du = ln u + C u sin u du = cos u + C cos u du = sin u + C sec u tan u du = sec u + C e u du = e u + C

NEW APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA

First derivative analysis

Calculus concepts derivatives

Exercise 1. Sketch the graph of the following function. (x 2

Differentiation of Exponential Functions

Differential Equations

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula

Finite Element Analysis

Schematic of a mixed flow reactor (both advection and dispersion must be accounted for)

a 1and x is any real number.

Probability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers Roy D. Yates and David J.

Indeterminate Forms and L Hôpital s Rule. Indeterminate Forms

dx equation it is called a second order differential equation.

MATH 1080 Test 2-SOLUTIONS Spring

A Propagating Wave Packet Group Velocity Dispersion

The second condition says that a node α of the tree has exactly n children if the arity of its label is n.

[1] (20 points) Find the general solutions of y y 2y = sin(t) + e t. Solution: y(t) = y c (t) + y p (t). Complementary Solutions: y

Grade 12 (MCV4UE) AP Calculus Page 1 of 5 Derivative of a Function & Differentiability

Addition of angular momentum

INTEGRATION BY PARTS

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory

MA1506 Tutorial 2 Solutions. Question 1. (1a) 1 ) y x. e x. 1 exp (in general, Integrating factor is. ye dx. So ) (1b) e e. e c.

64. A Conic Section from Five Elements.

Modeling with first order equations (Sect. 2.3).

Differential Equations

1973 AP Calculus AB: Section I


Function Spaces. a x 3. (Letting x = 1 =)) a(0) + b + c (1) = 0. Row reducing the matrix. b 1. e 4 3. e 9. >: (x = 1 =)) a(0) + b + c (1) = 0

Limiting value of higher Mahler measure

Basic Polyhedral theory

Addition of angular momentum

6.1 Integration by Parts and Present Value. Copyright Cengage Learning. All rights reserved.

ECE602 Exam 1 April 5, You must show ALL of your work for full credit.

A. Limits and Horizontal Asymptotes ( ) f x f x. f x. x "±# ( ).

Y 0. Standing Wave Interference between the incident & reflected waves Standing wave. A string with one end fixed on a wall

( ) Differential Equations. Unit-7. Exact Differential Equations: M d x + N d y = 0. Verify the condition

Math 34A. Final Review

On the irreducibility of some polynomials in two variables

VTU NOTES QUESTION PAPERS NEWS RESULTS FORUMS

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.

The graph of y = x (or y = ) consists of two branches, As x 0, y + ; as x 0, y +. x = 0 is the

(1) Then we could wave our hands over this and it would become:

Solution: APPM 1360 Final (150 pts) Spring (60 pts total) The following parts are not related, justify your answers:

Multiple Short Term Infusion Homework # 5 PHA 5127

SPH4U Electric Charges and Electric Fields Mr. LoRusso

y = 2xe x + x 2 e x at (0, 3). solution: Since y is implicitly related to x we have to use implicit differentiation: 3 6y = 0 y = 1 2 x ln(b) ln(b)

2008 AP Calculus BC Multiple Choice Exam

Notes on Differential Geometry

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

Instructions for Section 1

PROBLEM SET Problem 1.

Lie Groups HW7. Wang Shuai. November 2015

Calculus Revision A2 Level

Quasi-Classical States of the Simple Harmonic Oscillator

DIFFERENTIAL EQUATION

1. I = 2 3. I = 4 5. I = I = 5 2

Calculus II (MAC )

CS 491 G Combinatorial Optimization

Construction of asymmetric orthogonal arrays of strength three via a replacement method

Appendix 2.3 General Solutions for the Step Response of Third- and Fourth-Order Systems (with some unpleasant surprises!)

A RELATIVISTIC LAGRANGIAN FOR MULTIPLE CHARGED POINT-MASSES

Deift/Zhou Steepest descent, Part I

Problem Set 6 Solutions

MATH 319, WEEK 15: The Fundamental Matrix, Non-Homogeneous Systems of Differential Equations

EXST Regression Techniques Page 1

The Transfer Function. The Transfer Function. The Transfer Function. The Transfer Function. The Transfer Function. The Transfer Function

That is, we start with a general matrix: And end with a simpler matrix:

Chapter 1. Chapter 10. Chapter 2. Chapter 11. Chapter 3. Chapter 12. Chapter 4. Chapter 13. Chapter 5. Chapter 14. Chapter 6. Chapter 7.

Where k is either given or determined from the data and c is an arbitrary constant.

Engineering 323 Beautiful HW #13 Page 1 of 6 Brown Problem 5-12

Walk Like a Mathematician Learning Task:

Derangements and Applications

COMPUTER GENERATED HOLOGRAMS Optical Sciences 627 W.J. Dallas (Monday, April 04, 2005, 8:35 AM) PART I: CHAPTER TWO COMB MATH.

1 1 1 p q p q. 2ln x x. in simplest form. in simplest form in terms of x and h.

Propositional Logic. Combinatorial Problem Solving (CPS) Albert Oliveras Enric Rodríguez-Carbonell. May 17, 2018

BSc Engineering Sciences A. Y. 2017/18 Written exam of the course Mathematical Analysis 2 August 30, x n, ) n 2

Partial Derivatives: Suppose that z = f(x, y) is a function of two variables.

Sec 2.3 Modeling with First Order Equations

Bifurcation Theory. , a stationary point, depends on the value of α. At certain values

10. Limits involving infinity

1 General boundary conditions in diffusion

cycle that does not cross any edges (including its own), then it has at least

Section 11.6: Directional Derivatives and the Gradient Vector

Abstract Interpretation: concrete and abstract semantics

Note If the candidate believes that e x = 0 solves to x = 0 or gives an extra solution of x = 0, then withhold the final accuracy mark.

2.3 Matrix Formulation

Integration by Parts

AP PHYSICS C: ELECTRICITY AND MAGNETISM 2015 SCORING GUIDELINES

Introduction to Condensed Matter Physics

Homogeneous Constant Matrix Systems, Part I

Chapter 10. The singular integral Introducing S(n) and J(n)

Transcription:

First orr iffrntial quation Linar quation; Mtho of intgrating factors Exampl 1: Rwrit th lft han si as th rivativ of th prouct of y an som function by prouct rul irctly. Solving th first orr iffrntial quation From calculus for iffrntiating a prouct, w hav (4 + t ) y + ty = 4t. (0.1) (4 + t ) y + ty = ((4 + t )y); it follows that th quation (0.1) can b rwrit as ((4 + t )y) = 4t. (0.) Thus, vn though y is known, w can intgrat both sis of quation (0.) with rspct to t, thrby obtaining (4 + t )y = t + c, whr c is an arbitrary constant of intgration. Solving for y, w fin that y = t 4 + t + c 4 + t. This is th gnral solution of quation (0.1). Unfortunatly, most first orr linar iffrntial quation cannot b solv as illustrat in Exampl 1 bcaus thir lft han si ar not th rivativ of th prouct of y an som othr function. Howvr, Libniz iscovr that if th iffrntial quation is multipli by a crtain function µ(t), thn th quation is convrt into on that is immiatly intgrabl by using th prouct rul of rivativs, just as Exampl 1. Th function µ(t) is call an intgrating factor an our main task in this sction is to trmin how to fin it for a givn quation. Exampl : Multipli µ(t) an thn appli prouct rul. y + 1 y = 1 t 3. (0.3) Th first stp is to multiply quation (0.3) by a function µ(t), as yt untrmin; thus µ(t) y + 1 µ(t)y = 1 µ(t)t/3. (0.4) Th qustion now is whthr w can choos µ(t) such that th lft han si of abov quation is th rivativ of th prouct µ(t)y. For any iffrntial function µ(t) w hav (µ(t)y) = µ(t)y + µ(t) y. (0.5) Thus th lft han si of quation (0.4) an th right han si of quation (0.5) ar intical, provi w choos µ(t) to satisfy µ(t) = 1 µ(t). (0.6) 1

Our sarch for an intgrating factor will b succssful if w can fin a solution of quation (0.6). Rform quation (0.6) to µ µ = 1, intgrat on both si, w can gt that which givs log µ(t) = 1 t + c, µ(t) = c t/. (0.7) Th function µ(t) givn by (0.7) is an intgrating factor for quation (0.3). Sinc w o not n th most gnral intgrating factor, w will choos c to b 1 in quation (0.7) an us µ(t) = t/. Now w rturn to quation (0.3), multiply it by th intgrating factor t/, an obtain t/ y + 1 t/ y = 1 5t/6. (0.8) By th choic w hav ma of th intgrating factor, th lft han si of quation (0.8) is th rivativ of t/ y, so that quation (0.8) bcoms (t/ y) = 1 5t/6. By intgrating both sis of th abov quation, w obtain t/ y = 3 5 5t/6 + c, solving for y, w gt th gnral solution of quation (0.3), namly, y = 3 5 t/3 + c t/. Lt us xtn th mtho of intgrating factors to quations of th form y + ay = g(t). (0.9) Whr a is a givn constant an g(t) is a givn function. Procing as in Exampl, w fin that th intgrating factor µ(t) must satisfy µ = aµ. Thus th intgrating factor is µ(t) = at. Multiply th quation (0.9) by µ(t), w obtain at y + aat y = at g(t), or (at y) = at g(t). Intgrating on both sis of abov quation, w fin that at y = at g(t) + c,

whr c is an arbitrary constant. For many simpl function g(t), w can valuat th intgral in th abov quation an xprss th solution y in trms of lmntary functions, as in Exampl. Howvr, for mor complicat function g(t), it is ncssary to lav th solution in intgral form. In this cas, y = at t t 0 as g(s)s + c at, not that in quation abov w hav us s to not th intgration variabl to istinguish it from th inpnnt variabl t an w hav chosn convnint valu t 0 as th lowr limit if intgration. Th choic of t 0 trmins th spcific valu of th constant c but os not chang th solution. For xampl, plugging t = t 0 into th gnral solution formula shows that c = y(t 0 ) at0. Exampl 3: Fin th gnral solution of th iffrntial quation y y = 4 t (0.10) an plot th graphs of svral solutions. Discuss th bhavior of solutions as t. Equation (0.10) if of th form with a = ; thrfor th intgrating factor is µ(t) = t. Multiply th iffrntial quation (0.10) by µ(t), w obtian t y t y = 4 t t t, or ( t y) = 4 t t t. (0.11) Thn, by intgrating both sis of this quation, w hav t y = t + 1 t t + 1 4 t + c, whr w hav us intgration by parts on th last trm in quation. Thus th gnral solution of quation (0.10) is y = 7 4 + 1 t + ct. (0.1) Th bhavior of th solution for larg valus of t is trmin by th trm c t. If c 0, thn th solution grows xponntially larg in magnitu with th sam sign as c itslf. Thus th solutions ivrg as t bcoms larg. Th bounary btwn solutions that ultimatly grow positivly an thos that grow ngativly occurs whn c = 0. If c = 0 an thn st t = 0, w gt that y = 7 4 is th sparation point on th y axis. Not that for this initial valu, th solution is y = 7 4 + 1 t; it grows positivly, but linarly rathr xponntially. Th gnral first-orr linar iffrntial quation. y + p(t)y = g(t), whr p, g ar givn functions. To trmin an appropriat factor, w multiply quation by an as yt untrmin function µ(t), obtaining µ(t) y + p(t)µ(t)y = µ(t)g(t). (0.13) Follow th sam lin of vlopmnt as in Exampl, w s that h lft han si of quation (0.13) is th rivativ of th prouct µ(t)y, provi that µ(t) satisfy th quation µ(t) = p(t)µ(t). (0.14) 3

If w assum tmporarily that µ(t) is positiv, thn w hav 1 µ(t) = p(t) µ(t) an consquntly log µ(t) = p(t) + c. By choosing th arbitrary constant c to b zro, w obtain th simplst possibl function for µ, namly, µ(t) = xp( p(t)). Not that µ(t) is positiv for all t, as w assum. Rturning to quation (0.13), w hav Hnc (µ(t)y) = µ(t)y. (0.15) µ(t)y = µ(t)g(t) + c, (0.16) whr c is an arbitrary constant. Somtims th intgral in quation (0.16) can b valuat in trms of lmntary functions. Howvr, in gnral this is not possibl, so th gnral solution of quation is y = 1 t µ(t) ( µ(s)g(s)s + c), (0.17) t 0 whr again t 0 is som convnint lowr limit of intgration. Obsrv that quation (0.17) involvs two intgrations, on to obtain µ(t) an th othr to trmin y from quation (0.17). Exampl 4: Solv th initial valu problm { ty + y = 4t, y(1) =. (0.18) In orr to trmin p(t) an g(t) corrctly, w must first rwrit quation in th stanar form. Thus w hav y + y = 4t, (0.19) t so p(t) = /t an g(t) = 4t. To solv quation (0.19), w first comput th intgrating factor µ(t) µ(t) = xp( t ) = log t = t. On multiplying quation (0.19) by µ(t) = t, w obtain t y + ty = (t y) = 4t 3, an thrfor t y = 4t 3 = t 4 + c, whr c is an arbitrary constant. It follows that, for t > 0, is th gnral solution of quation from (0.18). y = t + c t (0.0) 4

To satisfy initial conition, st t = 1 an y = in quation (0.1), so c = 1; thus is th solution of th initial valu problm (0.18). It is important to not that whil th function y = t + 1 t, t > 0 (0.1) y = t + 1 t, t < 0 is part of th gnral solution of quation, it is not part of th solution of this initial valu problm. This is th first xampl in which th solution fails to xist for som valus of t. Again, this is u to th infinit iscontinuity in p(t) at t = 0, which rstrict th solution to th intrval 0 < t <. Draw th pictur for y = t + c t for iffrnt valu c. Exampl 5: Solv th initial valu problm y + ty =, (0.) y(0) = 1. (0.3) To convrt th iffrntial quation (0.) to th stanar form, w hav y + t y = 1. (0.4) Thus p(t) = t an g(t) = 1, th intgrating factor is µ(t) = xp(t /4). Thn multiply quation (0.4) by µ(t), so that t /4 y + t t /4 y = t /4. (0.5) hnc Intgrating on both sis, w obtain (t /4 y) = t /4. t /4 y = t /4 + c. (0.6) Th intgral on th right han si cannot b valuat by lmntary function, so w lav th intgral unvaluat. By choosing th lowr limit of th intgration as th initial point t = 0, w can rplac th quation (0.6) by t /4 y = t 0 s /4 s + c. (0.7) whr c is an arbitrary constant. It thn follows that th gnral solution of quation (0.) is givn by y = t /4 t By th initial conition, w can trmin that c = 1. 0 s /4 s + c t /4. (0.8) 5