ENGI Linear Approximation (2) Page Linear Approximation to a System of Non-Linear ODEs (2)

Size: px
Start display at page:

Download "ENGI Linear Approximation (2) Page Linear Approximation to a System of Non-Linear ODEs (2)"

Transcription

1 ENGI Linear Approximation () Page Linear Approximation to a System of Non-Linear ODEs () From sections 4.0 and 4.0, the non-linear system dx dy = x = P( x, y), = y = Q( x, y) () with critical point at (0, 0) may be expressed as x = ax + by + P ( x, y) y = cx + d y + Q x, y () where a, b, c, d are all constants and P( x, y) Q( x, y) lim = 0 and lim ( xy, ) ( 0,0 ) ( xy, ) ( 0,0) x + y x + y = 0. Near the critical point (0, 0), this system may be approximated by the linear system x = ax + by y = cx + d y () Effect of Small Perturbations Small perturbations in the values of the coefficients a, b, c, d are reflected in small changes in the eigenvalues λ. If the eigenvalues are a pure imaginary pair, λ = ± jv, then the critical point is a centre. The effect of small changes in the coefficients will change the eigenvalues to the complex conjugate pair λ' = u' ± jv', where u' is small in magnitude and v' is close to v. The trajectories of the new system are likely to be spirals. The critical point will be an asymptotically stable focus if u' < 0, a stable (but not asymptotically stable) centre if u' = 0, but it will be an unstable focus if u' > 0. Therefore small perturbations in a linear system with pure imaginary eigenvalues are likely to result in radical changes in the trajectories (orbits) and may change the stable system into an unstable system.

2 ENGI Linear Approximation () Page 4. If the eigenvalues are a real equal pair, λ = λ, then a slight perturbation is likely to separate the roots into distinct values. If those values are still real, then the critical point remains a node. If the perturbed eigenvalues are a complex conjugate pair, then the nature of the trajectories will change into spirals and the critical point changes from a node into a focus. However, in both cases, an asymptotically stable critical point remains asymptotically stable after a small perturbation, while an unstable critical point remains unstable. In all other cases, a slight perturbation leaves the sign of the real part of both eigenvalues unchanged and affects neither the type of critical point nor the overall type of the orbits. These results are summarized in the two following theorems. Poincaré s Theorem: The singularities of the non-linear system () are identical to the singularities of the linear system (), except for the cases D < 0 and a + d = 0, which is a centre in the linear case, but may be a centre or a focus in the non-linear case; and D = 0, which is either a node or a focus in the non-linear case. Theorem on stability of the singularity at (0, 0): Linear approximation asymptotically stable unstable stable but not asymptotically stable Non-linear system asymptotically stable unstable any (unstable, stable or asymptotically stable)

3 ENGI Linear Approximation () Page 4.4 Example Perform a stability analysis on the system dx dy = x x + xy, = y xy 6y Find the critical points: dx dy x( x+ y) = 0 and y( x 6y) = 0 This generates four solutions: x = 0 and y = 0 x, y = 0,0 = = 0 x x + xy = 0 and y xy y = x = 0 and x 6y = 0 x, y = 0, x+ y = 0 and y = 0 x, y =,0 8 x+ y = 0 and x 6y = 0 x, y =, Linearize the system near each critical point Near the critical point (0, 0), it is obvious that the linear approximation to the system is dx dy = x, = y This is an uncoupled system of ODEs, whose general solution is quickly found (by direct integration or by the results on page 4.0) to be t t x t, y t = ce, c e The critical point (0, 0) is an unstable node. Also note that y = k x, so that the solutions are parabolas all sharing the same vertex at the origin, except for the lines x = 0 (from c = 0) and y = 0 (from c = 0). The orbits all begin at the origin and move away, (a hallmark of an unstable node). ( () ()) ( ) This diagram is valid for the linear approximation everywhere, but is valid only in the immediate neighbourhood of (0, 0) for the non-linear system.

4 ENGI Linear Approximation () Page 4.5 Example (continued) P x, y = x x + xy, Q x, y = y xy 6 y P P = x + y = x x y Q Q and = y = x y x y Near a critical point (a, b), the linear approximation to the system is P P x x ( ab, ) y ( ab, ) x a a+ b a x a = = y Q Q y b b a b y b x ( ab, ) y ( ab, ) Near the critical point ( 0, ), the linear approximation to the system is 4 x + 0 x 0 x = = y y y xz, = xy, leads to A simple change of variables to 4 0 x x = z z The matrix is triangular, allowing us to identify the eigenvalues immediately. The eigenvalues are of opposite sign (+4/ and ). Therefore the critical point at ( 0, ) is an unstable saddle point. From page 4.0, 4 4 ( 0 ) 0 D = a d + bc = + + = The eigenvector corresponding to eigenvalue λ = is 0 0 α a d D 0 = any non-zero multiple of β = =, c The eigenvector corresponding to eigenvalue λ = +4/ is 0 0 α + a d + D 0 = any non-zero multiple of β = = c

5 ENGI Linear Approximation () Page 4.6 Example (continued) The asymptotes are therefore x = 0 (corresponding to c = 0; inward because λ < 0 ) and x = 0( y ), (corresponding to c = 0; outward because λ > 0 ). This is sufficient to sketch the phase portrait, even without the general solution. The general solution near the critical point ( 0, ) is 4 4 ( (), ()) = ( (), () ) = 0, + t t t x t z t x t y t c e c e c e where c, c are arbitrary constants. All trajectories (except for c = 0) therefore come in from infinity near the asymptote x = 0 (where c = 0) and all trajectories (except for c = 0) return to infinity near the asymptote y = 0 x (where c = 0). The diagram is valid for the linear approximation everywhere, but is valid only in the immediate neighbourhood of ( 0, ) for the non-linear system.

6 ENGI Linear Approximation () Page 4.7 Example (continued) Near the critical point (, 0), the linear approximation to the system is x x x = = y 0 y 0 y Change variables to ( wy, ) = ( x, y). Then w w = y 0 y The matrix is triangular, with eigenvalues of opposite sign ( and +). Therefore the critical point at (, 0) is also a[n unstable] saddle point. D = a d + 4bc = + 0 = 4 a d D α = any non-zero multiple of β = =, 0 c 0 c = 0, so the alternative form is needed for the other eigenvector: b α = any non-zero multiple of β ( d a) + D = + = The general solution near the critical point (, 0) is ( () ()) () () t t t ( ) w t, y t = x t, y t = c e + c e, c e where c, c are arbitrary constants. All trajectories (except for c = 0) therefore come in from infinity near the asymptote y = 0 (where c = 0) and all trajectories (except for c = 0) return to infinity near the asymptote y = (x ) (where c = 0). Again, the diagram is valid for the linear approximation everywhere, but is valid only in the immediate neighbourhood of (, 0) for the non-linear system.

7 ENGI Linear Approximation () Page 4.8 Example (continued) Near the critical point ( 8 7, 7 ), the linear approximation to the system is x x x = = y y y Change variables to ( wz, ) = ( x 8, y ) w 8 8 w = z 7 6 z D = ( a d) + 4bc = + 4 = = D < 0 and (a + d) < 0 the critical point is an asymptotically stable focus λ 7 The eigenvalues, which are the solutions to det = 0, are λ ( a+ d) ± D 4 ± 8 7 λ = = = ± j Using the formula on page 4.0, u =, v= =, u d = + = ( wt, zt ) xt, yt = = 7 7 t t t t t t t e c cos sin c4 cos sin, ccos c sin Redefining the arbitrary constants as c5 = c and c6 = c 4, 7 7 the general solution near the critical point 8 ( 7 7) ( wt (), zt ()) = 8 ( xt (), yt () ) = 7 7, is t t t t t t e c5 cos + 7 sin + c6 7 cos + sin, c5cos + c6 sin where c 5, c 6 are arbitrary constants. t 7 All trajectories spiral in to the critical point [a phase portrait for the linear approximation is on the next page].

8 ENGI Linear Approximation () Page 4.9 Example (continued) Maple produces the following direction field plots: One can clearly see trajectories flowing outward from the unstable node at the origin in all directions. The natures of the other critical points are somewhat less obvious. Zooming in to the neighbourhood of one of the saddle points and to the neighbourhood of the stable focus produces the next pair of diagrams:

9 ENGI Linear Approximation () Page 4.40 Example (continued) Maple can also superimpose some trajectories on these phase portraits: Maple commands for this plot: > with(detools): > phaseportrait( [diff(x(t),t) = x(t) - x(t)^ + x(t)*y(t), diff(y(t),t) = *y(t) - x(t)*y(t) - 6*y(t)^], [x(t), y(t)], t=-0..0, [[x(0)=0., y(0)=0.0], [x(0)=0.8, y(0)=-0.0], [x(0)=., y(0)=0.0], [x(0)=0.05, y(0)=0.], [x(0)=0.05, y(0)=0.4], [x(0)=., y(0)=-0.0]], x=-0...4, y= , stepsize=.0, colour=red, linecolour=[blue, cyan, magenta, sienna, orange, black], title=`example Non-Linear Solution`);

10 ENGI Linear Approximation () Page 4.4 Example Perform a stability analysis on the system dx dy = y + x x y, = x + y x y Find the critical points: dx dy Clearly (x, y) = (0, 0) satisfies = = 0. dx dy = = x = 0 y = 0 dx dy = = y = 0 x= 0 If x 0 and y 0, then at any critical point ( ) ( ) 0 ( ) and ( ) y + x x y = x + y x y = y = x x y x = y x y x y y ( x y ) = = = y x x which has no solution for real (x, y). Therefore the only critical point is (0, 0). The linear approximation to the non-linear system is x x = y y D = a d + 4bc = + 4 = 4 < 0 (a + d) = > 0 the critical point is an unstable focus. Using the formula on page 4.0, a+ d a d D u = =, u d = = 0, v = = x t, y t = e t c sin t c cos t, c cost + c sin t ( () ()) ( 4 4 ) and c, c 4 are [real] arbitrary constants.

11 ENGI Linear Approximation () Page 4.4 Example (continued) Solution in the neighbourhood of the only critical point (0, 0): Now consider the distance r of any point (x, y) from the critical point (0, 0): dr dx dy r = x + y r = x + y [chain rule] From the original non-linear system: dr r = x( y + x( x y )) + y( x + y( x y )) = xy + x r xy + y r = r r ( r ) ( r ) dr < 0 > = r( r ) > 0 < Therefore solutions starting closer than one unit to the critical point spiral out, but solutions starting further away than one unit approach the critical point. A solution on the circle r = never changes its distance from the origin and stays on that circle, but is not stationary (because the only critical point is at (0, 0), not on that circle). Note that x + y = is a solution to the non-linear equation: dx dy = y + x( x y ) = y and = x + y( x y ) = x dy dy dx x = = dx y dy dy x But x + y = x + y = 0 = dx dx y No other solution may cross this solution. Therefore all non-trivial solutions approach the unit circle in the limit as t. The unit circle here is an example of a limit cycle.

12 ENGI Linear Approximation () Page 4.4 Example (continued) Further notes: Let u = r, then du = r dr and r dr = r ( r ) du = u( u) du = u u ( ) + du = ln u ln ( u) = t + C u u u u t + C t ln = t + C = e = ce u u u ce t uce t ce t t = + u = ce t ce + u =, where c = c u = r() t = t + ce + ce t If r(0) = r o, then t u ro ro ce = c = c = u ro ro Therefore r() t = r o + e ro t

13 ENGI Linear Approximation () Page 4.44 Example (continued) Also note that, from the polar coordinate system, ( xy, ) = ( rcos θ, r sinθ ) xy, = r cosθ rsin θθ, r sinθ + rcosθθ ( θ ) θ ( θ) θ θ ( θ) θ ( θ) xy yx = rcos r sin + rcos rcos rsin r cos + rsin rsinθ θ ( cos sin ) = r θ + θ θ = r θ But the non-linear system is x, y = y + x r, x + y r ( ) ( ) xy yx = x + xy r y xy r = x + y = r Therefore = = = + C r θ r θ θ t All paths move clockwise with constant angular speed. In Cartesian coordinates all orbits can therefore be described by x() t, y() t = cos t θo, sin t θo r o t + e ro ( ) where (x(0), y(0)) = r o (cos θ o, sin θ o )

14 ENGI Linear Approximation () Page 4.45 Example (A more challenging and tedious case, for reference only) Find and determine the nature of all critical points of the system dx x dy = e y, = y sinx () dx dy x = = 0 y = e and y = sinx x Critical points occur where the graphs of y = e and y = sin x intersect. The critical points are (0, 0) and (x i, y i ), where ( 4 π ) δ, ( 4 π ) ε ; (,,, ) x = i i + i = i i i Linearize for the critical point at (0, 0): near ( 0,0) ( 0,0) ( 0,0) 0 x P P e y x + y = e x + y x y Q Q y sin x x + y = ( cos 0) x + y x y near ( 0,0) ( 0,0) ( 0,0) Therefore the linear system that models the non-linear system () near (0, 0) is x x = y y Finding the eigenvalues: λ det ( A λi) = = ( + λ)( λ) ( )( ) = λ λ ( A I) det λ = 0 λ = λ = ± () The eigenvalues are real and of opposite sign. The critical point (0, 0) of () [and therefore also of ()] is an unstable saddle point.

15 ENGI Linear Approximation () Page 4.46 Example (continued) Using the results on page 4.0, D = ( a d) + 4bc = ( ) + 4( )( ) = 4+ 4 = 8 The eigenvectors are α ( a d) D = any non-zero multiple of β = c for λ = and α ( a d) + D + = any non-zero multiple of β = c for λ = + Choose a multiple of in both cases. The general solution of () is therefore () () () () ( t t t ) ( () ()) t xt, yt = c + e + c e, ce + ce y t lim = = = > 0 and t x t + y t + lim = = = < 0 t + x t The trajectories (except for c = 0) therefore come in from infinity along the asymptote y = x. The trajectories (except for c = 0) return to infinity along the asymptote y = + x. The phase space diagram for this solution (completely valid only for the linear system ()) is The phase space for the non-linear system () resembles this diagram only in the immediate neighbourhood of the critical point (0, 0).

16 ENGI Linear Approximation () Page 4.47 Example (continued) At other critical points (k, l), dx dy = = 0 l = sink e sink = 0 Linearizing (Taylor s series for P(x, y) about (x, y) = (k, l): x P P ( e y ) ( x k) + ( y l) = ( e )( x k) + ( )( y l near ( kl, ) x y near ( kl, ) ( kl) ( kl) ( kl),, Q Q y sin x ( x k) + y l = cos k x k + x y ( kl),, ( y l) Therefore the linear system that models the non-linear system () near (k, l) is x e x k = y cos k y l Finding the eigenvalues: e λ det ( A λi) = = ( e + λ)( λ) ( cos k)( ) cos k λ ( e k ) λ ( e cosk) 0 ( ) ( = λ + = λ = e ± e + 4 e + cosk ( e ) ( e ) ± + + 4cosk λ = Now k > 0 0< e < 0< e < Recall that k = xi = ( 4i π ) δi, ( 4i π ) + εi ; ( i =,,, ) Examining the right-hand critical point in each pair, k = ( 4i π ) + ε, i ( 0< εi ) cos k = cos 4i π + ε = cos ε π = cos π ε = sin ε ε > 0 ( ) Therefore e 4( e cosk) ( e ) ) i i i i i + + > > 0 and the two eigenvalues are real and of opposite sign. These critical points are therefore all unstable saddle points. ) ()

17 ENGI Linear Approximation () Page 4.48 Example (continued) Examining the left-hand critical point in each pair, k = 4i π δ, 0< δ i ( i ) ( δi) ( δi ) ( δi ) cos k = cos 4i π = cos π = cos + π = sinδ k But k is the solution to e sin k = 0 ( ) e = + sin k = + sin 4i π δi = + sin δ i π = sin δ π i + = cosδ i δ i e + cos k ( cosδi) sinδi δi < 0 ( k e ) 4( e cosk) ( + + < e e ± e + 4 e + cosk) But λ = Therefore the eigenvalues are a real distinct positive pair and the singularity is an unstable node. The locations and nature of the first five critical points are listed here. k (δ i is small) x y λ λ Type unstable saddle point unstable node unstable saddle point unstable node unstable saddle point Here is a Maple session to create direction field plots for the first three critical points of the non-linear system. > with(detools): > DEplot([diff(x(t),t) = -y(t) - + exp(-x(t)), diff(y(t),t) = y(t) - sin(x(t))], [x(t),y(t)], t=-.., x= , y= , title=`example Exact Solution`); > DEplot([diff(x(t),t) = -y(t) - + exp(-x(t)), diff(y(t),t) = y(t) - sin(x(t))], [x(t),y(t)], t=-.., x=4.5..5, y= , title=`example Exact Solution`); i

18 ENGI Linear Approximation () Page 4.49 Example (continued) Saddle Point at (0, 0) Unstable Node near (4.57, 0.99), Saddle Point near (4.84, 0.99) Solution curves can be traced by following the arrows (which at every location point in the direction of dy dx ).

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.0 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential equations

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

Nonlinear Autonomous Systems of Differential

Nonlinear Autonomous Systems of Differential Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such

More information

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12 Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider

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

ENGI Duffing s Equation Page 4.65

ENGI Duffing s Equation Page 4.65 ENGI 940 4. - Duffing s Equation Page 4.65 4. Duffing s Equation Among the simplest models of damped non-linear forced oscillations of a mechanical or electrical system with a cubic stiffness term is Duffing

More information

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below. 54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and

More information

Math 3301 Homework Set Points ( ) ( ) I ll leave it to you to verify that the eigenvalues and eigenvectors for this matrix are, ( ) ( ) ( ) ( )

Math 3301 Homework Set Points ( ) ( ) I ll leave it to you to verify that the eigenvalues and eigenvectors for this matrix are, ( ) ( ) ( ) ( ) #7. ( pts) I ll leave it to you to verify that the eigenvalues and eigenvectors for this matrix are, λ 5 λ 7 t t ce The general solution is then : 5 7 c c c x( 0) c c 9 9 c+ c c t 5t 7 e + e A sketch of

More information

Math 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations

Math 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of

More information

Section 9.3 Phase Plane Portraits (for Planar Systems)

Section 9.3 Phase Plane Portraits (for Planar Systems) Section 9.3 Phase Plane Portraits (for Planar Systems) Key Terms: Equilibrium point of planer system yꞌ = Ay o Equilibrium solution Exponential solutions o Half-line solutions Unstable solution Stable

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. Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous.

Autonomous systems. Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous. Autonomous equations Autonomous systems Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous. i f i(x 1, x 2,..., x n ) for i 1,..., n As you

More information

Math 266: Phase Plane Portrait

Math 266: Phase Plane Portrait Math 266: Phase Plane Portrait Long Jin Purdue, Spring 2018 Review: Phase line for an autonomous equation For a single autonomous equation y = f (y) we used a phase line to illustrate the equilibrium solutions

More information

MATH 215/255 Solutions to Additional Practice Problems April dy dt

MATH 215/255 Solutions to Additional Practice Problems April dy dt . For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the

More information

Linear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form

Linear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form Linear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form d x x 1 = A, where A = dt y y a11 a 12 a 21 a 22 Here the entries of the coefficient matrix

More information

Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d

Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d July 6, 2009 Today s Session Today s Session A Summary of This Session: Today s Session A Summary of This

More information

Copyright (c) 2006 Warren Weckesser

Copyright (c) 2006 Warren Weckesser 2.2. PLANAR LINEAR SYSTEMS 3 2.2. Planar Linear Systems We consider the linear system of two first order differential equations or equivalently, = ax + by (2.7) dy = cx + dy [ d x x = A x, where x =, and

More information

Department of Mathematics IIT Guwahati

Department of Mathematics IIT Guwahati Stability of Linear Systems in R 2 Department of Mathematics IIT Guwahati A system of first order differential equations is called autonomous if the system can be written in the form dx 1 dt = g 1(x 1,

More information

2.10 Saddles, Nodes, Foci and Centers

2.10 Saddles, Nodes, Foci and Centers 2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one

More information

Math 273 (51) - Final

Math 273 (51) - Final Name: Id #: Math 273 (5) - Final Autumn Quarter 26 Thursday, December 8, 26-6: to 8: Instructions: Prob. Points Score possible 25 2 25 3 25 TOTAL 75 Read each problem carefully. Write legibly. Show all

More information

MATH 44041/64041 Applied Dynamical Systems, 2018

MATH 44041/64041 Applied Dynamical Systems, 2018 MATH 4441/6441 Applied Dynamical Systems, 218 Answers to Assessed Coursework 1 1 Method based on series expansion of the matrix exponential Let the coefficient matrix of α the linear system be A, then

More information

A plane autonomous system is a pair of simultaneous first-order differential equations,

A plane autonomous system is a pair of simultaneous first-order differential equations, Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium

More information

Homogeneous Constant Matrix Systems, Part II

Homogeneous Constant Matrix Systems, Part II 4 Homogeneous Constant Matrix Systems, Part II Let us now expand our discussions begun in the previous chapter, and consider homogeneous constant matrix systems whose matrices either have complex eigenvalues

More information

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point

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

AP Calculus (BC) Chapter 10 Test No Calculator Section. Name: Date: Period:

AP Calculus (BC) Chapter 10 Test No Calculator Section. Name: Date: Period: AP Calculus (BC) Chapter 10 Test No Calculator Section Name: Date: Period: Part I. Multiple-Choice Questions (5 points each; please circle the correct answer.) 1. The graph in the xy-plane represented

More information

Solutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.

Solutions to Dynamical Systems 2010 exam. Each question is worth 25 marks. Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the

More information

Phase portraits in two dimensions

Phase portraits in two dimensions Phase portraits in two dimensions 8.3, Spring, 999 It [ is convenient to represent the solutions to an autonomous system x = f( x) (where x x = ) by means of a phase portrait. The x, y plane is called

More information

B. Differential Equations A differential equation is an equation of the form

B. Differential Equations A differential equation is an equation of the form B Differential Equations A differential equation is an equation of the form ( n) F[ t; x(, xʹ (, x ʹ ʹ (, x ( ; α] = 0 dx d x ( n) d x where x ʹ ( =, x ʹ ʹ ( =,, x ( = n A differential equation describes

More information

4 Second-Order Systems

4 Second-Order Systems 4 Second-Order Systems Second-order autonomous systems occupy an important place in the study of nonlinear systems because solution trajectories can be represented in the plane. This allows for easy visualization

More information

Solutions Chapter 9. u. (c) u(t) = 1 e t + c 2 e 3 t! c 1 e t 3c 2 e 3 t. (v) (a) u(t) = c 1 e t cos 3t + c 2 e t sin 3t. (b) du

Solutions Chapter 9. u. (c) u(t) = 1 e t + c 2 e 3 t! c 1 e t 3c 2 e 3 t. (v) (a) u(t) = c 1 e t cos 3t + c 2 e t sin 3t. (b) du Solutions hapter 9 dode 9 asic Solution Techniques 9 hoose one or more of the following differential equations, and then: (a) Solve the equation directly (b) Write down its phase plane equivalent, and

More information

Nonlinear differential equations - phase plane analysis

Nonlinear differential equations - phase plane analysis Nonlinear differential equations - phase plane analysis We consider the general first order differential equation for y(x Revision Q(x, y f(x, y dx P (x, y. ( Curves in the (x, y-plane which satisfy this

More information

Solutions to Final Exam Sample Problems, Math 246, Spring 2011

Solutions to Final Exam Sample Problems, Math 246, Spring 2011 Solutions to Final Exam Sample Problems, Math 246, Spring 2 () Consider the differential equation dy dt = (9 y2 )y 2 (a) Identify its equilibrium (stationary) points and classify their stability (b) Sketch

More information

Created by T. Madas LINE INTEGRALS. Created by T. Madas

Created by T. Madas LINE INTEGRALS. Created by T. Madas LINE INTEGRALS LINE INTEGRALS IN 2 DIMENSIONAL CARTESIAN COORDINATES Question 1 Evaluate the integral ( x + 2y) dx, C where C is the path along the curve with equation y 2 = x + 1, from ( ) 0,1 to ( )

More information

First Midterm Exam Name: Practice Problems September 19, x = ax + sin x.

First Midterm Exam Name: Practice Problems September 19, x = ax + sin x. Math 54 Treibergs First Midterm Exam Name: Practice Problems September 9, 24 Consider the family of differential equations for the parameter a: (a Sketch the phase line when a x ax + sin x (b Use the graphs

More information

APPM 2360: Final Exam 10:30am 1:00pm, May 6, 2015.

APPM 2360: Final Exam 10:30am 1:00pm, May 6, 2015. APPM 23: Final Exam :3am :pm, May, 25. ON THE FRONT OF YOUR BLUEBOOK write: ) your name, 2) your student ID number, 3) lecture section, 4) your instructor s name, and 5) a grading table for eight questions.

More information

Calculus for the Life Sciences II Assignment 6 solutions. f(x, y) = 3π 3 cos 2x + 2 sin 3y

Calculus for the Life Sciences II Assignment 6 solutions. f(x, y) = 3π 3 cos 2x + 2 sin 3y Calculus for the Life Sciences II Assignment 6 solutions Find the tangent plane to the graph of the function at the point (0, π f(x, y = 3π 3 cos 2x + 2 sin 3y Solution: The tangent plane of f at a point

More information

Classification of Phase Portraits at Equilibria for u (t) = f( u(t))

Classification of Phase Portraits at Equilibria for u (t) = f( u(t)) Classification of Phase Portraits at Equilibria for u t = f ut Transfer of Local Linearized Phase Portrait Transfer of Local Linearized Stability How to Classify Linear Equilibria Justification of the

More information

Chapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems

Chapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.

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

Project 6.3 Predator-Prey and Your Own Game Preserve

Project 6.3 Predator-Prey and Your Own Game Preserve Project 6.3 Predator-Prey and Your Own Game Preserve The closed trajectories in the figure below represent periodic solutions of a typical predator-prey system, but provide no information as to the actual

More information

Def. (a, b) is a critical point of the autonomous system. 1 Proper node (stable or unstable) 2 Improper node (stable or unstable)

Def. (a, b) is a critical point of the autonomous system. 1 Proper node (stable or unstable) 2 Improper node (stable or unstable) Types of critical points Def. (a, b) is a critical point of the autonomous system Math 216 Differential Equations Kenneth Harris kaharri@umich.edu Department of Mathematics University of Michigan November

More information

ODE, part 2. Dynamical systems, differential equations

ODE, part 2. Dynamical systems, differential equations ODE, part 2 Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Dynamical systems, differential equations Consider a system of n first order equations du dt = f(u, t),

More information

Even-Numbered Homework Solutions

Even-Numbered Homework Solutions -6 Even-Numbered Homework Solutions Suppose that the matric B has λ = + 5i as an eigenvalue with eigenvector Y 0 = solution to dy = BY Using Euler s formula, we can write the complex-valued solution Y

More information

Dynamical Systems Solutions to Exercises

Dynamical Systems Solutions to Exercises Dynamical Systems Part 5-6 Dr G Bowtell Dynamical Systems Solutions to Exercises. Figure : Phase diagrams for i, ii and iii respectively. Only fixed point is at the origin since the equations are linear

More information

Half of Final Exam Name: Practice Problems October 28, 2014

Half of Final Exam Name: Practice Problems October 28, 2014 Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half

More information

C3 papers June 2007 to 2008

C3 papers June 2007 to 2008 physicsandmathstutor.com June 007 C3 papers June 007 to 008 1. Find the exact solutions to the equations (a) ln x + ln 3 = ln 6, (b) e x + 3e x = 4. *N6109A04* physicsandmathstutor.com June 007 x + 3 9+

More information

Math 331 Homework Assignment Chapter 7 Page 1 of 9

Math 331 Homework Assignment Chapter 7 Page 1 of 9 Math Homework Assignment Chapter 7 Page of 9 Instructions: Please make sure to demonstrate every step in your calculations. Return your answers including this homework sheet back to the instructor as a

More information

Lecture 13 - Wednesday April 29th

Lecture 13 - Wednesday April 29th Lecture 13 - Wednesday April 29th jacques@ucsdedu Key words: Systems of equations, Implicit differentiation Know how to do implicit differentiation, how to use implicit and inverse function theorems 131

More information

Exam 1 Review SOLUTIONS

Exam 1 Review SOLUTIONS 1. True or False (and give a short reason): Exam 1 Review SOLUTIONS (a) If the parametric curve x = f(t), y = g(t) satisfies g (1) = 0, then it has a horizontal tangent line when t = 1. FALSE: To make

More information

Section 5. Graphing Systems

Section 5. Graphing Systems Section 5. Graphing Systems 5A. The Phase Plane 5A-1. Find the critical points of each of the following non-linear autonomous systems. x = x 2 y 2 x = 1 x + y a) b) y = x xy y = y + 2x 2 5A-2. Write each

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

ENGI 9420 Lecture Notes 2 - Matrix Algebra Page Matrix operations can render the solution of a linear system much more efficient.

ENGI 9420 Lecture Notes 2 - Matrix Algebra Page Matrix operations can render the solution of a linear system much more efficient. ENGI 940 Lecture Notes - Matrix Algebra Page.0. Matrix Algebra A linear system of m equations in n unknowns, a x + a x + + a x b (where the a ij and i n n a x + a x + + a x b n n a x + a x + + a x b m

More information

Hello everyone, Best, Josh

Hello everyone, Best, Josh Hello everyone, As promised, the chart mentioned in class about what kind of critical points you get with different types of eigenvalues are included on the following pages (The pages are an ecerpt from

More information

Practice Problems for Final Exam

Practice Problems for Final Exam Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,

More information

Lecture 38. Almost Linear Systems

Lecture 38. Almost Linear Systems Math 245 - Mathematics of Physics and Engineering I Lecture 38. Almost Linear Systems April 20, 2012 Konstantin Zuev (USC) Math 245, Lecture 38 April 20, 2012 1 / 11 Agenda Stability Properties of Linear

More information

Nonlinear dynamics & chaos BECS

Nonlinear dynamics & chaos BECS Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes

More information

Volumes of Solids of Revolution Lecture #6 a

Volumes of Solids of Revolution Lecture #6 a Volumes of Solids of Revolution Lecture #6 a Sphereoid Parabaloid Hyperboloid Whateveroid Volumes Calculating 3-D Space an Object Occupies Take a cross-sectional slice. Compute the area of the slice. Multiply

More information

Math 4B Notes. Written by Victoria Kala SH 6432u Office Hours: T 12:45 1:45pm Last updated 7/24/2016

Math 4B Notes. Written by Victoria Kala SH 6432u Office Hours: T 12:45 1:45pm Last updated 7/24/2016 Math 4B Notes Written by Victoria Kala vtkala@math.ucsb.edu SH 6432u Office Hours: T 2:45 :45pm Last updated 7/24/206 Classification of Differential Equations The order of a differential equation is the

More information

Calculus and Differential Equations II

Calculus and Differential Equations II MATH 250 B Second order autonomous linear systems We are mostly interested with 2 2 first order autonomous systems of the form { x = a x + b y y = c x + d y where x and y are functions of t and a, b, c,

More information

FINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y.

FINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y. FINAL EXAM MATH0 Theory of Ordinary Differential Equations There are 5 problems on 2 pages. Spring 2009. 25 points Consider the linear plane autonomous system x + y x y. Find a fundamental matrix of the

More information

Math 1272 Solutions for Fall 2005 Final Exam

Math 1272 Solutions for Fall 2005 Final Exam Math 7 Solutions for Fall 5 Final Exam ) This fraction appears in Problem 5 of the undated-? exam; a solution can be found in that solution set. (E) ) This integral appears in Problem 6 of the Fall 4 exam;

More information

Lyapunov functions and stability problems

Lyapunov functions and stability problems Lyapunov functions and stability problems Gunnar Söderbacka, Workshop Ghana, 29.5-10.5, 2013 1 Introduction In these notes we explain the power of Lyapunov functions in determining stability of equilibria

More information

8.1 Bifurcations of Equilibria

8.1 Bifurcations of Equilibria 1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations

More information

ANSWERS Final Exam Math 250b, Section 2 (Professor J. M. Cushing), 15 May 2008 PART 1

ANSWERS Final Exam Math 250b, Section 2 (Professor J. M. Cushing), 15 May 2008 PART 1 ANSWERS Final Exam Math 50b, Section (Professor J. M. Cushing), 5 May 008 PART. (0 points) A bacterial population x grows exponentially according to the equation x 0 = rx, where r>0is the per unit rate

More information

1. Find and classify the extrema of h(x, y) = sin(x) sin(y) sin(x + y) on the square[0, π] [0, π]. (Keep in mind there is a boundary to check out).

1. Find and classify the extrema of h(x, y) = sin(x) sin(y) sin(x + y) on the square[0, π] [0, π]. (Keep in mind there is a boundary to check out). . Find and classify the extrema of hx, y sinx siny sinx + y on the square[, π] [, π]. Keep in mind there is a boundary to check out. Solution: h x cos x sin y sinx + y + sin x sin y cosx + y h y sin x

More information

MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. 3 2, 5 2 C) - 5 2

MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. 3 2, 5 2 C) - 5 2 Test Review (chap 0) Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Solve the problem. ) Find the point on the curve x = sin t, y = cos t, -

More information

C3 A Booster Course. Workbook. 1. a) Sketch, on the same set of axis the graphs of y = x and y = 2x 3. (3) b) Hence, or otherwise, solve the equation

C3 A Booster Course. Workbook. 1. a) Sketch, on the same set of axis the graphs of y = x and y = 2x 3. (3) b) Hence, or otherwise, solve the equation C3 A Booster Course Workbook 1. a) Sketch, on the same set of axis the graphs of y = x and y = 2x 3. b) Hence, or otherwise, solve the equation x = 2x 3 (3) (4) BlueStar Mathematics Workshops (2011) 1

More information

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland 4 May 2012 Because the presentation of this material

More information

ENGI Second Order Linear ODEs Page Second Order Linear Ordinary Differential Equations

ENGI Second Order Linear ODEs Page Second Order Linear Ordinary Differential Equations ENGI 344 - Second Order Linear ODEs age -01. Second Order Linear Ordinary Differential Equations The general second order linear ordinary differential equation is of the form d y dy x Q x y Rx dx dx Of

More information

Math 216 Final Exam 24 April, 2017

Math 216 Final Exam 24 April, 2017 Math 216 Final Exam 24 April, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that

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

Solved Examples. Given are two sets A {1, 2, -2, 3} and B = {1, 2, 3, 5}. Is the function f(x) = 2x - 1 defined from A to B?

Solved Examples. Given are two sets A {1, 2, -2, 3} and B = {1, 2, 3, 5}. Is the function f(x) = 2x - 1 defined from A to B? Solved Examples Example 1: Given are two sets A {1, 2, -2, 3} and B = {1, 2, 3, 5}. Is the function f(x) = 2x - 1 defined from A to B? Solution : Out of all the ordered pairs, the ordered pairs which are

More information

Find the general solution of the system y = Ay, where

Find the general solution of the system y = Ay, where Math Homework # March, 9..3. Find the general solution of the system y = Ay, where 5 Answer: The matrix A has characteristic polynomial p(λ = λ + 7λ + = λ + 3(λ +. Hence the eigenvalues are λ = 3and λ

More information

Linearization of Differential Equation Models

Linearization of Differential Equation Models Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking

More information

10.1 Review of Parametric Equations

10.1 Review of Parametric Equations 10.1 Review of Parametric Equations Recall that often, instead of representing a curve using just x and y (called a Cartesian equation), it is more convenient to define x and y using parametric equations

More information

Vectors, matrices, eigenvalues and eigenvectors

Vectors, matrices, eigenvalues and eigenvectors Vectors, matrices, eigenvalues and eigenvectors 1 ( ) ( ) ( ) Scaling a vector: 0.5V 2 0.5 2 1 = 0.5 = = 1 0.5 1 0.5 ( ) ( ) ( ) ( ) Adding two vectors: V + W 2 1 2 + 1 3 = + = = 1 3 1 + 3 4 ( ) ( ) a

More information

Non-Linear Dynamics Homework Solutions Week 6

Non-Linear Dynamics Homework Solutions Week 6 Non-Linear Dynamics Homework Solutions Week 6 Chris Small March 6, 2007 Please email me at smachr09@evergreen.edu with any questions or concerns reguarding these solutions. 6.8.3 Locate annd find the index

More information

Parametric Curves. Calculus 2 Lia Vas

Parametric Curves. Calculus 2 Lia Vas Calculus Lia Vas Parametric Curves In the past, we mostly worked with curves in the form y = f(x). However, this format does not encompass all the curves one encounters in applications. For example, consider

More information

36. Double Integration over Non-Rectangular Regions of Type II

36. Double Integration over Non-Rectangular Regions of Type II 36. Double Integration over Non-Rectangular Regions of Type II When establishing the bounds of a double integral, visualize an arrow initially in the positive x direction or the positive y direction. A

More information

ENGI Parametric Vector Functions Page 5-01

ENGI Parametric Vector Functions Page 5-01 ENGI 3425 5. Parametric Vector Functions Page 5-01 5. Parametric Vector Functions Contents: 5.1 Arc Length (Cartesian parametric and plane polar) 5.2 Surfaces of Revolution 5.3 Area under a Parametric

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

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

Phase Portraits of Nonlinear Differential Equations

Phase Portraits of Nonlinear Differential Equations ODE4-net.nb Phase Portraits of Nonlinear Differential Equations Nonlinear Differential Equations: x' = f (x, y) (1) Consider the system where f and g are functions of two y' = g(x, y) () variables x and

More information

Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories

Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 6, 2013 Outline

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

MathQuest: Differential Equations

MathQuest: Differential Equations MathQuest: Differential Equations Geometry of Systems 1. The differential equation d Y dt = A Y has two straight line solutions corresponding to [ ] [ ] 1 1 eigenvectors v 1 = and v 2 2 = that are shown

More information

You may use a calculator, but you must show all your work in order to receive credit.

You may use a calculator, but you must show all your work in order to receive credit. Math 2410-010/015 Exam II April 7 th, 2017 Name: Instructions: Key Answer each question to the best of your ability. All answers must be written clearly. Be sure to erase or cross out any work that you

More information

Week Quadratic forms. Principal axes theorem. Text reference: this material corresponds to parts of sections 5.5, 8.2,

Week Quadratic forms. Principal axes theorem. Text reference: this material corresponds to parts of sections 5.5, 8.2, Math 051 W008 Margo Kondratieva Week 10-11 Quadratic forms Principal axes theorem Text reference: this material corresponds to parts of sections 55, 8, 83 89 Section 41 Motivation and introduction Consider

More information

= e t sin 2t. s 2 2s + 5 (s 1) Solution: Using the derivative of LT formula we have

= e t sin 2t. s 2 2s + 5 (s 1) Solution: Using the derivative of LT formula we have Math 090 Midterm Exam Spring 07 S o l u t i o n s. Results of this problem will be used in other problems. Therefore do all calculations carefully and double check them. Find the inverse Laplace transform

More information

+ i. cos(t) + 2 sin(t) + c 2.

+ i. cos(t) + 2 sin(t) + c 2. MATH HOMEWORK #7 PART A SOLUTIONS Problem 7.6.. Consider the system x = 5 x. a Express the general solution of the given system of equations in terms of realvalued functions. b Draw a direction field,

More information

Homogeneous Constant Matrix Systems, Part II

Homogeneous Constant Matrix Systems, Part II 4 Homogeneous Constant Matri Systems, Part II Let us now epand our discussions begun in the previous chapter, and consider homogeneous constant matri systems whose matrices either have comple eigenvalues

More information

In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation.

In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation. 1 2 Linear Systems In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation 21 Matrix ODEs Let and is a scalar A linear function satisfies Linear superposition ) Linear

More information

Sample Solutions of Assignment 10 for MAT3270B

Sample Solutions of Assignment 10 for MAT3270B Sample Solutions of Assignment 1 for MAT327B 1. For the following ODEs, (a) determine all critical points; (b) find the corresponding linear system near each critical point; (c) find the eigenvalues of

More information

Vector Field Topology. Ronald Peikert SciVis Vector Field Topology 8-1

Vector Field Topology. Ronald Peikert SciVis Vector Field Topology 8-1 Vector Field Topology Ronald Peikert SciVis 2007 - Vector Field Topology 8-1 Vector fields as ODEs What are conditions for existence and uniqueness of streamlines? For the initial value problem i x ( t)

More information

Nonlinear Control Lecture 2:Phase Plane Analysis

Nonlinear Control Lecture 2:Phase Plane Analysis Nonlinear Control Lecture 2:Phase Plane Analysis Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 r. Farzaneh Abdollahi Nonlinear Control Lecture 2 1/53

More information

Math 308 Final Exam Practice Problems

Math 308 Final Exam Practice Problems Math 308 Final Exam Practice Problems This review should not be used as your sole source for preparation for the exam You should also re-work all examples given in lecture and all suggested homework problems

More information

Sample Solutions of Assignment 9 for MAT3270B

Sample Solutions of Assignment 9 for MAT3270B Sample Solutions of Assignment 9 for MAT370B. For the following ODEs, find the eigenvalues and eigenvectors, and classify the critical point 0,0 type and determine whether it is stable, asymptotically

More information

Solutions to Final Exam Sample Problems, Math 246, Spring 2018

Solutions to Final Exam Sample Problems, Math 246, Spring 2018 Solutions to Final Exam Sample Problems, Math 46, Spring 08 () Consider the differential equation dy dt = (9 y )y. (a) Find all of its stationary points and classify their stability. (b) Sketch its phase-line

More information

SECTION A. f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes.

SECTION A. f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes. SECTION A 1. State the maximal domain and range of the function f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes. 2. By evaluating f(0),

More information