Exam in TMA4195 Mathematical Modeling Solutions
|
|
- Maximilian McGee
- 5 years ago
- Views:
Transcription
1 Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Exam in TMA495 Mathematical Modeling Solutions Problem a Here x, y are two populations varying with time t. The equation for ẋ consists of: rx( x K logistic growth term of rate r and capacity K, axy c + x death term of rate a and dependent on both x and y, while the equation for ẏ consists of: my death term, of rate a, + bxy c + x growth term of rate b and dependent on both x and y. Possible models include: x prey, y predators, x humans, y bacteria, x lemmings, y foxes, x fish, y fishermen. b We have [ẋ ] ẏ [ x( x y = F (x, y := ] +x y( + x +x ( Equilibrium points (x e, y e are constant solutions of ( and therefore the solutions of ( x x y = 0, + x ( y + x = 0. + x
2 TMA495 Mathematical Modeling, 6..07, solutions Page of 9 The second equation is satisfied exactly when y = 0 or x =. For y = 0, the first equation is satisfied exactly when x = 0 or x =, while if x = the first equation is satisfied only for y =. Thus the equilibrium points are given by: (0, 0, (, 0, (,. OBS: By carelessly multiplying the equations with ( + x, one might conclude that (, 0 is also an equilibrium point, which is not the case (check!. To calculate the stability of the points we look at the Jacobian of F : DF (x, y := [ ( x y y (+x (+x x +x + x +x We calculate the eigenvalues λ, λ of DF in the different equilibrium points: [ ] 0 DF (0, 0 = where λ 0 =, λ =, [ ] /3 DF (, 0 = where λ 0 /3 =, λ = /3, [ ] / / DF (, = where λ 0 = i 4, λ = 7 4 i 4. The eigenvalues of the last matrix can be find by solving det(df (, λi = 0 ( / λ( λ ( / = 0 λ + λ + = 0. An equilibrium point is stable if max{reλ,reλ } < 0 and unstable if max{reλ,reλ } > 0. Consequently we see that (0, 0 and (, 0 are unstable while (, is stable. ] Problem Performing regular perturbation we write y as y(x = y 0 (x + εy (x + O(ε, where y 0, y are not dependent on ε. As f is smooth we can Taylor expand f about 0 to obtain f(εy = f(0 + f (0εy + O(ε, which in turn can be written as f(εy = f(0 + f (0ε ( y 0 (t + εy (t + O(ε + O(ε = f(0 + f (0εy 0 (t + O(ε.
3 TMA495 Mathematical Modeling, 6..07, solutions Page 3 of 9 Inserting these two representations in the differential equation we obtain while the initial conditions become y 0 + εy + f(0 + f (0εy 0 + O(ε = 0, y 0 (0 + εy (0 + O(ε = 0, y 0 ( + εy ( + O(ε = ε. Equating terms of different order in ε, we obtain the system of equations: O( : y 0 + f(0 = 0, y 0 (0 = 0, y 0 ( = 0, O(ε : y + f (0y 0 = 0, y 0 (0 = 0, y ( =. Solving first for y 0 we obtain y 0 (x = f(0x + Ax + B, where we must take A = f(0 and B = 0 to satisfy the boundary conditions for y 0. We can now insert for y 0 in the expression for y to obtain y + f (0f(0(x x = 0, = y (x = f (0f(0 ( x 4 x3 6 + Cx + D, and the boundary conditions gives C = + 4 f (0f(0 and D = 0. We conclude that y(x = y 0 (x + εy (x + O(ε = f(0(x x + ε ( x + 4 f (0f(0(x 4 x 3 + x + O(ε. Problem 3 We set z(t = ρ(x(t, t and find that the characteristic equations in this case are given by ẋ = j (z = e z, x(0 = x 0 and ż = 0, z(0 = ρ 0 (x 0. Solving we get z(t = ρ 0 (x 0 and x(t = e ρ 0(x 0 t + x 0. a In this first case, we obtain the two family of (x-characteristics x(t = x 0 + et, x 0 < 0, x(t = x 0 + t, x 0 > 0.
4 TMA495 Mathematical Modeling, 6..07, solutions Page 4 of 9 The left characteristics overtake the right characteristics (collisions in the region t < x < et, so the physical solution ρ is a shock solution, x < s(t, ρ(x, t = 0, x > s(t, where the shock curve s(t starts at s = 0 and satisfies the Rankine-Hugoniot condition: ṡ = j( j(0 0 = e, s(0 = 0 = s(t = (e t. b In the second case, the two families of characteristics are given by x(t = x 0 + t, x 0 < 0, x(t = x 0 + et, x 0 > 0. The left characteristics are slower than the right characteristics so there is a region t < x < et not reached by any characteristic (a dead sector. Hence the physical solution is given by a rarefaction wave ρ(x, t = ϕ(x/t. From the PDE for ρ we find that x t ϕ + e ϕ t ϕ = 0 = φ 0 We then get the solution: e ϕ = x t 0, x < t, ρ(x, t = ln(x/t, t < x < et,, et < x. = ϕ(x/t = ln(x/t. Problem 4 a For simplicity we drop the -notation in this problem and let x, y denote points in R 3. We fix a point x in the rock formation and consider the ball B r := B(x, r, centered at x with radius r > 0, whose volume will be denoted B r. By conservation of mass we have d φρ dy = (j ˆndσ. ( dt B r B r Since we have a bounded smooth domain and smooth integrands, we may change the order of differentiation and integration and use the divergence theorem to get d φρ dy = φρ t dy, dt B r B r (j ˆndσ = ( j dy. B r B r
5 TMA495 Mathematical Modeling, 6..07, solutions Page 5 of 9 Consequently, ( can be rewritten as B r (φρ t + j dy = 0. (3 Let f := φρ t + j, and divide by B r, and add and subtract f(x in (3 to get f(x + (f(y f(x dy = 0. (4 B r B r Since f is continuous f(y f(x for y B r and r, and hence for r, B (f(y f(x dy 0 r dy B r B r B r = φρ t (x, t + j(x, t = f(x 0. Sending r 0 we get equality and the first equation for ρ in the problem. [A rigorous argument which is not required for this exam, it given below: f(y f(x dy B r = ( f(x + s(y x (y x ds dy, B r B r B r 0 max f r ds dy = B r 0 R 3 B r B r max f r 0 as r 0. ] R 3 B r Combining Darcy s law and ideal gas law, we get j = ρ k µ p = ρk RT ρ, µ where we used that R, T are constants. Inserting this into the equation for ρ, we obtain ρ t = krt µφ b Let x = Xx, y = Y y, z = Zz, t = tt, ρ = ρρ, and (ρ ρ =: κ (ρ ρ. (5 x, y, z, t, ρ, ρ x, ρ xx, ρ y, ρ yy, ρ z, ρ zz (scaling assumption. By ideal gas law and max p = p, max ρ = max p RT = p RT. Thus ρ = p RT is a natural scale for ρ. Good scales for x, y, z are X = L, Y = Z = 500 m = εl.
6 TMA495 Mathematical Modeling, 6..07, solutions Page 6 of 9 The time scale is determined introducing the scaled quatities into (5 and balancing, ( ρ ρ t ρ t = κ = ρ t = κρt L L (ρρ x x + ρ ε L (ρρ y y + ρ ε L (ρρ z z ( (ρρ x x + ε [(ρρ y y + (ρρ z z ]. To arrive at the desired equation, we set t = L /κρ. c When x > s(t, ρ F = 0, and we immediately get (ρ F t = 0 = (ρ F (ρ F x x. When x < s(t, ρ F = ct 3 6 x t, and we calculate to conclude that (ρ F t = 3 ct x t, (ρ F x = 3 xt, (ρ F xx = 3 t, (ρ F (ρ F x x = ((ρ F x + ρ F (ρ F xx = 9 x t 3 ct x t = 3 ct x t = (ρ F t. Now we will show that ρ F solves the conservation law in integral form on any 0 < a < s(t < b: The time change of mass inside [a, b] is equal the influx at x = a (i.e. +j(a, t plus the influx at x = b (i.e. j(b, t, and the (scaled flux here is j = ρρ x (see also the hint. Hence we have d b ρ dx = j(b, t + j(a, t = (ρρ x (b, t (ρρ x (a, t. (6 dt a Since ρ F is a (classical solution for x (a, s(t and ρ = 0 for x [s(t, b, we get d ρ F dx = d ( s(t b ρ F dx + ρ F dx dt [a,b] dt a s(t }{{} = ρ F (s(t, ts (t + }{{} = =0 s(t a =0 s(t a (ρ F (ρ F x x dx = ρρ x s(t a = ρρ x b a, (ρ F t dx (Leibniz rule
7 TMA495 Mathematical Modeling, 6..07, solutions Page 7 of 9 and ρ F satisfies (6. Since ρ F = 0 for x > s(t, we first note that s(t s(t ρ F dx = ρ F dx =. Then since also ρ F 0 and lim t 0 s(t = 0, for any continuous function f, s(t ρ F (x, tf(x dx f(0 = ρ F (x, t ( f(x f(0 dx s(t s(t = max f(x f(0 ρ F (x, t dx x s(t s(t = max f(x f(0 0 as t 0. x s(t The function ρ F is a fundamental solution of equation (5 on the exam, if it is a solution with initial condition equal to the delta function. The meaning of the initial condition is that lim t 0 ρ F (x, tf(x = f(0 for all continuous functions f. Both requirements are satisfied by the previous parts of problem c, and hence ρ F is a fundamental solution. d The idea of the method of intermediate asymptotics is to rescale the problem for large (but finite! times in such a way that the rescaled intitial data approximates the δ- function. Then the solution of the rescaled problem will be close to the fundamental solution of this problem. Going back to original variables, we get an approximation valid for large times. Let t = 3 months and t = 5 days. Here we want to rescale the initial value problem such that ρ t + κ(ρ ρ x x = 0, ρ (x, t = ρ 0(x, (i the time scale is t = 3 months (large compared to t, (ii the scaled equation becomes equation (5 on the exam (to be able to use part c, (iii the scaled initial data ρ(x, 0 δ(x, the δ-function. That is, we seek a scaled problem of the form ρ t + (ρρ x x = 0, ρ(x, 0 δ(x, (7 If we can do that, we can conclude (by stability of the PDE that the scaled solution ρ ρ F, where ρ F is the fundamental solution given in part c. Going back to original variables, we find an approximation of the solution ρ at t = t = 3 months.
8 TMA495 Mathematical Modeling, 6..07, solutions Page 8 of 9 Let us now show that (7 holds after scaling. Take ρ = ρρ, t = T t, and x = Xx. The time scale T = t, and we choose X and ρ such that (7 holds. From the equation, we find that ρ t + κ(ρ ρ x x = 0 = ρ t ρ t + κρ X (ρρ x x = 0 = ρ t + κρ t X (ρρ x x = 0, and then we get (7 if κρ t X =. Next we look at the initial condition. Recall that ρ 0 = 0 for x > l := 5 m and of total mass M := ρ dx. We want to choose a scaling such that ρ dx = and ρ(x, 0 = 0 for x > ε for some 0 < ε, or in other words, ρ(x, 0 δ(x. The integral conditions is satisfied if M = ρ dx = ρx ρ dx = ρx. We now have two equations for two unknowns X and ρ, and the solution is ( M 3 ρ =, X = ( κm t 3. κ t Note that ρ(x, t t = ρ ρ (Xx, t = 0 for Xx l or x l X. Since t = 3months s, a quick calculation shows that X = ( 3 κmt ( m = 000 m 5 m = l. Hence l X, and since also t t, it follows that ρ(x, 0 ρ(x, t t δ(x. We conclude that an approximation of ρ at t = t = 3 months, is given by ( x ρ (x, t ρρ F X, t t ( M ( 3 x = ρ F κ t (κm t 3,. As the rock formation is surrounded by impermeable rock, there can be no flux at the boundaries, x = ± L. Thus the boundary conditions are 0 = j (± L, T = ρ (± L, T ρ x (± L, T.
9 TMA495 Mathematical Modeling, 6..07, solutions Page 9 of 9 This is indeed satisfied by our approximation since L X = L (κm t 3 =.5 >.65 (6c = s(, and then since ρ F (x, t = 0 for x > s(t, ρ (± L, t = ρρ F ( L X, = 0.
Examination paper for TMA4195 Mathematical Modeling
Department of Mathematical Sciences Examination paper for TMA4195 Mathematical Modeling Academic contact during examination: Espen R. Jakobsen Phone: 73 59 35 12 Examination date: December 16, 2017 Examination
More informationTMA4195 Mathematical modelling Autumn 2012
Norwegian University of Science an Technology Department of Mathematical Sciences TMA495 Mathematical moelling Autumn 202 Solutions to exam December, 202 Dimensional matrix A: τ µ u r m - - s 0-2 - - 0
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationSTABILITY. 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 informationDimensional Analysis - Concepts
Dimensional Analysis - Concepts Physical quantities: R j = v(r j )[R j ] = value unit, j = 1,..., m. Units: Dimension matrix of R 1,, R m : A = Change of units change of values: [R j ] = F a 1j 1 F a nj
More informationEG4321/EG7040. Nonlinear Control. Dr. Matt Turner
EG4321/EG7040 Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear [System Analysis] and Control Dr. Matt
More informationEECS 700: Exam # 1 Tuesday, October 21, 2014
EECS 700: Exam # 1 Tuesday, October 21, 2014 Print Name and Signature The rules for this exam are as follows: Write your name on the front page of the exam booklet. Initial each of the remaining pages
More informationCoupled 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 informationFinal Exam. Monday March 19, 3:30-5:30pm MAT 21D, Temple, Winter 2018
Name: Student ID#: Section: Final Exam Monday March 19, 3:30-5:30pm MAT 21D, Temple, Winter 2018 Show your work on every problem. orrect answers with no supporting work will not receive full credit. Be
More informationInteractions. 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 informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More informationMATH 220: MIDTERM OCTOBER 29, 2015
MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
More information154 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 informationExamination paper for TMA4195 Mathematical Modeling
Department of Mathematical Sciences Examination paper for TMA4195 Mathematical Modeling Academic contact during examination: Elena Celledoni Phone: 48238584, 73593541 Examination date: 11th of December
More informationDynamical 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 informationMath 126 Final Exam Solutions
Math 126 Final Exam Solutions 1. (a) Give an example of a linear homogeneous PE, a linear inhomogeneous PE, and a nonlinear PE. [3 points] Solution. Poisson s equation u = f is linear homogeneous when
More informationSOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More informationMath 46, Applied Math (Spring 2008): Final
Math 46, Applied Math (Spring 2008): Final 3 hours, 80 points total, 9 questions, roughly in syllabus order (apart from short answers) 1. [16 points. Note part c, worth 7 points, is independent of the
More informationMATH 173: PRACTICE MIDTERM SOLUTIONS
MATH 73: PACTICE MIDTEM SOLUTIONS This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve all of them. Write your solutions to problems and in blue book #, and your
More informationEntrance 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 informationMath Exam III - Spring
Math 3 - Exam III - Spring 8 This exam contains 5 multiple choice questions and hand graded questions. The multiple choice questions are worth 5 points each and the hand graded questions are worth a total
More informationPreliminary Exam 2018 Solutions to Morning Exam
Preliminary Exam 28 Solutions to Morning Exam Part I. Solve four of the following five problems. Problem. Consider the series n 2 (n log n) and n 2 (n(log n)2 ). Show that one converges and one diverges
More informationx + ye z2 + ze y2, y + xe z2 + ze x2, z and where T is the
1.(8pts) Find F ds where F = x + ye z + ze y, y + xe z + ze x, z and where T is the T surface in the pictures. (The two pictures are two views of the same surface.) The boundary of T is the unit circle
More informationPractice 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 informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationM340 HW 2 SOLUTIONS. 1. For the equation y = f(y), where f(y) is given in the following plot:
M340 HW SOLUTIONS 1. For the equation y = f(y), where f(y) is given in the following plot: (a) What are the critical points? (b) Are they stable or unstable? (c) Sketch the solutions in the ty plane. (d)
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
More informationSolutions 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 informationExam 4 SCORE. MA 114 Exam 4 Spring Section and/or TA:
Exam 4 Name: Section and/or TA: Last Four Digits of Student ID: Do not remove this answer page you will return the whole exam. You will be allowed two hours to complete this test. No books or notes may
More informationMATH H53 : Final exam
MATH H53 : Final exam 11 May, 18 Name: You have 18 minutes to answer the questions. Use of calculators or any electronic items is not permitted. Answer the questions in the space provided. If you run out
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More informationGeometry and Motion Selected answers to Sections A and C Dwight Barkley 2016
MA34 Geometry and Motion Selected answers to Sections A and C Dwight Barkley 26 Example Sheet d n+ = d n cot θ n r θ n r = Θθ n i. 2. 3. 4. Possible answers include: and with opposite orientation: 5..
More informationChapter #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 informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationDifferential equations
Differential equations Math 27 Spring 2008 In-term exam February 5th. Solutions This exam contains fourteen problems numbered through 4. Problems 3 are multiple choice problems, which each count 6% of
More informationMath Partial Differential Equations 1
Math 9 - Partial Differential Equations Homework 5 and Answers. The one-dimensional shallow water equations are h t + (hv) x, v t + ( v + h) x, or equivalently for classical solutions, h t + (hv) x, (hv)
More informationProblem 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 informationENGI Partial Differentiation Page y f x
ENGI 344 4 Partial Differentiation Page 4-0 4. Partial Differentiation For functions of one variable, be found unambiguously by differentiation: y f x, the rate of change of the dependent variable can
More informationDirection of maximum decrease = P
APPM 35 FINAL EXAM PING 15 INTUTION: Electronic devices, books, and crib sheets are not permitted. Write your name and your instructor s name on the front of your bluebook. Work all problems. how your
More informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
More informationModeling with differential equations
Mathematical Modeling Lia Vas Modeling with differential equations When trying to predict the future value, one follows the following basic idea. Future value = present value + change. From this idea,
More informationIntegration - Past Edexcel Exam Questions
Integration - Past Edexcel Exam Questions 1. (a) Given that y = 5x 2 + 7x + 3, find i. - ii. - (b) ( 1 + 3 ) x 1 x dx. [4] 2. Question 2b - January 2005 2. The gradient of the curve C is given by The point
More informationMATH 32A: MIDTERM 2 REVIEW. sin 2 u du z(t) = sin 2 t + cos 2 2
MATH 3A: MIDTERM REVIEW JOE HUGHES 1. Curvature 1. Consider the curve r(t) = x(t), y(t), z(t), where x(t) = t Find the curvature κ(t). 0 cos(u) sin(u) du y(t) = Solution: The formula for curvature is t
More information16.2 Line Integrals. Lukas Geyer. M273, Fall Montana State University. Lukas Geyer (MSU) 16.2 Line Integrals M273, Fall / 21
16.2 Line Integrals Lukas Geyer Montana State University M273, Fall 211 Lukas Geyer (MSU) 16.2 Line Integrals M273, Fall 211 1 / 21 Scalar Line Integrals Definition f (x) ds = lim { s i } N f (P i ) s
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationSalmon: Lectures on partial differential equations
6. The wave equation Of the 3 basic equations derived in the previous section, we have already discussed the heat equation, (1) θ t = κθ xx + Q( x,t). In this section we discuss the wave equation, () θ
More informationPH.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 informationMultivariable Calculus Notes. Faraad Armwood. Fall: Chapter 1: Vectors, Dot Product, Cross Product, Planes, Cylindrical & Spherical Coordinates
Multivariable Calculus Notes Faraad Armwood Fall: 2017 Chapter 1: Vectors, Dot Product, Cross Product, Planes, Cylindrical & Spherical Coordinates Chapter 2: Vector-Valued Functions, Tangent Vectors, Arc
More informationMS 2001: Test 1 B Solutions
MS 2001: Test 1 B Solutions Name: Student Number: Answer all questions. Marks may be lost if necessary work is not clearly shown. Remarks by me in italics and would not be required in a test - J.P. Question
More informationOrdinary 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 informationHomework 2. Due Friday, July We studied the logistic equation in class as a model of population growth. It is given by dn dt = rn 1 N
Problem 1 (10 points) Homework Due Friday, July 7 017 We studied the logistic equation in class as a model of population growth. It is given by dn dt = rn 1 N, (1) K with N(0) = N 0. (a) Make the substitutions
More informationWorkshop on Theoretical Ecology and Global Change March 2009
2022-3 Workshop on Theoretical Ecology and Global Change 2-18 March 2009 Stability Analysis of Food Webs: An Introduction to Local Stability of Dynamical Systems S. Allesina National Center for Ecological
More informationPractice Problems for Exam 3 (Solutions) 1. Let F(x, y) = xyi+(y 3x)j, and let C be the curve r(t) = ti+(3t t 2 )j for 0 t 2. Compute F dr.
1. Let F(x, y) xyi+(y 3x)j, and let be the curve r(t) ti+(3t t 2 )j for t 2. ompute F dr. Solution. F dr b a 2 2 F(r(t)) r (t) dt t(3t t 2 ), 3t t 2 3t 1, 3 2t dt t 3 dt 1 2 4 t4 4. 2. Evaluate the line
More informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
More informationVariational Methods & Optimal Control
Variational Methods & Optimal Control lecture 02 Matthew Roughan Discipline of Applied Mathematics School of Mathematical Sciences University of Adelaide April 14, 2016
More informationFIRST-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 informationHomework #3 Solutions
Homework #3 Solutions Math 82, Spring 204 Instructor: Dr. Doreen De eon Exercises 2.2: 2, 3 2. Write down the heat equation (homogeneous) which corresponds to the given data. (Throughout, heat is measured
More informationMTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.
MTH4101 CALCULUS II REVISION NOTES 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) 1.1 Introduction Types of numbers (natural, integers, rationals, reals) The need to solve quadratic equations:
More informationThe functions in all models depend on two variables: time t and spatial variable x, (x, y) or (x, y, z).
Review of Multi-variable calculus: The functions in all models depend on two variables: time t and spatial variable x, (x, y) or (x, y, z). The spatial variable represents the environment where the species
More informationTMA 4195 Mathematical Modeling, 30 November 2006 Solution with additional comments
TMA 495 Mathematical Modeling, 30 November 006 Solution with additional comments Problem A cookbook states a cooking time of 3.5 hours at the oven temperature T oven For a 0 kg (American) turkey, 7 hours
More informationSolution of Additional Exercises for Chapter 4
1 1. (1) Try V (x) = 1 (x 1 + x ). Solution of Additional Exercises for Chapter 4 V (x) = x 1 ( x 1 + x ) x = x 1 x + x 1 x In the neighborhood of the origin, the term (x 1 + x ) dominates. Hence, the
More information1 (t + 4)(t 1) dt. Solution: The denominator of the integrand is already factored with the factors being distinct, so 1 (t + 4)(t 1) = A
Calculus Topic: Integration of Rational Functions Section 8. # 0: Evaluate the integral (t + )(t ) Solution: The denominator of the integrand is already factored with the factors being distinct, so (t
More informationMath 142, Final Exam. 12/7/10.
Math 4, Final Exam. /7/0. No notes, calculator, or text. There are 00 points total. Partial credit may be given. Write your full name in the upper right corner of page. Number the pages in the upper right
More information1.2. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy
.. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy = f(x, y). In this section we aim to understand the solution
More informationHalf 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 informationStreamline calculations. Lecture note 2
Streamline calculations. Lecture note 2 February 26, 2007 1 Recapitulation from previous lecture Definition of a streamline x(τ) = s(τ), dx(τ) dτ = v(x,t), x(0) = x 0 (1) Divergence free, irrotational
More information1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system:
1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system: ẋ = y x 2, ẏ = z + xy, ż = y z + x 2 xy + y 2 + z 2 x 4. (ii) Determine
More information2D-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 informationTHE WAVE EQUATION. d = 1: D Alembert s formula We begin with the initial value problem in 1 space dimension { u = utt u xx = 0, in (0, ) R, (2)
THE WAVE EQUATION () The free wave equation takes the form u := ( t x )u = 0, u : R t R d x R In the literature, the operator := t x is called the D Alembertian on R +d. Later we shall also consider the
More informationAPPM 2350 FINAL EXAM FALL 2017
APPM 25 FINAL EXAM FALL 27. ( points) Determine the absolute maximum and minimum values of the function f(x, y) = 2 6x 4y + 4x 2 + y. Be sure to clearly give both the locations and values of the absolute
More informationFinal exam (practice) UCLA: Math 31B, Spring 2017
Instructor: Noah White Date: Final exam (practice) UCLA: Math 3B, Spring 207 This exam has 8 questions, for a total of 80 points. Please print your working and answers neatly. Write your solutions in the
More informationDO NOT BEGIN THIS TEST UNTIL INSTRUCTED TO START
Math 265 Student name: KEY Final Exam Fall 23 Instructor & Section: This test is closed book and closed notes. A (graphing) calculator is allowed for this test but cannot also be a communication device
More informationFinal Exam Due on Sunday 05/06
Final Exam Due on Sunday 05/06 The exam should be completed individually without collaboration. However, you are permitted to consult with the textbooks, notes, slides and even internet resources. If you
More information11 Chaos in Continuous Dynamical Systems.
11 CHAOS IN CONTINUOUS DYNAMICAL SYSTEMS. 47 11 Chaos in Continuous Dynamical Systems. Let s consider a system of differential equations given by where x(t) : R R and f : R R. ẋ = f(x), The linearization
More informationDynamical Systems and Chaos Part II: Biology Applications. Lecture 6: Population dynamics. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part II: Biology Applications Lecture 6: Population dynamics Ilya Potapov Mathematics Department, TUT Room TD325 Living things are dynamical systems Dynamical systems theory
More informationExercises for Multivariable Differential Calculus XM521
This document lists all the exercises for XM521. The Type I (True/False) exercises will be given, and should be answered, online immediately following each lecture. The Type III exercises are to be done
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Spring 208 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationExercise Set 4. D s n ds + + V. s dv = V. After using Stokes theorem, the surface integral becomes
Exercise Set Exercise - (a) Let s consider a test volum in the pellet. The substract enters the pellet by diffusion and some is created and disappears due to the chemical reaction. The two contribute to
More informationDRAFT - Math 101 Lecture Note - Dr. Said Algarni
2 Limits 2.1 The Tangent Problems The word tangent is derived from the Latin word tangens, which means touching. A tangent line to a curve is a line that touches the curve and a secant line is a line that
More informationAn Application of Perturbation Methods in Evolutionary Ecology
Dynamics at the Horsetooth Volume 2A, 2010. Focused Issue: Asymptotics and Perturbations An Application of Perturbation Methods in Evolutionary Ecology Department of Mathematics Colorado State University
More informationSolutions to Section 1.1
Solutions to Section True-False Review: FALSE A derivative must involve some derivative of the function y f(x), not necessarily the first derivative TRUE The initial conditions accompanying a differential
More informationSystems of Ordinary Differential Equations
Systems of Ordinary Differential Equations Systems of ordinary differential equations Last two lectures we have studied models of the form y (t) F (y), y(0) y0 (1) this is an scalar ordinary differential
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationMath 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 information1. (30 points) In the x-y plane, find and classify all local maxima, local minima, and saddle points of the function. f(x, y) = 3y 2 2y 3 3x 2 + 6xy.
APPM 35 FINAL EXAM FALL 13 INSTUTIONS: Electronic devices, books, and crib sheets are not permitted. Write your name and your instructor s name on the front of your bluebook. Work all problems. Show your
More informationPH.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 informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More information2.152 Course Notes Contraction Analysis MIT, 2005
2.152 Course Notes Contraction Analysis MIT, 2005 Jean-Jacques Slotine Contraction Theory ẋ = f(x, t) If Θ(x, t) such that, uniformly x, t 0, F = ( Θ + Θ f x )Θ 1 < 0 Θ(x, t) T Θ(x, t) > 0 then all solutions
More informationGlobal Stability Analysis on a Predator-Prey Model with Omnivores
Applied Mathematical Sciences, Vol. 9, 215, no. 36, 1771-1782 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.512 Global Stability Analysis on a Predator-Prey Model with Omnivores Puji Andayani
More informationNonlinear Control. Nonlinear Control Lecture # 2 Stability of Equilibrium Points
Nonlinear Control Lecture # 2 Stability of Equilibrium Points Basic Concepts ẋ = f(x) f is locally Lipschitz over a domain D R n Suppose x D is an equilibrium point; that is, f( x) = 0 Characterize and
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationLAURENTIAN 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 informationMath 5588 Final Exam Solutions
Math 5588 Final Exam Solutions Prof. Jeff Calder May 9, 2017 1. Find the function u : [0, 1] R that minimizes I(u) = subject to u(0) = 0 and u(1) = 1. 1 0 e u(x) u (x) + u (x) 2 dx, Solution. Since the
More informationSOLUTIONS, PROBLEM SET 11
SOLUTIONS, PROBLEM SET 11 1 In this problem we investigate the Lagrangian formulation of dynamics in a rotating frame. Consider a frame of reference which we will consider to be inertial. Suppose that
More informationThe death of an epidemic
LECTURE 2 Equilibrium Stability Analysis & Next Generation Method The death of an epidemic In SIR equations, let s divide equation for dx/dt by dz/ dt:!! dx/dz = - (β X Y/N)/(γY)!!! = - R 0 X/N Integrate
More informationMath 220A - Fall 2002 Homework 5 Solutions
Math 0A - Fall 00 Homework 5 Solutions. Consider the initial-value problem for the hyperbolic equation u tt + u xt 0u xx 0 < x 0 u t (x, 0) ψ(x). Use energy methods to show that the domain of dependence
More information