Solution Manual for: Numerical Computing with MATLAB by Cleve B. Moler
|
|
- Angelica Hancock
- 5 years ago
- Views:
Transcription
1 Solution Manual for: Numerical Computing with MATLAB by Cleve B. Moler John L. Weatherwax July 5, 007 Chapter 7 (Ordinary Differential Equations) Problem 7.1 Defining the vector y as y Then we have for its time derivative the following v dt u v ü v u v u y y 1 (1) y 3 y 4 + cos(r) 1+t sin(r) 1+t () Where r given in terms of components of y is r y4 + y 3 (3) Given the initial conditions in terms of u, v, u, and v we have 0 y (4) 0 wax@alum.mit.edu 1
2 Problem 7. See the matlab code in the file prob 7.m. Problem 7.3 The algorithm BS3 can be represented by the following sequence of steps s 1 f(t n, y n ) (5) s f(t n + h, y n + h s 1) (6) s 3 f(t n h, y n s ) (7) y n+1 y n + h 9 (s 1 + 3s + 4s 3 ) (8) Part (a): The specific ODE s (and their exact solutions) to consider are given by dt dt dt dt 1 y t + C (9) t y t + C (10) t y t3 3 + C (11) t 3 y t4 4 + C (1) For the first system, /dt 1, so f(t, y) 1 then the BS3 procedure gives the following s 1 1 (13) s 1 (14) s 3 1 (15) y n+1 y n + h 9 ( ) y n + h (16) Does this equal the true solution at t n + h? We evaluate the above to obtain y(t n + h) (t n + h) + C t n + C + h y n + h (17) and the answer is yes. For the second function f(t, y) t we see that the BS3 procedure is s 1 t n (18) s t n + h (19) s 3 t n h (0) y n+1 y n + h (t n + 3(t n + h ) + 4(t n + 3 h)) (1) 4
3 Simplifying this last equation we have y n + h 9 (9t n + ( 3 + 3)h) () y n + ht n + h 9 (9 )h (3) y n + ht n + h Again we ask, does this equal the true solution at t n+1? Lets check by evaluating (4) y(t n+1 ) y(t n + h) 1 (t n + h) + C t n + C + ht n + h y(t n) + ht n + h and we see that the answer is yes. For the second function f(t, y) t, then the BS3 gives (5) s 1 t n (6) s (t n + h ) (7) s 3 (t n h) (8) y n+1 y n + h 9 (t n + 3(t n + h ) + 4(t n h) ) (9) Now simplifying the last equation we have y n+1 y n + h 9 (t n + (3 + 6)ht n + ( )h ) (30) y n + h 9 (9t n + 9 ht n h ) (31) y n + ht n + h t n + h Does this equal the exact solution at t n + h? To decide we compute (3) y(t n + h) (t n + h) 3 + C (t3 n + 3t n h + 3t nh + h 3 ) + C (33) t3 n 3 + C + (t n h + t nh + h3 3 ) (34) and we see that the answer is yes. Finally, for the function f(t, y) t 3, the BS3 algorithm gives s 1 t 3 n (35) s (t n + h )3 t 3 n + 3 t n h t nh h3 (36) s 3 (t n h)3 t 3 n t n h t nh h3 (37) y n+1 y n + h 9 (t3 n + 3(t3 n + 3 t n h t nh h3 ) (38) +4(t 3 n t n h t nh h3 )) (39)
4 Simplifying this last term gives y n+1 y n + ht 3 n + 3 t n h + t n h 3 + h4 6 The true solution at t n + h however is equal to (40) y(t n + h) (t n + h) 4 + C (t4 n + 4t3 n h + 6t n h + 4t n h 3 + h 4 ) + C (41) 1 4 t4 n + C + t 3 nh + 3 t nh + t n h 3 + h4 4 (4) y n + ht 3 n + 3 t nh + t n h 3 + h4 4. (43) Since this differs from the the approximate solution by ( ) h 4 we see that this method is not exact for f(t) t 3. Part (b): The BS3 error estimate is given by e n+1 h 7 ( 5s 1 + 6s + 8s 3 9s 4 ) (44) where s i for i 1,, 3, 4 is computed as above. When one computes s 1, s, s 3, and s 4 for f(t) 1 and f(t) t one obtains a zero estimation error. Considering the function f(t) t we have an estimated error given by e n+1 h 7 ( 5t n + 6(t n + h ) + 8(t n h) 9(t n + h) ) (45) when one expands the above expression one obtains e n+1 h3 4 which is not zero. Thus the estimated error is exact for f(t) 1 and f(t) t only. (46) Problem 7.4 See the matlab code in the file prob 7 4.m. Problem 7.5 Part (a): See the matlab code in the file myrk4.m. Part (b): The local truncation error of RK4 is given by O(h 4+1 ), so the global truncation error should be O(h 4 ). When you decrease the step size h by this should reduce the global error by 1/ 4 1/16. An experiment illustrating this behaviour can be found in
5 the routine/code prob 7 5.m. There we sequentially decrease the stepsize parameter h and compute the global error between the final timestep and the known true solution value. We then plot the log of stepsize parameter v.s. the log of the global error. On this plot one can see that in this space we have a linear decrease in error representing a powerlaw. A linear function is fitted to to this data with the Matlab function polyfit. The slope of this line is given by 4 empirically displaying the known relation between the error and the stepsize of error h 4 Part (c): For this problem formulation our differential equation is given by ÿ y. The system formulation for this equation is given by defining [ ] [ ] y1 (t) y(t) Y (47) y (t) ẏ(t) then d dt Y (t) d [ y(t) dt ẏ(t) ] [ ẏ(t) ÿ(t) ] so in total in terms of y 1 and y the differential equation is [ ] [ ] d y1 y dt y y 1 [ y (t) y 1 (t) ] (48) (49) This is implemented in the Matlab code prob 7 5 fn.m. In addition, in the driver routine prob 7 5.m we have implemented the numerical experiments mentioned in the book and have verified the statements there. Problem 7.6 Part (a): The differential equation ẏ 1000y sin(t) + cos(t) with y(0) 1 (50) has a homogeneous solution given by y(t) Ce 1000t and a particular solution given by a function of the form A sin(t)+b cos(t). By putting this particular solution into the equation above we can solve for the coefficients A and B. When we do this we get the following A cos(t) B sin(t) 1000A sin(t) 1000B cos(t) sin(t) + cos(t) (51) Grouping the coefficients of cos(t) and sin(t) we obtain (A B 1) cos(t) + ( B A 1000) sin(t) 0. (5) When we require that the coefficients of cos(t) and sin(t) vanish we obtain the following system of equations A B 1 0 (53) B A (54)
6 The solution to the first equation in terms of B is given by A B. When substituted into the second equation we obtain 1000(1 1000B) B 1000 (55) Which has as its solution B 0. Thus A 1 giving in total the solution y of y(t) sin(t) + Ce 1000t (56) using the initial condition y(0) 1 we obtain C 1 and the total analytical solution is given by y(t) e 1000t + sin(t) (57) I should mention that heuristically this problem is midly stiff because of the various time scales present in the solution. The exponential term has a timescale O(1/1000) while the sin term has a timescale like O(1) or effectivly three orders of magnitude different. Part (b)-(f): See the routines prob 7 6.m and prob 7 6 fn.m for an implementation of this problem. Problem 7.7 Part (a): Each equation has as its solution y(t) sin(t). Part (b): The corresponding first order system for the third given differential equations is (defining y 1 (t) y(t)) [ ] [ ] [ ] [ ] d y1 y 0 1 y1 (58) dt y y y The corresponding system for the fourth differential equation is given by [ ] [ ] [ ] [ ] [ ] d y1 y 0 1 y1 0 + dt y sin(t) 0 0 y sin(t) Note that the first representation above is autonomous while the second is not. Part (c): The Jacobian for each of the functions f involved in the definitions of each differential equation are given by f 1 y f y f 3 y f 4 y 0 y 1 y [ ] [ ] From the above we can see that as t π that WWX finish (59)
7 Problem 7.9 Consider the matrix A given by A 0 η 0 σ σ η ρ 1 (60) Expanding the determinant of the above matrix we have or or A σ σ ρ 1 η 0 η σ σ (61) A (σ σρ) η(ση) (6) A η ρ + (63) Setting this equal to zero to find the values where the matrix is singular we see that η ± (ρ 1). (64) Which was the expression we were looking to find. To find the vectors in the nullity of A when η is one of these special values, insert this into the matrix above and reduce it to row-echelon form. Inserting this expression into A means we are looking for a vector v such that Av 0 ± (ρ 1) 0 σ σ (ρ 1) ρ 1 v 1 v v 3 0 (65) dividing the first row in the above expression by and the second row by σ we obtain (ρ 1) 1 0 v v 0 (66) (ρ 1) ρ 1 v 3 Now multiplying the first row by ± (ρ 1) and adding it to the third row we obtain (ρ 1) ρ 1 (ρ 1) (ρ 1) v 1 v v 3 0 (67) which simplifies to (ρ 1) ρ ρ v 1 v v 3 0. (68)
8 which demonstrates that the third row is indeed a multiple of the second row. Eliminating this row we obtain the following requirements among the coefficients of v (ρ 1) v 1 ± v 3 (69) v v 3 (70) Letting v 3 η (ρ 1) and assuming that ρ 1 > 0 we see that the above null vector has components given by as desired. v 1 ρ 1 (71) v η (7) v 3 η (73) Problem 7.10 (the Jacobian of the Lorenz equations) The Lorenz equations are given by ẏ Ay, with y 1 (t) y(t) y (t) and A y 3 (t) 0 y (t) 0 σ σ y (t) ρ 1 Performing the matrix multiplication we have for the full non-linear system representing the Lorenz equations, the following y 1 y 1 (t) + y (t)y 3 (t) y σy (t) + σy 3 (t) y 3 y 1 (t)y (t) + ρy (t) y 3 (t). This expression makes it easier to compute the required Jacobian of our mapping f( ). We find that y 3 (t) y (t) J 0 σ σ. y (t) y 1 (t) + ρ 1 The fixed points of the Lorentz equations are given by solving for y in f(y) 0, which gives the following system of equations y 1 + y y 3 0 σy + σy 3 0 y 1 y + ρy y 3 0. The second equation above has the simple solution of y y 3, which when used in the first and third equaions above gives y 1 + y 0 y 1 y + ρy y 1 0.,
9 The last equation above can be factored as y ( y 1 + ρ 1) 0, which has as solutions y 0, and y 1 ρ 1. Following the consequences of each of these choices, if y 0, then y 1 0, while if y 1 ρ 1 then the first equation above gives (ρ 1) + y 0, which has a solution given by y ± (ρ 1). Thus in summary then the three fixed points of the Jacobian of the Lorentz equations are (y 1, y, y 3 ) (0, 0, 0) (y 1, y, y 3 ) (ρ 1, (ρ 1), (ρ 1)) (y 1, y, y 3 ) (0, 0, 0).
Ordinary Differential Equations I
Ordinary Differential Equations I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations I 1 / 32 Theme of Last Few
More informationOrdinary Differential Equations I
Ordinary Differential Equations I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations I 1 / 27 Theme of Last Three
More informationAPPM 2360: Midterm exam 1 February 15, 2017
APPM 36: Midterm exam 1 February 15, 17 On the front of your Bluebook write: (1) your name, () your instructor s name, (3) your recitation section number and () a grading table. Text books, class notes,
More informationOrdinary Differential Equations I
Ordinary Differential Equations I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Doug James (and Justin Solomon) CS 205A: Mathematical Methods Ordinary Differential Equations I 1 / 35
More informationThe 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 informationChapter 6 - Ordinary Differential Equations
Chapter 6 - Ordinary Differential Equations 7.1 Solving Initial-Value Problems In this chapter, we will be interested in the solution of ordinary differential equations. Ordinary differential equations
More informationMath 322. Spring 2015 Review Problems for Midterm 2
Linear Algebra: Topic: Linear Independence of vectors. Question. Math 3. Spring Review Problems for Midterm Explain why if A is not square, then either the row vectors or the column vectors of A are linearly
More informationOrdinary Differential Equations
Ordinary Differential Equations Introduction: first order ODE We are given a function f(t,y) which describes a direction field in the (t,y) plane an initial point (t 0,y 0 ) We want to find a function
More informationSection 7.4 Runge-Kutta Methods
Section 7.4 Runge-Kutta Methods Key terms: Taylor methods Taylor series Runge-Kutta; methods linear combinations of function values at intermediate points Alternatives to second order Taylor methods Fourth
More informationThe family of Runge Kutta methods with two intermediate evaluations is defined by
AM 205: lecture 13 Last time: Numerical solution of ordinary differential equations Today: Additional ODE methods, boundary value problems Thursday s lecture will be given by Thomas Fai Assignment 3 will
More informationHigher-order differential equations
Higher-order differential equations Peyam Tabrizian Wednesday, November 16th, 2011 This handout is meant to give you a couple more example of all the techniques discussed in chapter 6, to counterbalance
More informationNonhomogeneous Equations and Variation of Parameters
Nonhomogeneous Equations Variation of Parameters June 17, 2016 1 Nonhomogeneous Equations 1.1 Review of First Order Equations If we look at a first order homogeneous constant coefficient ordinary differential
More informationNumerical Methods for the Solution of Differential Equations
Numerical Methods for the Solution of Differential Equations Markus Grasmair Vienna, winter term 2011 2012 Analytical Solutions of Ordinary Differential Equations 1. Find the general solution of the differential
More informationMath 23: Differential Equations (Winter 2017) Midterm Exam Solutions
Math 3: Differential Equations (Winter 017) Midterm Exam Solutions 1. [0 points] or FALSE? You do not need to justify your answer. (a) [3 points] Critical points or equilibrium points for a first order
More informationAnswer Key b c d e. 14. b c d e. 15. a b c e. 16. a b c e. 17. a b c d. 18. a b c e. 19. a b d e. 20. a b c e. 21. a c d e. 22.
Math 20580 Answer Key 1 Your Name: Final Exam May 8, 2007 Instructor s name: Record your answers to the multiple choice problems by placing an through one letter for each problem on this answer sheet.
More information1. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal
. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal 3 9 matrix D such that A = P DP, for A =. 3 4 3 (a) P = 4, D =. 3 (b) P = 4, D =. (c) P = 4 8 4, D =. 3 (d) P
More informationThe 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 informationLecture 1. Scott Pauls 1 3/28/07. Dartmouth College. Math 23, Spring Scott Pauls. Administrivia. Today s material.
Lecture 1 1 1 Department of Mathematics Dartmouth College 3/28/07 Outline Course Overview http://www.math.dartmouth.edu/~m23s07 Matlab Ordinary differential equations Definition An ordinary differential
More informationMIDTERM 1 PRACTICE PROBLEM SOLUTIONS
MIDTERM 1 PRACTICE PROBLEM SOLUTIONS Problem 1. Give an example of: (a) an ODE of the form y (t) = f(y) such that all solutions with y(0) > 0 satisfy y(t) = +. lim t + (b) an ODE of the form y (t) = f(y)
More informationMath K (24564) - Homework Solutions 02
Math 39100 K (24564) - Homework Solutions 02 Ethan Akin Office: NAC 6/287 Phone: 650-5136 Email: ethanakin@earthlink.net Spring, 2018 Contents Reduction of Order, B & D Chapter 3, p. 174 Constant Coefficient
More informationMath 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 informationSECTION 3.3. PROBLEM 22. The null space of a matrix A is: N(A) = {X : AX = 0}. Here are the calculations of AX for X = a,b,c,d, and e. =
SECTION 3.3. PROBLEM. The null space of a matrix A is: N(A) {X : AX }. Here are the calculations of AX for X a,b,c,d, and e. Aa [ ][ ] 3 3 [ ][ ] Ac 3 3 [ ] 3 3 [ ] 4+4 6+6 Ae [ ], Ab [ ][ ] 3 3 3 [ ]
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationMathQuest: Differential Equations
MathQuest: Differential Equations Solutions to Linear Systems. Consider the linear system given by dy dt = 4 True or False: Y e t t = is a solution. c False, but I am not very confident Y.. Consider the
More information. For each initial condition y(0) = y 0, there exists a. unique solution. In fact, given any point (x, y), there is a unique curve through this point,
1.2. Direction Fields: Graphical Representation of the ODE and its Solution Section Objective(s): Constructing Direction Fields. Interpreting Direction Fields. Definition 1.2.1. A first order ODE of the
More informationSolutions to Homework 5 - Math 3410
Solutions to Homework 5 - Math 34 (Page 57: # 489) Determine whether the following vectors in R 4 are linearly dependent or independent: (a) (, 2, 3, ), (3, 7,, 2), (, 3, 7, 4) Solution From x(, 2, 3,
More informationMath 341 Fall 2008 Friday December 12
FINAL EXAM: Differential Equations Math 341 Fall 2008 Friday December 12 c 2008 Ron Buckmire 1:00pm-4:00pm Name: Directions: Read all problems first before answering any of them. There are 17 pages in
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationSolution: In standard form (i.e. y + P (t)y = Q(t)) we have y t y = cos(t)
Math 380 Practice Final Solutions This is longer than the actual exam, which will be 8 to 0 questions (some might be multiple choice). You are allowed up to two sheets of notes (both sides) and a calculator,
More informationWe begin exploring Euler s method by looking at direction fields. Consider the direction field below.
Emma Reid- MA 66, Lesson 9 (SU 17) Euler s Method (.7) So far in this course, we have seen some very special types of first order ODEs. We ve seen methods to solve linear, separable, homogeneous, Bernoulli,
More informationAPPM 2360 Section Exam 3 Wednesday November 19, 7:00pm 8:30pm, 2014
APPM 2360 Section Exam 3 Wednesday November 9, 7:00pm 8:30pm, 204 ON THE FRONT OF YOUR BLUEBOOK write: () your name, (2) your student ID number, (3) lecture section, (4) your instructor s name, and (5)
More informationSolutions to Homework 3
Solutions to Homework 3 Section 3.4, Repeated Roots; Reduction of Order Q 1). Find the general solution to 2y + y = 0. Answer: The charactertic equation : r 2 2r + 1 = 0, solving it we get r = 1 as a repeated
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 information1 Differential Equations
Reading [Simon], Chapter 24, p. 633-657. 1 Differential Equations 1.1 Definition and Examples A differential equation is an equation involving an unknown function (say y = y(t)) and one or more of its
More informationSecond Order Linear Equations
October 13, 2016 1 Second And Higher Order Linear Equations In first part of this chapter, we consider second order linear ordinary linear equations, i.e., a differential equation of the form L[y] = d
More informationSolutions 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 informationHomework Solutions: , plus Substitutions
Homework Solutions: 2.-2.2, plus Substitutions Section 2. I have not included any drawings/direction fields. We can see them using Maple or by hand, so we ll be focusing on getting the analytic solutions
More informationMath 323 Exam 2 Sample Problems Solution Guide October 31, 2013
Math Exam Sample Problems Solution Guide October, Note that the following provides a guide to the solutions on the sample problems, but in some cases the complete solution would require more work or justification
More informationOld Math 330 Exams. David M. McClendon. Department of Mathematics Ferris State University
Old Math 330 Exams David M. McClendon Department of Mathematics Ferris State University Last updated to include exams from Fall 07 Contents Contents General information about these exams 3 Exams from Fall
More informationy0 = 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 informationPart IB Numerical Analysis
Part IB Numerical Analysis Definitions Based on lectures by G. Moore Notes taken by Dexter Chua Lent 206 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after
More informationFirst Order Linear Ordinary Differential Equations
First Order Linear Ordinary Differential Equations The most general first order linear ODE is an equation of the form p t dy dt q t y t f t. 1 Herepqarecalledcoefficients f is referred to as the forcing
More informationSection 2.4 Linear Equations
Section 2.4 Linear Equations Key Terms: Linear equation Homogeneous linear equation Nonhomogeneous (inhomogeneous) linear equation Integrating factor General solution Variation of parameters A first-order
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 informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationUniversity of Alberta ENGM 541: Modeling and Simulation of Engineering Systems Laboratory #5
University of Alberta ENGM 54: Modeling and Simulation of Engineering Systems Laboratory #5 M.G. Lipsett, Updated 00 Integration Methods with Higher-Order Truncation Errors with MATLAB MATLAB is capable
More informationAPPM 2360: Midterm exam 3 April 19, 2017
APPM 36: Midterm exam 3 April 19, 17 On the front of your Bluebook write: (1) your name, () your instructor s name, (3) your lecture section number and (4) a grading table. Text books, class notes, cell
More informationSection 9.8 Higher Order Linear Equations
Section 9.8 Higher Order Linear Equations Key Terms: Higher order linear equations Equivalent linear systems for higher order equations Companion matrix Characteristic polynomial and equation A linear
More informationMATH 307: Problem Set #3 Solutions
: Problem Set #3 Solutions Due on: May 3, 2015 Problem 1 Autonomous Equations Recall that an equilibrium solution of an autonomous equation is called stable if solutions lying on both sides of it tend
More informationSimple ODE Solvers - Derivation
Simple ODE Solvers - Derivation These notes provide derivations of some simple algorithms for generating, numerically, approximate solutions to the initial value problem y (t =f ( t, y(t y(t 0 =y 0 Here
More informationChapter 5 Exercises. (a) Determine the best possible Lipschitz constant for this function over 2 u <. u (t) = log(u(t)), u(0) = 2.
Chapter 5 Exercises From: Finite Difference Methods for Ordinary and Partial Differential Equations by R. J. LeVeque, SIAM, 2007. http://www.amath.washington.edu/ rjl/fdmbook Exercise 5. (Uniqueness for
More informationINTRODUCTION TO COMPUTER METHODS FOR O.D.E.
INTRODUCTION TO COMPUTER METHODS FOR O.D.E. 0. Introduction. The goal of this handout is to introduce some of the ideas behind the basic computer algorithms to approximate solutions to differential equations.
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationLinear algebra and differential equations (Math 54): Lecture 20
Linear algebra and differential equations (Math 54): Lecture 20 Vivek Shende April 7, 2016 Hello and welcome to class! Last time We started discussing differential equations. We found a complete set of
More information1 Systems of Differential Equations
March, 20 7- Systems of Differential Equations Let U e an open suset of R n, I e an open interval in R and : I R n R n e a function from I R n to R n The equation ẋ = ft, x is called a first order ordinary
More informationChap. 20: Initial-Value Problems
Chap. 20: Initial-Value Problems Ordinary Differential Equations Goal: to solve differential equations of the form: dy dt f t, y The methods in this chapter are all one-step methods and have the general
More informationA: Brief Review of Ordinary Differential Equations
A: Brief Review of Ordinary Differential Equations Because of Principle # 1 mentioned in the Opening Remarks section, you should review your notes from your ordinary differential equations (odes) course
More informationNumerical solution of ODEs
Numerical solution of ODEs Arne Morten Kvarving Department of Mathematical Sciences Norwegian University of Science and Technology November 5 2007 Problem and solution strategy We want to find an approximation
More informationRicharson Extrapolation for Runge-Kutta Methods
Richarson Extrapolation for Runge-Kutta Methods Zahari Zlatevᵃ, Ivan Dimovᵇ and Krassimir Georgievᵇ ᵃ Department of Environmental Science, Aarhus University, Frederiksborgvej 399, P. O. 358, 4000 Roskilde,
More informationOrdinary Differential Equations. Monday, October 10, 11
Ordinary Differential Equations Monday, October 10, 11 Problems involving ODEs can always be reduced to a set of first order differential equations. For example, By introducing a new variable z, this can
More informationACCELERATION OF RUNGE-KUTTA INTEGRATION SCHEMES
ACCELERATION OF RUNGE-KUTTA INTEGRATION SCHEMES PHAILAUNG PHOHOMSIRI AND FIRDAUS E. UDWADIA Received 24 November 2003 A simple accelerated third-order Runge-Kutta-type, fixed time step, integration scheme
More informationSect2.1. Any linear equation:
Sect2.1. Any linear equation: Divide a 0 (t) on both sides a 0 (t) dt +a 1(t)y = g(t) dt + a 1(t) a 0 (t) y = g(t) a 0 (t) Choose p(t) = a 1(t) a 0 (t) Rewrite it in standard form and ḡ(t) = g(t) a 0 (t)
More informationHomogeneous Equations with Constant Coefficients
Homogeneous Equations with Constant Coefficients MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 General Second Order ODE Second order ODEs have the form
More informationSOLVING DIFFERENTIAL EQUATIONS. Amir Asif. Department of Computer Science and Engineering York University, Toronto, ON M3J1P3
SOLVING DIFFERENTIAL EQUATIONS Amir Asif Department of Computer Science and Engineering York University, Toronto, ON M3J1P3 ABSTRACT This article reviews a direct method for solving linear, constant-coefficient
More informationPhysics 584 Computational Methods
Physics 584 Computational Methods Introduction to Matlab and Numerical Solutions to Ordinary Differential Equations Ryan Ogliore April 18 th, 2016 Lecture Outline Introduction to Matlab Numerical Solutions
More informationChapter 1: Introduction
Chapter 1: Introduction Definition: A differential equation is an equation involving the derivative of a function. If the function depends on a single variable, then only ordinary derivatives appear and
More information8.1 Introduction. Consider the initial value problem (IVP):
8.1 Introduction Consider the initial value problem (IVP): y dy dt = f(t, y), y(t 0)=y 0, t 0 t T. Geometrically: solutions are a one parameter family of curves y = y(t) in(t, y)-plane. Assume solution
More informationNumerical Integration of Ordinary Differential Equations for Initial Value Problems
Numerical Integration of Ordinary Differential Equations for Initial Value Problems Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@me.pdx.edu These slides are a
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 33 Almost Done! Last
More informationFirst-Order ODEs. Chapter Separable Equations. We consider in this chapter differential equations of the form dy (1.1)
Chapter 1 First-Order ODEs We consider in this chapter differential equations of the form dy (1.1) = F (t, y), where F (t, y) is a known smooth function. We wish to solve for y(t). Equation (1.1) is called
More informationFINAL EXAM SOLUTIONS, MATH 123
FINAL EXAM SOLUTIONS, MATH 23. Find the eigenvalues of the matrix ( 9 4 3 ) So λ = or 6. = λ 9 4 3 λ = ( λ)( 3 λ) + 36 = λ 2 7λ + 6 = (λ 6)(λ ) 2. Compute the matrix inverse: ( ) 3 3 = 3 4 ( 4/3 ) 3. Let
More informationOrdinary differential equation II
Ordinary Differential Equations ISC-5315 1 Ordinary differential equation II 1 Some Basic Methods 1.1 Backward Euler method (implicit method) The algorithm looks like this: y n = y n 1 + hf n (1) In contrast
More informationFINAL 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 information3.5 Undetermined Coefficients
3.5. UNDETERMINED COEFFICIENTS 153 11. t 2 y + ty + 4y = 0, y(1) = 3, y (1) = 4 12. t 2 y 4ty + 6y = 0, y(0) = 1, y (0) = 1 3.5 Undetermined Coefficients In this section and the next we consider the nonhomogeneous
More informationOrdinary Differential Equations
CHAPTER 8 Ordinary Differential Equations 8.1. Introduction My section 8.1 will cover the material in sections 8.1 and 8.2 in the book. Read the book sections on your own. I don t like the order of things
More informationExam 2 Solutions. (a) Is W closed under addition? Why or why not? W is not closed under addition. For example,
Exam 2 Solutions. Let V be the set of pairs of real numbers (x, y). Define the following operations on V : (x, y) (x, y ) = (x + x, xx + yy ) r (x, y) = (rx, y) Check if V together with and satisfy properties
More informationMA 527 first midterm review problems Hopefully final version as of October 2nd
MA 57 first midterm review problems Hopefully final version as of October nd The first midterm will be on Wednesday, October 4th, from 8 to 9 pm, in MTHW 0. It will cover all the material from the classes
More informationUC Berkeley Math 10B, Spring 2015: Midterm 2 Prof. Sturmfels, April 9, SOLUTIONS
UC Berkeley Math 10B, Spring 2015: Midterm 2 Prof. Sturmfels, April 9, SOLUTIONS 1. (5 points) You are a pollster for the 2016 presidential elections. You ask 0 random people whether they would vote for
More informationLecture 2. Classification of Differential Equations and Method of Integrating Factors
Math 245 - Mathematics of Physics and Engineering I Lecture 2. Classification of Differential Equations and Method of Integrating Factors January 11, 2012 Konstantin Zuev (USC) Math 245, Lecture 2 January
More informationThe integrating factor method (Sect. 1.1)
The integrating factor method (Sect. 1.1) Overview of differential equations. Linear Ordinary Differential Equations. The integrating factor method. Constant coefficients. The Initial Value Problem. Overview
More informationFinal Exam Sample Problems, Math 246, Spring 2018
Final Exam Sample Problems, Math 246, Spring 2018 1) Consider the differential equation dy dt = 9 y2 )y 2. a) Find all of its stationary points and classify their stability. b) Sketch its phase-line portrait
More informationChapter 9b: Numerical Methods for Calculus and Differential Equations. Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers
Chapter 9b: Numerical Methods for Calculus and Differential Equations Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers Acceleration Initial-Value Problems Consider a skydiver
More informationSystems 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 informationODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class
ODE Homework Due Wed. 9 August 2009; At the beginning of the class. (a) Solve Lẏ + Ry = E sin(ωt) with y(0) = k () L, R, E, ω are positive constants. (b) What is the limit of the solution as ω 0? (c) Is
More informationOrdinary Differential Equations: Initial Value problems (IVP)
Chapter Ordinary Differential Equations: Initial Value problems (IVP) Many engineering applications can be modeled as differential equations (DE) In this book, our emphasis is about how to use computer
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 29 Almost Done! No
More informationInitial value problems for ordinary differential equations
AMSC/CMSC 660 Scientific Computing I Fall 2008 UNIT 5: Numerical Solution of Ordinary Differential Equations Part 1 Dianne P. O Leary c 2008 The Plan Initial value problems (ivps) for ordinary differential
More informationReduction to the associated homogeneous system via a particular solution
June PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 5) Linear Algebra This study guide describes briefly the course materials to be covered in MA 5. In order to be qualified for the credit, one
More informationCS520: numerical ODEs (Ch.2)
.. CS520: numerical ODEs (Ch.2) Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca people.cs.ubc.ca/ ascher/520.html Uri Ascher (UBC) CPSC 520: ODEs (Ch. 2) Fall
More information8. Qualitative analysis of autonomous equations on the line/population dynamics models, phase line, and stability of equilibrium points (corresponds
c Dr Igor Zelenko, Spring 2017 1 8. Qualitative analysis of autonomous equations on the line/population dynamics models, phase line, and stability of equilibrium points (corresponds to section 2.5) 1.
More informationSection 3: Summary. Section 4. Evaluating Convolution Integrals. Convolution integral (memorise this):
Section 3: Summary Convolution integral (memorise this): f(t) input g(t) impulse response y(t) output y(t) g(t τ) f(τ) dτ Way to find the output of a linear system, described by a differential equation,
More informationThis is a closed book exam. No notes or calculators are permitted. We will drop your lowest scoring question for you.
Math 54 Fall 2017 Practice Final Exam Exam date: 12/14/17 Time Limit: 170 Minutes Name: Student ID: GSI or Section: This exam contains 9 pages (including this cover page) and 10 problems. Problems are
More informationHigher Order Taylor Methods
Higher Order Taylor Methods Marcelo Julio Alvisio & Lisa Marie Danz May 6, 2007 Introduction Differential equations are one of the building blocks in science or engineering. Scientists aim to obtain numerical
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 information40h + 15c = c = h
Chapter One Linear Systems I Solving Linear Systems Systems of linear equations are common in science and mathematics. These two examples from high school science [Onan] give a sense of how they arise.
More informationHigher Order Linear Equations
C H A P T E R 4 Higher Order Linear Equations 4.1 1. The differential equation is in standard form. Its coefficients, as well as the function g(t) = t, are continuous everywhere. Hence solutions are valid
More informationIntroduction to numerical simulations for Stochastic ODEs
Introduction to numerical simulations for Stochastic ODEs Xingye Kan Illinois Institute of Technology Department of Applied Mathematics Chicago, IL 60616 August 9, 2010 Outline 1 Preliminaries 2 Numerical
More informationODEs Cathal Ormond 1
ODEs Cathal Ormond 2 1. Separable ODEs Contents 2. First Order ODEs 3. Linear ODEs 4. 5. 6. Chapter 1 Separable ODEs 1.1 Definition: An ODE An Ordinary Differential Equation (an ODE) is an equation whose
More informationSolutions 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: This is an isocline associated with a slope
More informationDETERMINANTS DEFINED BY ROW OPERATIONS
DETERMINANTS DEFINED BY ROW OPERATIONS TERRY A. LORING. DETERMINANTS DEFINED BY ROW OPERATIONS Determinants of square matrices are best understood in terms of row operations, in my opinion. Most books
More information