Size: px
Start display at page:

Download ""

Transcription

1 Part 1: Overview of Ordinary Dierential Equations 1

2

3 Chapter 1 Basic Concepts and Problems 1.1 Problems Leading to Ordinary Dierential Equations Many scientic and engineering problems are modeled by systems of ordinary dierential equations (ODEs). Some examples [] follow. 1. Mechanical Vibrations. Let y(t) denote the displacement at time t of a block of mass m that is connected to a spring of stiness k, a damper of resistance c, and an oscillator f(t) (Figure 1.1.1). If the system is released from position y 0 velocity y 0 then its subsequent motion satises 0 with a m d y dt + c dy dt + ky = f(t) t>0 y(0) = y 0 dy(0) dt = y 0 0 : This is an example of an initial value problem (IVP) for a second-order linear ODE. The variable y(t) is called: The variable t is called: The ODE is second order because: It's linear because: This is an IVP because:

4 k m c f(t) y(t) Figure 1.1.1: Motion of a spring-mass-damper system.. Ecology. Consider a population of predators and prey living in an ecological niche. The predators survive by eating the prey and the prey exist on an independent source of food. A classical situation involves foxes and rabbits. Let P (t) and p(t), respectively, denote the populations of predators and prey at time t. Given the initial populations P 0 and p 0 of predators and prey, their subsequent populations satisfy the Lotka-Volterra equations dp dt = p(a ; P ) t>0 p(0) = p 0 dp dt = P (;c + p) t>0 P (0) = P 0 where a, c,, and are positive constants corresponding to the prey's natural growth rate, the predator's natural death rate, the prey's death rate upon coming into contact with predators, and the predator's growth rate upon coming into contact with prey. This example involves an IVP for a system of two rst-order nonlinear ODEs. The ODEs are nonlinear because: The system is rst order because:. Column Buckling. The lateral displacement y(t) at position t of a clamped-hinged bar of length l that is subjected to a load P (Figure 1.1.) may be approximated by the Euler-Bernoulli equations d 4 y dt 4 with the boundary conditions y(0) = dy(0) dt + d y =0 0 <t<l dt =0 y(l) = d y(l) dt =0: 4

5 The parameter = P=EI, where EI is the exural rigidity ofthe bar. This is an example of a boundary value problem (BVP) for a fourth-order linear ODE. This is a boundary value problem because: Observe that y(t) = 0 is a solution of this problem for all values of and l. A more interesting problem is to determine y and those values of and l for which non-trivial (y(t) 6= 0) solutions exist. As such, this BVP is also a dierential eigenvalue problem. Nontrivial solutions y(t) are called eigen functions and their corresponding values of are eigenvalues. P t y(t) P Figure 1.1.: Buckling of an elastic column. 4. Pendelum Oscillations ([1], Section 1.). The position (x(t) y(t)) at time t of a particle of mass m oscillating on a pendelum of length l (Figure 1.1.) is m d x dt = ;T sin = ; T l y x md = mg ; T cos = mg ; T y t>0 dt l where g is the acceleration of gravity, T (t) is the tension in the string, and (t) is the angle of the pendelum relative to the vertical at time t (Figure 1.1.). These equations are, however, insucient to guarantee that the particle stays on the string. To ensure that this is so, we must supplement the ODEs by the algebraic constraint x + y = l : The two second-order dierential equations and the constraint comprise a system of ve dierential algebraic equations (DAEs) for the unknowns x(t), y(t), and T (t). Initial conditions specify x(0), dx(0)=dt, y(0), and dy(0)=dt, but not T (0). 5

6 Of course, for this problem, it's easy to eliminate the constraint byintroducing the change of variables x = l sin, y = l cos. This would reduce the DAEs to the second-order ODE d = ; g sin : dt l Such simplications would not be possible with more dicult systems. θ y T x mg Figure 1.1.: Oscillations of a simple Pendelum. IVPs will comprise our initial study (Part of these notes). We'll take up BVPs next (Part ) and conclude with a study of DAEs (Part 4). In almost all cases, it will suce to develop and analyze numerical methods for IVPs and BVPs that are written as rst-order vector systems of ODEs in the explicit form y 0 (t) =f(t y) t>0 (1.1.1a) where y(t) = 6 4 y 1 (t) y (t). y m (t) f(t y) = f 1 (t y 1 ::: y m ) f (t y 1 ::: y m ). f m (t y 1 ::: y m ) 5 (1.1.1b) and y 0 := dy=dt. Data for IVPs would is specied as y(0) = y 0 = 6 4 y 10 y 0. 5 : (1.1.1c) y m0 6

7 Boundary data is more complex. The most general boundary conditions specify a nonlinear relationship between y(0) and y(l) oftheform g 1 (y(0) y(l)) g (y(0) y(l)) g(y(0) y(l)) = 6 = 0: (1.1.1d) 4. 5 g m (y(0) y(l)) In many cases the extension of a scalar to a vector IVP will be obvious and we can study methods for a scalar IVP y 0 (t) =f(t y) t>0 y(0) = y 0 (1.1.) Example Higher-order ODEs can be written as rst-order systems. Consider the m th order equation z (m) = g(t z z 0 ::: z (m;1) ) where z (i) := d i z=dt i, and introduce the new variables y 1 = z y = z 0 y = z 00 ::: y m = z (m;1) : Then, we have a system in the form (1.1.1) with y 0 1 = y y 0 = y ::: y 0 m = g(t y 1 y ::: y m ): For an IVP, that data would prescribe as z (i) (0) = c i i =0 1 ::: m; 1 with the understanding that z (0) = z. Written in terms of y, we have y 1 (0) = c 0 y (0) = c 1 ::: y m (0) = c m;1 : With DAEs and sometimes with IVPs and BVPs, it will be necessary to consider implicit dierential systems F(t y y 0 )=0: (1.1.) If the 0 is nonsingular then, by the implicit function theorem, (1.1.) is equivalent to (1.1.1) however, it may be natural to solve (1.1.) in its implicit form to maintain, e.g., sparsity. DAEs involve systems 0 is singular.

8 Example Let us write the oscillating pendelum problem as a rst-order system by letting y 1 = x y = x 0 y = y y 4 = y 0 : Then, we have y 0 1 = y y 0 = ; T ml y 1 y 0 = y 4 y 0 4 = g ; T ml y y 5 = T with the constraint y 1 + y = l : This can be written in the form (1.1.) with y 0 ; y 1 y 0 + y 1y 5 F(t y y 0 ml )= 6 y 0 ; y 4 4 y 0 ; g + y y 5 4 ml y + y ; l 1 5 = 0: The Jacobian is clearly singular. F y = Problems 5 1. Letting y 1 = p and y = P, write the predator-prey model dp dt = p(a ; P ) t>0 p(0) = p 0 dp dt = P (;c + p) t>0 P (0) = P 0 in the vector form (1.1.1).. Write the mechanical vibration problem m d y dt in the form (1.1.1a-c) + c dy dt + ky = f(t) t>0 y(0) = y 0 8 dy(0) dt = y 0 0

9 . Write the column buckling problem d 4 y dt 4 + d y =0 0 <t<l y(0) = dy(0) y(l) =0 y(l) = d =0 dt dt dt in the form (1.1.1a,b,d) 1. Existence, Uniqueness, and Stability of IVPs Before considering numerical procedures for (1.1.1, 1.1.), let us review some results from ODE theory that ensure the existence of unique solutions which depend continuously on the initial data. For simplicity, letusfocus on a scalar IVP having the form (1.1.). Denition A function f(t y) satises a Lipschitz condition in a domain D if there exists a non-negative constant L such that jf(t y) ; f(t z)j Ljy ; zj 8(t y) (t z) D: (1..4) Satisfying a Lipschitz condition guarantees that solutions y(t) of (1.1.) are unique as expressed by the following theorem. Theorem Let f(t y) be continuous and satisfy a Lipschitz condition on f(t y) j 0 t T ;1 <y<1g: Then the IVP (1.1.) has a unique continuously dierentiable solution on [0 T] for all y 0 (;1 1). Proof. The proof appears in most books on ODE theory, e.g., []. Remark 1. A Lipschitz conditions guarantees unique solutions and is not needed for esistence. For example, the IVP y 0 =y = y(0) = 0 has the two solution y(t) = 0 and y(t) = t. This f(t y) = y = Lipschitz condition. does not satisfy a 9

10 Remark. Theorem 1..1 applies when f satisfys a Lipschitz condition on a compact, rather than an unbounded, domain D as long as the solution y(t) remains in D []. In addition to existence and uniqueness, we will want to know something about the stability of solutions of the IVP. In particular, we will usually be interested in the sensitivity of the solution to small changes in the data. Perturbations arise naturally in numerical computation due to discretization and round o errors. A formal study of sensitivity would lead us to the following notion of a well-posed system. Denition 1... The IVP (1.1.) is well-posed if there exists positive constants k and ^ such that, for any ^, the perturbed IVP z 0 = f(t z)+(t) z(0) = y (1..5) satises jy(t) ; z(t)j k (1..6) whenever j 0 j <and j(t)j <for t [0 T]. Again, a Lipschitz condition ensures that we are dealing with a well-posed IVP. Theorem 1... If f(t y) satises a Lipschitz condition on f(t y) j 0 t T ;1 <y<1g then y 0 = f(t y) is well posed on [0 T] with respect to all initial data. Proof. Let (t) =z(t) ; y(t) and subtract (1.1.) from (1..5) to obtain 0 (t) =f(t z) ; f(t y)+(t) (0) = 0 : Taking an absolute value and using the Lipschitz condition (1..4) j 0 (t)j Lj(t)j + j(t)j j(0)j = 0 : 10

11 We easily see that j(t)j 0 j 0 (t)j provided that j(t)j 0 exists and it is fairly easy to show that j(t)j 0 exists. Additionally, by assumption, max j(t)j < 0tT max j 0 j < 0tT so j(t)j 0 Lj(t)j + j(0)j <: Multiply by the integrating factor e ;Lt to obtain (e ;Lt j(t)j) 0 e ;Lt j(0)j <: Integrating j(t)j Since the denominator is smallest when t = 0 jz(t) ; y(t)j thus, y 0 = f(t y) iswell posed on [0 T]. L[(L +1)e Lt ; 1] : L = k 8t [0 T] The notion of awell-posed system is related to the more common notion of stability as indicated by the following denition. Denition 1... Consider the dierential equation y 0 = f(t y) and (without loss of generality) let the origin (0 0) be anequilibrium point, i.e., f(0 0) = 0. Then, the origin is: stable if a perturbation of the initial condition jy(0)j <grows no larger than for subsequent times, i.e., if jy(t)j <for t>0. asymptotically stable if it is stable and jy(0)j <implies that lim t!1 jy(t)j =0. unstable if it is not stable. Remark. This denition could, like Denition 1.., also involve perturbations of f(t y). We have omitted these for simplicity. 11

12 An autonomous system is one where f(t y) does not explicitly depend on t, i.e., f(t y) =f(y). If (0 0) is an equilibrium point then, in this case, f(0) = 0. Expanding the solution in a Taylor's series in y, we have y 0 (t) =f(y) =f(0) + f y (0)y(t)+O(y ): Since f(0) = 0, y 0 (t) =f y (0)y(t)+O(y ) where f y =@y. Recall that a function g(y) iso(y p ) if there exists a constant C > 0 such that jg(y)j Cy p as y! 0. Letting = f y (0), we see that the stability of an autonomous system is related to that of the simple linear IVP y 0 = y t>0 y(0) = : (1..) Additional details on autonomous and non-autonomous systems appear in Birkho and Rota [] or Hairer et al. [5], Section 1.1. The solution of the linear IVP (1..) is y(t) =e t hence, (1..) is stable when Re() 0, asymptotically stable when Re() < 0, and unstable when Re() > 0. These conclusions remain true for the original nonlinear autonomous problem when Re() 6= 0 and y is small enough for the O(y ) term to be negligible relative toy. This cannot happen in the stable case (Re() = 0) hence, it requires a more careful analysis. Remark 4. The analysis of non-autonomous systems is similar [1]. Analyzing the stability of autonomous vector systems y 0 = f(y(t)) (1..8) 1

13 is more complicated. Once again, assume that y = 0 is an equilibrium point and expand f(0) in a series y 0 = f(0)+f y (0)y + O(kyk )=f y (0)y + O(kyk ): We need a brief digression for a few denitions. Denition The Jacobian matrix of a vector-valued function f(y) with respect to y is the matrix f = @ym Denition The norm of a vector y is a scalar kyk such that 1. kyk 0 and kyk =0if and only if y = 0,. kyk = jjkyk for any scalar, and. kx + yk kxk + kyk. We'll identify specic vector norms when they are needed. 5 : (1..9) For the moment, let's return to our stability analysis and let A = f y (0). Then, the stability of the nonlinear autonomous system (1..8) is related to that of the linear system The solution of this system is y 0 = Ay y(0) = y 0 ky 0 k: (1..10a) y(t) =e At y 0 (1..10b) where the matrix exponential is dened by the series Often, A can be diagonalized as e At = I + ta + t + t + : (1..10c)! A! A 1 T ;1 AT = = m 5 (1..11) 1

14 where 1 ::: m are the eigenvalues of A and the columns of T are the corresponding eigenvectors. In this case, we may easily verify that the solution of (1..10) is y(t) =Te t T ;1 y 0 : Thus, (1..10) is stable if all of the eigenvalues have non-positive real parts, asymptotically stable if all of the eigenvalues have negative real parts, and unstable otherwise. Unfortunately, not all matrices are diagonalizable. However, A can always be reduced to the Jordan canonical form 1 T ;1 AT = 6 (1..1a) l where each Jordan block has the form i = 6 4 i 1 i (1..1b) i The dimension of the Jordan block i corresponds to the multiplicity ofthe eigenvalue i. Thus, if i is simple, the block is a scalar. With this, it is relatively easy to show [] that y = 0 is stable when either Re( i ) < 0orRe( i )=0and i is simple, i =1 ::: m, asymptotically stable when Re( i ) < 0, i =1 ::: m, and unstable otherwise. Example Consider the predator-prey model of Section 1.1 y 0 1 = y 1 (a ; y ) y 0 = y (;c + y 1 ) where a,, c, and are positive constants. This autonomous problem has two equilibrium points: y1 y = c= a= :

15 The point [0 0] T corresponds to extinction of predators and prey, whereas [c= a=] T implies a co-existence of both species. The stability of[0 0] T the linear system 0 y1 a 0 = y 0 ;c y1 y : may be analyzed by using Since the system matrix is diagonal, its eigenvalues are 1 = a, = ;c. Thus, 1 > 0 and [0 0] T is unstable. In order to analyze the stability of the second equilibrium point, we introduce the change of variables z 1 = y 1 ; c z = y ; a to move the equilibrium point to the origin. Substituting this transformation into the predator-prey equations gives z 0 1 =(z 1 + c )[a ; (z + a )] = ;z ( c + z 1) z 0 =(z + a )[;c + (z 1 + c )] = z 1( a + z ): The linearized system is now 0 z1 = z 0 ;c= a= 0 z1 z : The eigenvalues of this system are 1 = p i c=a. Since the eigenvalues are imaginary, the linear system is stable however, the stability of the nonlinear system is undetermined without additional analysis. Problems 1. =@y is continuous on a closed convex domain D. Show that f(t y) satises a Lipschitz condition on D with (Hint: use the Mean Value y) L = max j j: (t Are the following IVPs well posed ([4], p. )? Explain. y 0 = p 1 ; y y(0) = 0 y 0 = p y ; 1 y(0) = : 15

16 16

17 Bibliography [1] U.M. Ascher and L.R. Petzold. Computer Methods for Ordinary Dierential Equations and Dierential-Algebraic Equations. SIAM, Philadelphia, [] G. Birkho and G.-C. Rota. Ordinary Dierential Equations. John Wiley and Sons, New York, third edition, 198. [] W.E. Boyce and R.C. DiPrima. Elementary Dierential Equations and Boundary Value Problems. John Wiley and Sons, New York, third edition, 19. [4] C.W. Gear. Numerical Initial Value Problems in Ordinary Dierential Equations. Prentice Hall, Englewood Clis, 191. [5] E. Hairer, S.P. Norsett, and G. Wanner. Solving Ordinary Dierential Equations I: Nonsti Problems. Springer-Verlag, Berlin, second edition,

Math 266: Ordinary Differential Equations

Math 266: Ordinary Differential Equations Math 266: Ordinary Differential Equations Long Jin Purdue University, Spring 2018 Basic information Lectures: MWF 8:30-9:20(111)/9:30-10:20(121), UNIV 103 Instructor: Long Jin (long249@purdue.edu) Office

More information

APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS

APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS LECTURE 10 APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS Ordinary Differential Equations Initial Value Problems For Initial Value problems (IVP s), conditions are specified

More information

MATH 4B Differential Equations, Fall 2016 Final Exam Study Guide

MATH 4B Differential Equations, Fall 2016 Final Exam Study Guide MATH 4B Differential Equations, Fall 2016 Final Exam Study Guide GENERAL INFORMATION AND FINAL EXAM RULES The exam will have a duration of 3 hours. No extra time will be given. Failing to submit your solutions

More information

Dynamical Systems. August 13, 2013

Dynamical Systems. August 13, 2013 Dynamical Systems Joshua Wilde, revised by Isabel Tecu, Takeshi Suzuki and María José Boccardi August 13, 2013 Dynamical Systems are systems, described by one or more equations, that evolve over time.

More information

ORDINARY DIFFERENTIAL EQUATIONS

ORDINARY DIFFERENTIAL EQUATIONS ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 4884 NOVEMBER 9, 7 Summary This is an introduction to ordinary differential equations We

More information

! 4 4! o! +! h 4 o=0! ±= ± p i And back-substituting into the linear equations gave us the ratios of the amplitudes of oscillation:.»» = A p e i! +t»»

! 4 4! o! +! h 4 o=0! ±= ± p i And back-substituting into the linear equations gave us the ratios of the amplitudes of oscillation:.»» = A p e i! +t»» Topic 6: Coupled Oscillators and Normal Modes Reading assignment: Hand and Finch Chapter 9 We are going to be considering the general case of a system with N degrees of freedome close to one of its stable

More information

3.3. SYSTEMS OF ODES 1. y 0 " 2y" y 0 + 2y = x1. x2 x3. x = y(t) = c 1 e t + c 2 e t + c 3 e 2t. _x = A x + f; x(0) = x 0.

3.3. SYSTEMS OF ODES 1. y 0  2y y 0 + 2y = x1. x2 x3. x = y(t) = c 1 e t + c 2 e t + c 3 e 2t. _x = A x + f; x(0) = x 0. .. SYSTEMS OF ODES. Systems of ODEs MATH 94 FALL 98 PRELIM # 94FA8PQ.tex.. a) Convert the third order dierential equation into a rst oder system _x = A x, with y " y" y + y = x = @ x x x b) The equation

More information

Math Ordinary Differential Equations

Math Ordinary Differential Equations Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x

More information

SYSTEMS OF ODES. mg sin ( (x)) dx 2 =

SYSTEMS OF ODES. mg sin ( (x)) dx 2 = SYSTEMS OF ODES Consider the pendulum shown below. Assume the rod is of neglible mass, that the pendulum is of mass m, and that the rod is of length `. Assume the pendulum moves in the plane shown, and

More information

The Plan. Initial value problems (ivps) for ordinary differential equations (odes) Review of basic methods You are here: Hamiltonian systems

The Plan. Initial value problems (ivps) for ordinary differential equations (odes) Review of basic methods You are here: Hamiltonian systems Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit V Solution of Differential Equations Part 2 Dianne P. O Leary c 2008 The Plan

More information

1 The pendulum equation

1 The pendulum equation Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating

More information

NUMERICAL SOLUTION OF ODE IVPs. Overview

NUMERICAL SOLUTION OF ODE IVPs. Overview NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental

More information

ORDINARY DIFFERENTIAL EQUATIONS

ORDINARY DIFFERENTIAL EQUATIONS ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 48824. JANUARY 3, 25 Summary. This is an introduction to ordinary differential equations.

More information

Local Stability Analysis of a Mathematical Model of the Interaction of Two Populations of Differential Equations (Host-Parasitoid)

Local Stability Analysis of a Mathematical Model of the Interaction of Two Populations of Differential Equations (Host-Parasitoid) Biology Medicine & Natural Product Chemistry ISSN: 089-6514 Volume 5 Number 1 016 Pages: 9-14 DOI: 10.1441/biomedich.016.51.9-14 Local Stability Analysis of a Mathematical Model of the Interaction of Two

More information

3 Stability and Lyapunov Functions

3 Stability and Lyapunov Functions CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such

More information

Lab #2 - Two Degrees-of-Freedom Oscillator

Lab #2 - Two Degrees-of-Freedom Oscillator Lab #2 - Two Degrees-of-Freedom Oscillator Last Updated: March 0, 2007 INTRODUCTION The system illustrated in Figure has two degrees-of-freedom. This means that two is the minimum number of coordinates

More information

Numerical Integration of Equations of Motion

Numerical Integration of Equations of Motion GraSMech course 2009-2010 Computer-aided analysis of rigid and flexible multibody systems Numerical Integration of Equations of Motion Prof. Olivier Verlinden (FPMs) Olivier.Verlinden@fpms.ac.be Prof.

More information

Entrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.

Entrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4. Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the

More information

A brief introduction to ordinary differential equations

A brief introduction to ordinary differential equations Chapter 1 A brief introduction to ordinary differential equations 1.1 Introduction An ordinary differential equation (ode) is an equation that relates a function of one variable, y(t), with its derivative(s)

More information

Systems of Linear ODEs

Systems of Linear ODEs P a g e 1 Systems of Linear ODEs Systems of ordinary differential equations can be solved in much the same way as discrete dynamical systems if the differential equations are linear. We will focus here

More information

Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)

Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas) Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in

More information

Coupled differential equations

Coupled differential equations Coupled differential equations Example: dy 1 dy 11 1 1 1 1 1 a y a y b x a y a y b x Consider the case with b1 b d y1 a11 a1 y1 y a1 ay dy y y e y 4/1/18 One way to address this sort of problem, is to

More information

Chapter 13 Internal (Lyapunov) Stability 13.1 Introduction We have already seen some examples of both stable and unstable systems. The objective of th

Chapter 13 Internal (Lyapunov) Stability 13.1 Introduction We have already seen some examples of both stable and unstable systems. The objective of th Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology 1 1 c Chapter

More information

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) 28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.

More information

2nd-Order Linear Equations

2nd-Order Linear Equations 4 2nd-Order Linear Equations 4.1 Linear Independence of Functions In linear algebra the notion of linear independence arises frequently in the context of vector spaces. If V is a vector space over the

More information

Math 232, Final Test, 20 March 2007

Math 232, Final Test, 20 March 2007 Math 232, Final Test, 20 March 2007 Name: Instructions. Do any five of the first six questions, and any five of the last six questions. Please do your best, and show all appropriate details in your solutions.

More information

Ordinary Differential Equations

Ordinary Differential Equations Chapter 13 Ordinary Differential Equations We motivated the problem of interpolation in Chapter 11 by transitioning from analzying to finding functions. That is, in problems like interpolation and regression,

More information

Ordinary Differential Equation Theory

Ordinary Differential Equation Theory Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit

More information

7 Planar systems of linear ODE

7 Planar systems of linear ODE 7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution

More information

Higher Order Averaging : periodic solutions, linear systems and an application

Higher Order Averaging : periodic solutions, linear systems and an application Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,

More information

Chapter 4. Systems of ODEs. Phase Plane. Qualitative Methods

Chapter 4. Systems of ODEs. Phase Plane. Qualitative Methods Chapter 4 Systems of ODEs. Phase Plane. Qualitative Methods Contents 4.0 Basics of Matrices and Vectors 4.1 Systems of ODEs as Models 4.2 Basic Theory of Systems of ODEs 4.3 Constant-Coefficient Systems.

More information

MCE693/793: Analysis and Control of Nonlinear Systems

MCE693/793: Analysis and Control of Nonlinear Systems MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear

More information

Coordinate Curves for Trajectories

Coordinate Curves for Trajectories 43 The material on linearizations and Jacobian matrices developed in the last chapter certainly expanded our ability to deal with nonlinear systems of differential equations Unfortunately, those tools

More information

N.G.Bean, D.A.Green and P.G.Taylor. University of Adelaide. Adelaide. Abstract. process of an MMPP/M/1 queue is not a MAP unless the queue is a

N.G.Bean, D.A.Green and P.G.Taylor. University of Adelaide. Adelaide. Abstract. process of an MMPP/M/1 queue is not a MAP unless the queue is a WHEN IS A MAP POISSON N.G.Bean, D.A.Green and P.G.Taylor Department of Applied Mathematics University of Adelaide Adelaide 55 Abstract In a recent paper, Olivier and Walrand (994) claimed that the departure

More information

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.01 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential

More information

2D-Volterra-Lotka Modeling For 2 Species

2D-Volterra-Lotka Modeling For 2 Species Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose

More information

Lecture 2: Review of Prerequisites. Table of contents

Lecture 2: Review of Prerequisites. Table of contents Math 348 Fall 217 Lecture 2: Review of Prerequisites Disclaimer. As we have a textbook, this lecture note is for guidance and supplement only. It should not be relied on when preparing for exams. In this

More information

21 Linear State-Space Representations

21 Linear State-Space Representations ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may

More information

STABILITY. Phase portraits and local stability

STABILITY. Phase portraits and local stability MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),

More information

Autonomous Systems and Stability

Autonomous Systems and Stability LECTURE 8 Autonomous Systems and Stability An autonomous system is a system of ordinary differential equations of the form 1 1 ( 1 ) 2 2 ( 1 ). ( 1 ) or, in vector notation, x 0 F (x) That is to say, an

More information

Interactions. Yuan Gao. Spring Applied Mathematics University of Washington

Interactions. Yuan Gao. Spring Applied Mathematics University of Washington Interactions Yuan Gao Applied Mathematics University of Washington yuangao@uw.edu Spring 2015 1 / 27 Nonlinear System Consider the following coupled ODEs: dx = f (x, y). dt dy = g(x, y). dt In general,

More information

Solutions to Math 53 Math 53 Practice Final

Solutions to Math 53 Math 53 Practice Final Solutions to Math 5 Math 5 Practice Final 20 points Consider the initial value problem y t 4yt = te t with y 0 = and y0 = 0 a 8 points Find the Laplace transform of the solution of this IVP b 8 points

More information

Ordinary Differential Equations

Ordinary Differential Equations Ordinary Differential Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro September 19, 2014 1 / 55 Motivation

More information

Lecture 6. Eigen-analysis

Lecture 6. Eigen-analysis Lecture 6 Eigen-analysis University of British Columbia, Vancouver Yue-Xian Li March 7 6 Definition of eigenvectors and eigenvalues Def: Any n n matrix A defines a LT, A : R n R n A vector v and a scalar

More information

(Refer Slide Time: 00:32)

(Refer Slide Time: 00:32) Nonlinear Dynamical Systems Prof. Madhu. N. Belur and Prof. Harish. K. Pillai Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture - 12 Scilab simulation of Lotka Volterra

More information

Solution of Linear State-space Systems

Solution of Linear State-space Systems Solution of Linear State-space Systems Homogeneous (u=0) LTV systems first Theorem (Peano-Baker series) The unique solution to x(t) = (t, )x 0 where The matrix function is given by is called the state

More information

Mathematics for chemical engineers. Numerical solution of ordinary differential equations

Mathematics for chemical engineers. Numerical solution of ordinary differential equations Mathematics for chemical engineers Drahoslava Janovská Numerical solution of ordinary differential equations Initial value problem Winter Semester 2015-2016 Outline 1 Introduction 2 One step methods Euler

More information

Lecture 20/Lab 21: Systems of Nonlinear ODEs

Lecture 20/Lab 21: Systems of Nonlinear ODEs Lecture 20/Lab 21: Systems of Nonlinear ODEs MAR514 Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth Coupled ODEs: Species

More information

Linear Algebra Review (Course Notes for Math 308H - Spring 2016)

Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Dr. Michael S. Pilant February 12, 2016 1 Background: We begin with one of the most fundamental notions in R 2, distance. Letting (x 1,

More information

A note on the uniform perturbation index 1

A note on the uniform perturbation index 1 Rostock. Math. Kolloq. 52, 33 46 (998) Subject Classification (AMS) 65L5 M. Arnold A note on the uniform perturbation index ABSTRACT. For a given differential-algebraic equation (DAE) the perturbation

More information

LAURENTIAN UNIVERSITY UNIVERSITÉ LAURENTIENNE

LAURENTIAN UNIVERSITY UNIVERSITÉ LAURENTIENNE Page 1 of 15 LAURENTIAN UNIVERSITY UNIVERSITÉ LAURENTIENNE Friday, December 14 th Course and No. 2007 MATH 2066 EL Date................................... Cours et no........................... Total no.

More information

MIDTERM REVIEW AND SAMPLE EXAM. Contents

MIDTERM REVIEW AND SAMPLE EXAM. Contents MIDTERM REVIEW AND SAMPLE EXAM Abstract These notes outline the material for the upcoming exam Note that the review is divided into the two main topics we have covered thus far, namely, ordinary differential

More information

Math 215/255 Final Exam (Dec 2005)

Math 215/255 Final Exam (Dec 2005) Exam (Dec 2005) Last Student #: First name: Signature: Circle your section #: Burggraf=0, Peterson=02, Khadra=03, Burghelea=04, Li=05 I have read and understood the instructions below: Please sign: Instructions:.

More information

Math 2a Prac Lectures on Differential Equations

Math 2a Prac Lectures on Differential Equations Math 2a Prac Lectures on Differential Equations Prof. Dinakar Ramakrishnan 272 Sloan, 253-37 Caltech Office Hours: Fridays 4 5 PM Based on notes taken in class by Stephanie Laga, with a few added comments

More information

The Exponential of a Matrix

The Exponential of a Matrix The Exponential of a Matrix 5-8- The solution to the exponential growth equation dx dt kx is given by x c e kt It is natural to ask whether you can solve a constant coefficient linear system x A x in a

More information

PH.D. PRELIMINARY EXAMINATION MATHEMATICS

PH.D. PRELIMINARY EXAMINATION MATHEMATICS UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem

More information

A(t) L T *(t) t/t t/t

A(t) L T *(t) t/t t/t he Real Story behind Floquet-Lyapunov Factorizations RAYMOND J. SPIERI Department of Mathematics and Statistics Acadia University Wolfville, NS BP X CANADA raymond.spiteri@acadiau.ca http//ace.acadiau.ca/~rspiteri

More information

Econ Lecture 14. Outline

Econ Lecture 14. Outline Econ 204 2010 Lecture 14 Outline 1. Differential Equations and Solutions 2. Existence and Uniqueness of Solutions 3. Autonomous Differential Equations 4. Complex Exponentials 5. Linear Differential Equations

More information

PH.D. PRELIMINARY EXAMINATION MATHEMATICS

PH.D. PRELIMINARY EXAMINATION MATHEMATICS UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem

More information

Our goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always

Our goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always October 5 Relevant reading: Section 2.1, 2.2, 2.3 and 2.4 Our goal is to solve a general constant coecient linear second order ODE a d2 y dt + bdy + cy = g (t) 2 dt where a, b, c are constants and a 0.

More information

Lecture Notes in Mathematics. Arkansas Tech University Department of Mathematics

Lecture Notes in Mathematics. Arkansas Tech University Department of Mathematics Lecture Notes in Mathematics Arkansas Tech University Department of Mathematics Introductory Notes in Ordinary Differential Equations for Physical Sciences and Engineering Marcel B. Finan c All Rights

More information

The Fundamental Theorem of Calculus: Suppose f continuous on [a, b]. 1.) If G(x) = x. f(t)dt = F (b) F (a) where F is any antiderivative

The Fundamental Theorem of Calculus: Suppose f continuous on [a, b]. 1.) If G(x) = x. f(t)dt = F (b) F (a) where F is any antiderivative 1 Calulus pre-requisites you must know. Derivative = slope of tangent line = rate. Integral = area between curve and x-axis (where area can be negative). The Fundamental Theorem of Calculus: Suppose f

More information

The goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T

The goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T 1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.

More information

Second-Order Linear ODEs

Second-Order Linear ODEs Second-Order Linear ODEs A second order ODE is called linear if it can be written as y + p(t)y + q(t)y = r(t). (0.1) It is called homogeneous if r(t) = 0, and nonhomogeneous otherwise. We shall assume

More information

LMI Methods in Optimal and Robust Control

LMI Methods in Optimal and Robust Control LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear

More information

Understand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.

Understand the existence and uniqueness theorems and what they tell you about solutions to initial value problems. Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics

More information

Differential Equations and Linear Algebra Supplementary Notes. Simon J.A. Malham. Department of Mathematics, Heriot-Watt University

Differential Equations and Linear Algebra Supplementary Notes. Simon J.A. Malham. Department of Mathematics, Heriot-Watt University Differential Equations and Linear Algebra Supplementary Notes Simon J.A. Malham Department of Mathematics, Heriot-Watt University Contents Chapter 1. Linear algebraic equations 5 1.1. Gaussian elimination

More information

Solutions to the Review Questions

Solutions to the Review Questions Solutions to the Review Questions Short Answer/True or False. True or False, and explain: (a) If y = y + 2t, then 0 = y + 2t is an equilibrium solution. False: (a) Equilibrium solutions are only defined

More information

Lecture I Introduction, concept of solutions, application

Lecture I Introduction, concept of solutions, application S. Ghorai Lecture I Introduction, concept of solutions, application Definition. A differential equation (DE) is a relation that contains a finite set of functions and their derivatives with respect to

More information

Final Exam Review. Review of Systems of ODE. Differential Equations Lia Vas. 1. Find all the equilibrium points of the following systems.

Final Exam Review. Review of Systems of ODE. Differential Equations Lia Vas. 1. Find all the equilibrium points of the following systems. Differential Equations Lia Vas Review of Systems of ODE Final Exam Review 1. Find all the equilibrium points of the following systems. (a) dx = x x xy (b) dx = x x xy = 0.5y y 0.5xy = 0.5y 0.5y 0.5xy.

More information

Econ 204 Differential Equations. 1 Existence and Uniqueness of Solutions

Econ 204 Differential Equations. 1 Existence and Uniqueness of Solutions Econ 4 Differential Equations In this supplement, we use the methods we have developed so far to study differential equations. 1 Existence and Uniqueness of Solutions Definition 1 A differential equation

More information

Physics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics

Physics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics Applications of nonlinear ODE systems: Physics: spring-mass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:

More information

MATH 353 LECTURE NOTES: WEEK 1 FIRST ORDER ODES

MATH 353 LECTURE NOTES: WEEK 1 FIRST ORDER ODES MATH 353 LECTURE NOTES: WEEK 1 FIRST ORDER ODES J. WONG (FALL 2017) What did we cover this week? Basic definitions: DEs, linear operators, homogeneous (linear) ODEs. Solution techniques for some classes

More information

AMATH 351 Mar 15, 2013 FINAL REVIEW. Instructor: Jiri Najemnik

AMATH 351 Mar 15, 2013 FINAL REVIEW. Instructor: Jiri Najemnik AMATH 351 Mar 15, 013 FINAL REVIEW Instructor: Jiri Najemni ABOUT GRADES Scores I have so far will be posted on the website today sorted by the student number HW4 & Exam will be added early next wee Let

More information

Sufficiency Conditions for Finite Escape Times in Systems of Quadratic Differential Equations

Sufficiency Conditions for Finite Escape Times in Systems of Quadratic Differential Equations /. Inst. Maths Applies (1977) 19, 377-383 Sufficiency Conditions for Finite Escape Times in Systems of Quadratic Differential Equations W. M. GETZ AND D. H. JACOBSON National Research Institute for Mathematical

More information

Ordinary Differential Equations

Ordinary Differential Equations Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005 Overview What are Ordinary Differential Equations

More information

MAT 22B - Lecture Notes

MAT 22B - Lecture Notes MAT 22B - Lecture Notes 4 September 205 Solving Systems of ODE Last time we talked a bit about how systems of ODE arise and why they are nice for visualization. Now we'll talk about the basics of how to

More information

Notes for Numerical Analysis Math 5466 by S. Adjerid Virginia Polytechnic Institute and State University (A Rough Draft) Contents Numerical Methods for ODEs 5. Introduction............................

More information

4. Higher Order Linear DEs

4. Higher Order Linear DEs 4. Higher Order Linear DEs Department of Mathematics & Statistics ASU Outline of Chapter 4 1 General Theory of nth Order Linear Equations 2 Homogeneous Equations with Constant Coecients 3 The Method of

More information

y0 = F (t0)+c implies C = y0 F (t0) Integral = area between curve and x-axis (where I.e., f(t)dt = F (b) F (a) wheref is any antiderivative 2.

y0 = F (t0)+c implies C = y0 F (t0) Integral = area between curve and x-axis (where I.e., f(t)dt = F (b) F (a) wheref is any antiderivative 2. Calulus pre-requisites you must know. Derivative = slope of tangent line = rate. Integral = area between curve and x-axis (where area can be negative). The Fundamental Theorem of Calculus: Suppose f continuous

More information

On The Numerical Solution of Differential-Algebraic Equations(DAES) with Index-3 by Pade Approximation

On The Numerical Solution of Differential-Algebraic Equations(DAES) with Index-3 by Pade Approximation Appl. Math. Inf. Sci. Lett., No. 2, 7-23 (203) 7 Applied Mathematics & Information Sciences Letters An International Journal On The Numerical Solution of Differential-Algebraic Equations(DAES) with Index-3

More information

THETA-LOGISTIC PREDATOR PREY

THETA-LOGISTIC PREDATOR PREY THETA-LOGISTIC PREDATOR PREY What are the assumptions of this model? 1.) Functional responses are non-linear. Functional response refers to a change in the rate of exploitation of prey by an individual

More information

Lecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems

Lecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems Lecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems Shawn D. Ryan Spring 2012 Last Time: We finished Chapter 9: Nonlinear Differential Equations and Stability. Now

More information

Introduction to Ordinary Differential Equations with Mathematica

Introduction to Ordinary Differential Equations with Mathematica ALFRED GRAY MICHAEL MEZZINO MARKA. PINSKY Introduction to Ordinary Differential Equations with Mathematica An Integrated Multimedia Approach %JmT} Web-Enhanced Includes CD-ROM TABLE OF CONTENTS Preface

More information

Boyce/DiPrima/Meade 11 th ed, Ch 1.1: Basic Mathematical Models; Direction Fields

Boyce/DiPrima/Meade 11 th ed, Ch 1.1: Basic Mathematical Models; Direction Fields Boyce/DiPrima/Meade 11 th ed, Ch 1.1: Basic Mathematical Models; Direction Fields Elementary Differential Equations and Boundary Value Problems, 11 th edition, by William E. Boyce, Richard C. DiPrima,

More information

The Unifying View on Ordinary Differential Equations and Automatic Differentiation, yet with a Gap to Fill Alexander Gofen.

The Unifying View on Ordinary Differential Equations and Automatic Differentiation, yet with a Gap to Fill Alexander Gofen. The Unifying View on Ordinary Differential Equations and Automatic Differentiation, yet with a Gap to Fill Alexander Gofen Conventions Differentiation is understood as for holomorphic functions in the

More information

Fundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad

Fundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian

More information

Lecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9

Lecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9 Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction

More information

Problem set 7 Math 207A, Fall 2011 Solutions

Problem set 7 Math 207A, Fall 2011 Solutions Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase

More information

Topic 6: Projected Dynamical Systems

Topic 6: Projected Dynamical Systems Topic 6: Projected Dynamical Systems John F. Smith Memorial Professor and Director Virtual Center for Supernetworks Isenberg School of Management University of Massachusetts Amherst, Massachusetts 01003

More information

Nonlinear System Analysis

Nonlinear System Analysis Nonlinear System Analysis Lyapunov Based Approach Lecture 4 Module 1 Dr. Laxmidhar Behera Department of Electrical Engineering, Indian Institute of Technology, Kanpur. January 4, 2003 Intelligent Control

More information

Lecture 19: Solving linear ODEs + separable techniques for nonlinear ODE s

Lecture 19: Solving linear ODEs + separable techniques for nonlinear ODE s Lecture 19: Solving linear ODEs + separable techniques for nonlinear ODE s Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth

More information

Non-Dimensionalization of Differential Equations. September 9, 2014

Non-Dimensionalization of Differential Equations. September 9, 2014 Non-Dimensionalization of Differential Equations September 9, 2014 1 Non-dimensionalization The most elementary idea of dimensional analysis is this: in any equation, the dimensions of the right hand side

More information

Scientific Computing: An Introductory Survey

Scientific Computing: An Introductory Survey Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign

More information

Math Assignment 2

Math Assignment 2 Math 2280 - Assignment 2 Dylan Zwick Spring 2014 Section 1.5-1, 15, 21, 29, 38, 42 Section 1.6-1, 3, 13, 16, 22, 26, 31, 36, 56 Section 2.1-1, 8, 11, 16, 29 Section 2.2-1, 10, 21, 23, 24 1 Section 1.5

More information

SMA 208: Ordinary differential equations I

SMA 208: Ordinary differential equations I SMA 208: Ordinary differential equations I First Order differential equations Lecturer: Dr. Philip Ngare (Contacts: pngare@uonbi.ac.ke, Tue 12-2 PM) School of Mathematics, University of Nairobi Feb 26,

More information

Differential Equations and Modeling

Differential Equations and Modeling Differential Equations and Modeling Preliminary Lecture Notes Adolfo J. Rumbos c Draft date: March 22, 2018 March 22, 2018 2 Contents 1 Preface 5 2 Introduction to Modeling 7 2.1 Constructing Models.........................

More information

Applications of eigenvalues/eigenvectors

Applications of eigenvalues/eigenvectors Example 1. Applications of eigenvalues/eigenvectors Predator-Prey models. Foxes and Rabbits share a large forest, with the foxes eating rabbits and rabbits eating the abundant vegetation. The sizes of

More information

Initial value problems for ordinary differential equations

Initial value problems for ordinary differential equations Initial value problems for ordinary differential equations Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 IVP

More information

M A : Ordinary Differential Equations

M A : Ordinary Differential Equations M A 2 0 5 1: Ordinary Differential Equations Essential Class Notes & Graphics D 19 * 2018-2019 Sections D07 D11 & D14 1 1. INTRODUCTION CLASS 1 ODE: Course s Overarching Functions An introduction to the

More information