2 Ordinary Differential Equations: Initial Value Problems
|
|
- Beatrix Cory West
- 5 years ago
- Views:
Transcription
1 Ordinar Differential Equations: Initial Value Problems Read sections 9., (9. for information), 9.3, 9.3., 9.3. (up to p. 396), Review questions 9.3, 9.4, 9.8, 9.9, Two Examples.. Foxes and Rabbits Assume that we have a forest where foxes and rabbits are living. The forest contains infinite resources of feed for the rabbits. So rabbits could live as the would if the would not meet foxes. As soon as a rabbit meets a fox the former will be eaten b the fox. On the other hand, rabbits are the onl feed for foxes. If a fox does not meet a rabbit sufficientl often, it will die. Let t be the number of rabbits at a certain time t. Similarl, let t be the number of foxes. Without the presence of foxes the number of rabbits will increase proportional to the number of individuals. The probabilit of a meeting between foxes and rabbits is proportional to the product of both populations. This gives rise to the following differential equation for the number of rabbits: d dt α rabbits reproduction β rabbits death Similarl, the differential equation for the number of foxes reads d dt γ foxes death δ foxes survival Such a sstem is called a Lotka-Volterra equation. This sstem was introduced in 95. In order to describe the two populations completel, we need to know the number of individuals at a certain point in time, sa t. This additional information is called initial values:.. Van der Pol Equation In the 9-ies van der Pol investigated a nonlinear amplifier circuit. The mathematical model he proposed is a differential equation of the form d dt ε The parameter ε describes some electrical properties of the amplifier circuit. Again, one needs initial values in order to single out a unique solution. Because the differential equation is of second order, two initial conditions are needed: α d dt β
2 . Generalities About Ordinar Differential Equations Differential equation are often obtained if one describes a continuous process. Most often the state of a sstem is described b a function of time t: t. A differential equation is a mathematical model which relates the values of a function and its derivatives and which must hold at ever point in time t. The notion ordinar indicates that there are onl derivatives with respect to one independent variable involved. The other independent variables are treated as parameters. The order of a differential equation is the order of the highest derivative which is present in the equation. The Lotka-Volterra equation is an example of a sstem of first order while the van der Pol equation is a second order equation. Assume that we are given a (scalar) differential equation of k-th order k f t k We define k new dependent variables u t t u t t u k. t k t The original k-th order differential equation is then equivalent to the following first-order sstem with k equations: u u. uk uk u u 3. u k f t u u k Hence, it is sufficient to consider onl sstems of first-order ordinar differential equations In more detail, this notation means t f t t t t d i f i t n t t i n dt i t i i n Note that this transformation is rather common. Almost all available codes in program libraries are designed to solve sstems of first-order ordinar differential equations. For instance, this is the approach taken in MATLAB.
3 .3 One-Step Methods In order to simplif the notation we will assume for the moment that we want to solve one scalar differential equation, onl, that means n. It is almost obvious to see how the considerations and methods can be generalized to sstems of equations. The differential equation at hand is the mathematical model of a process which the scientist or engineer is interested in. Depending on the application he might be interested in asking different questions: The engineer might be interested in knowing how the solution curve looks like in detail. The mathematical question would be: How does the solution behave on a given finite interval t T? Another engineer might be interested in a more qualitative picture, for example if the intended working point is reached if the sstem is switched on. The mathematical question would be: How does the solution behave on an infinite interval t? Of course, one cannot compute over infinite time intervals, but one wants to obtain a solution over a large interval such that the solution behavior in a vicinit of infinit becomes clear. Usuall, the mathematical model contains a lot of parameters, most of them with uncertainties. The engineer might wonder how sensitive the solution is to changes in those parameters. The mathematical question would be: What is the conditioning of the differential equation? Some states of the sstem ma be of special interest. In our introductor example the question might be when all animals are dead. The mathematical question would be: Does there exist a point in time where the solution has special well-defined properties? All these questions require special algorithms and specialized programs. The basic question, however, is if the ordinar differential equation can be solved reliabl and efficientl. There are onl ver few problems where an analtic (or, closed form) solution can be provided. Therefore, we are almost alwas required to use numerical methods. Most of the numerical methods for solving initial value problems for ordinar differential equations are based on a discretization method which is called the finite difference method. In contrast to an analtical solution, we do no longer attempt to provide an approximation of the solution for all points t t T. Instead we tr to compute onl a finite number (sufficientl man, of course) of approximations k t k to the true values t k for selected points t t t t k T such that t t t t k T Such a finite sequence of points is called a grid. The general idea consists of the following process: At t t, the initial value t is known. Then we tr to compute a value t b using the function f, the previous value, and the grid points t t. Our problem is now reduced to finding a solution to the problem t f t t t 3
4 subject to the initial condition t Obviousl, this problem can be handled in the same wa as before thus providing an approximation t at t. This process can be continued until T is reached. More general, if we know k, then we can compute k. The difference between two consecutive grid points is called the step size h k, h k t k t k In practical computations, one uses varing (or adaptive) step sizes. For our purposes it is sufficient to consider onl constant step sizes h k h such that.3. Euler s Method We start b writing down the Talor expansion for the true solution: t k t kh where s k t k t k t k t k t k h t k h f t k t k h s k h s k. This gives rise to the following recursive method: given k k h f t k for k Example.. We consider the following simple initial value problem for t The exact solution of this problem is t e t. Let us choose a step size of.5. The -th approximation is the initial guess. The next approximation is obtained for t t h 5 5. According to Euler s method, this becomes h f 5 5 The exact solution is 5 e Let us take one additional step. The next grid point is t t h. Euler s approximation gives This gives the following little table: h f
5 t k k t k error It is relativel simple to implement Euler s method in a MATLAB function. The input arguments to such a function should be the data of the problem at hand: the function f, the interval t T, the initial value, and the step size h. The output should include the computed grid points t k as well as the approximations k. This gives rise to the following code: function [t,] = Euler(f,interval,,h) t = interval(); T = interval(); % and f are assumed to provide column vectors = ; t = t; k = ; tk = t; while tk+h < T+h/ k = k+h*feval(f,tk,k); = [,k]; tk = tk+h; t = [t,tk]; end This function can now be used for solving our little example problem. The result is given in Figure. clear clf f = inline(, t, ); [t,] = Euler(f,[,],,.5); % plot the result plot(t,,t,, o ) xlabel( t ) label( ) title( A simple example ) Euler s method is said to be a one-step method because the computation of the new value k at t k depends onl on information at the previous grid point t k. Such a method can be written down more generall as k k hφ t k k h 5
6 .5 A simple example t Figure : A simple example.3. Accurac The most important question is how accurate the numerical results is and if the accurac becomes better and better if the step size h becomes smaller and smaller. As usual, there are two tpes of error sources: discretization errors which appear because the differential equation is onl approximated; rounding errors. We will in the following consider onl the discretization error because, in the context of the solution of ordinar differential equations, the rounding error is usuall negligible. What we are mostl interested in is the global error e k k t k Unfortunatel, this error is not available explicitl. Even reliable and practicall useful estimates are hard to obtain. Opposed to that, it appears to be rather eas to describe the local error. The latter is the error which arises if one takes onl one step of the method with an errorfree k t k : This statement does not hold for other algorithms! l k k t k t k hφ t k t k h t k 6
7 Example.. For Euler s method, we can compute l k t k t k h s k h f t k t k t k h f t k t k t k h t k h s k Hence, we can write l k O h. A method has the order of accurac (or simpl the order) p if l k O h p This means that Euler s method has the order. Wh is the order of a method important? One can show that the order of magnitude of the global error for one-step methods can be estimated b that of the local error, that is, for a method of order p it holds e k O h p if the step size h is sufficientl small. Example.3. What are the consequences of the order of a method? The user of a method is usuall interested in obtaining a solution with a given accurac with minimal computation time. Therefore, it is probabl a good idea to use step sizes as large as possible. If a method provides a result with a certain accurac for a given step size h, a reduction of the step size b a factor a leads to a reduction of the error b a factor of a p if we have a p-th order method. This is illustrated in the following table where we have used the factor a which amounts to step size halving. if the order p is 3 4 the error is reduced b So the general rule is to alwas appl a higher order method if there are no other reasons for not doing it..3.3 Runge-Kutta methods There are a number of higher order one-step methods available. Most of them belong to the huge class of Runge-Kutta methods. As examples, we provide two methods of this tpe which were invented more than ears ago. 7
8 Heun s Method. This method can be written rather compact in the form k h f t k k f t k k h f t k k It is convenient to write down the method in a more sstematic wa: k k k f t k k f t k h k h k k Heun s method has second order. Here we give a sample implementation of Heun s method in MATLAB. Do not use it for real application problems! hk function [t,] = Heun(f,interval,,h) t = interval(); T = interval(); % and f are assumed to provide column vectors = ; t = t; k = ; tk = t; while tk+h < T+h/ k = feval(f,tk,k); k = feval(f,tk+h,k+h*k); k = k+.5*h*(k+k); = [,k]; tk = tk+h; t = [t,tk]; end Example.4. Let us demonstrate the application of our code for solving the Lotka-Volterra equation. We use the following parameters: α β γ δ T 9 The initial values are 3 and 5. Let us choose for example h (which is rather large!). Then the following codes provides us with nice results (tr it out b ourself). 8
9 clear clf [t,] = Heun(@flotka,[,9],[3;5],.); plot(t,(,:),t,(,:), -- ) xlabel( time ) label( number of individuals ) title( Lotka-Volterra equation ) % % The function: Should appear in a separate file!!! function p = flotka(t,) p = [*()-.*()*();... -()+.*()*()]; Classical Runge-Kutta Method. This is probabl the best-known Runge-Kutta method. It should be noted, however, that this method is superseded b much more efficient methods. The classical method can be formulated as follows: This method is of order 4. k f t k k k f h t k h k k k 3 f h t k h k k k 4 f t k h k hk 3 h k k k 3 k 4 6 k Runge-Kutta Methods in MATLAB. MATLAB provides two (explicit) Runge-Kutta methods. ode45 is a fifth-order method proposed b Dormand and Prince. The other solver is ode3 having order three using a method b Bogacki and Shampine. Both methods use adaptive step size. According to our considerations about the consequences of different orders on the computational work it is proposed alwas to tr the Dormand-Prince method first..3.4 Automatic Step Size Control While examining the local error of Euler s method we found out that beside its dependence on h it depends essentiall on the second derivative of the solution. If one requires a certain bound on the local error, sa l k tol then the step size can be chosen larger where the second derivative is small, and smaller where the second derivative is larger. Such an idea can be realized. The resulting algorithm is called adaptive. It must be noted, however, that the control variable is the local error and not the global 9
10 error. This automatic step size control can onl work reliabl if the problem at hand is wellconditioned. A discussion of that can be found in the course book..4 Some Numerical Examples Here we show some applications. We will use the van der Pol equation as a running example. This can be a ver hard problem depending on the parameter ε. In order to simplif things we will choose ε. This preserves nice features while the problem is not too hard. There are different was of presenting the solution of a problem to the user. The worst one is obviousl to provide a table containing all results. The useful information ma be too much hidden in the numbers. Therefore, a graphical representation is usuall the method of choice. Before presenting a graph one should think ver carefull how the intended informations can be shown as informativel (and suggestive?) as possible. In the context of ordinar differential equations one ver often chooses among the following two possibilities: Plot trajectories. This means that some or all solution components i t are plotted as a function of the independent variable t. Plot a phase portrait. Here one plots two solution components i t j t in a i j - coordinate sstem. The independent variable t appears merel as a parameter of the curve. Our method of choice is MATLAB s ode45 solver. Figure contains a trajector plot and a phase portrait. The initial values are and Periodic solution 4 Periodic solution t Figure : Solution to the van der Pol equation The graphs in Figure 3 contain some representations which demonstrate different possibilities of data enhancement. We will use phase portraits. The first plot gives the raw data as the are provided b ode45. The second plot uses piecewise linear interpolation which is the standard algorithm of MATLAB s plot command. The final plot uses advanced interpolation techniques
11 which are built-in into ode45. Note that the latter interpolation provides a highl accurate interpolation even if the data points are rather sparsel distributed!. 4 Phase portrait: Raw data 4 Phase portrait: Linear interpolation Phase portrait: ode45 interpolation Figure 3: Interpolation of discrete data In the last plot (Figure 4), we compared the approximate solution obtained b Euler s method with a step size of and, respectivel, with the solution obtained b ode45.
12 5 Phase portrait: Euler method 4 3 ode45 Euler. Euler Figure 4: Euler solutions of the van der Pol equation
Ordinary Differential Equations n
Numerical Analsis MTH63 Ordinar Differential Equations Introduction Talor Series Euler Method Runge-Kutta Method Predictor Corrector Method Introduction Man problems in science and engineering when formulated
More informationGraded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo
2008-11-07 Graded Project #1 Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo This homework is due to be handed in on Wednesday 12 November 2008 before 13:00 in the post box of the numerical
More informationOn Solving Systems of ODEs Numerically
On Solving Sstems of ODEs Numericall Temple H. Fa, Stephan V. Joubert, and Andrew Mkolesia Department of Mathematical Technolog Tshwane Universit of Technolog Private Bag x680 Pretoria 0001 SOUTH AFRICA
More informationCh 3 Alg 2 Note Sheet.doc 3.1 Graphing Systems of Equations
Ch 3 Alg Note Sheet.doc 3.1 Graphing Sstems of Equations Sstems of Linear Equations A sstem of equations is a set of two or more equations that use the same variables. If the graph of each equation =.4
More informationInitial Value Problems for. Ordinary Differential Equations
Initial Value Problems for Ordinar Differential Equations INTRODUCTION Equations which are composed of an unnown function and its derivatives are called differential equations. It becomes an initial value
More informationLinear Transformations
inear Transformations 6 The main objects of stud in an course in linear algebra are linear functions: Definition A function : V! W is linear if V and W are vector spaces and (ru + sv) r(u)+s(v) for all
More informationy R T However, the calculations are easier, when carried out using the polar set of co-ordinates ϕ,r. The relations between the co-ordinates are:
Curved beams. Introduction Curved beams also called arches were invented about ears ago. he purpose was to form such a structure that would transfer loads, mainl the dead weight, to the ground b the elements
More informationIntroduction to Differential Equations. National Chiao Tung University Chun-Jen Tsai 9/14/2011
Introduction to Differential Equations National Chiao Tung Universit Chun-Jen Tsai 9/14/011 Differential Equations Definition: An equation containing the derivatives of one or more dependent variables,
More informationDifferential Equations FMNN10 Graded Project #1 c G Söderlind 2017
Differential Equations FMNN10 Graded Project #1 c G Söderlind 2017 Published 2017-10-30. Instruction in computer lab 2017-11-02/08/09. Project report due date: Monday 2017-11-13 at 10:00. Goals. The goal
More informationIntroduction to Vector Spaces Linear Algebra, Spring 2011
Introduction to Vector Spaces Linear Algebra, Spring 2011 You probabl have heard the word vector before, perhaps in the contet of Calculus III or phsics. You probabl think of a vector like this: 5 3 or
More information10 Back to planar nonlinear systems
10 Back to planar nonlinear sstems 10.1 Near the equilibria Recall that I started talking about the Lotka Volterra model as a motivation to stud sstems of two first order autonomous equations of the form
More informationOrdinary Differential Equations
Instructor: José Antonio Agapito Ruiz UCSC Summer Session II, 00 Math Ordinar Differential Equations Final Solution This report does not contain full answers to all questions of the final. Instead, ou
More informationCONTINUOUS SPATIAL DATA ANALYSIS
CONTINUOUS SPATIAL DATA ANALSIS 1. Overview of Spatial Stochastic Processes The ke difference between continuous spatial data and point patterns is that there is now assumed to be a meaningful value, s
More informationComputational Methods CMSC/AMSC/MAPL 460. Ordinary differential equations
/0/0 Computational Methods CMSC/AMSC/MAPL 460 Ordinar differential equations Ramani Duraiswami Dept. of Computer Science Several slides adapted from Profs. Dianne O Lear and Eric Sandt TAMU Higher Order
More information3.0 PROBABILITY, RANDOM VARIABLES AND RANDOM PROCESSES
3.0 PROBABILITY, RANDOM VARIABLES AND RANDOM PROCESSES 3.1 Introduction In this chapter we will review the concepts of probabilit, rom variables rom processes. We begin b reviewing some of the definitions
More informationCHAPTER 80 NUMERICAL METHODS FOR FIRST-ORDER DIFFERENTIAL EQUATIONS
CHAPTER 8 NUMERICAL METHODS FOR FIRST-ORDER DIFFERENTIAL EQUATIONS EXERCISE 33 Page 834. Use Euler s method to obtain a numerical solution of the differential equation d d 3, with the initial conditions
More informationMath 2930 Worksheet Equilibria and Stability
Math 2930 Worksheet Equilibria and Stabilit Week 3 September 7, 2017 Question 1. (a) Let C be the temperature (in Fahrenheit) of a cup of coffee that is cooling off to room temperature. Which of the following
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 information11.4 Polar Coordinates
11. Polar Coordinates 917 11. Polar Coordinates In Section 1.1, we introduced the Cartesian coordinates of a point in the plane as a means of assigning ordered pairs of numbers to points in the plane.
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 informationScientific Computing with Case Studies SIAM Press, Lecture Notes for Unit V Solution of
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 1 Dianne P. O Leary c 2008 1 The
More informationA First Course on Kinetics and Reaction Engineering Supplemental Unit S5. Solving Initial Value Differential Equations
Supplemental Unit S5. Solving Initial Value Differential Equations Defining the Problem This supplemental unit describes how to solve a set of initial value ordinary differential equations (ODEs) numerically.
More informationLab 5: Nonlinear Systems
Lab 5: Nonlinear Systems Goals In this lab you will use the pplane6 program to study two nonlinear systems by direct numerical simulation. The first model, from population biology, displays interesting
More informationSolution: (a) Before opening the parachute, the differential equation is given by: dv dt. = v. v(0) = 0
Math 2250 Lab 4 Name/Unid: 1. (25 points) A man bails out of an airplane at the altitute of 12,000 ft, falls freely for 20 s, then opens his parachute. Assuming linear air resistance ρv ft/s 2, taking
More informationEngineering Mathematics I
Engineering Mathematics I_ 017 Engineering Mathematics I 1. Introduction to Differential Equations Dr. Rami Zakaria Terminolog Differential Equation Ordinar Differential Equations Partial Differential
More informationSOLVING SOME PHYSICAL PROBLEMS USING THE METHODS OF NUMERICAL MATHEMATICS WITH THE HELP OF SYSTEM MATHEMATICA
SOLVING SOME PHYSICAL PROBLEMS USING THE METHODS OF NUMERICAL MATHEMATICS WITH THE HELP OF SYSTEM MATHEMATICA Pavel Trojovský Jaroslav Seibert Antonín Slabý ABSTRACT Ever future teacher of mathematics
More informationLESSON #12 - FORMS OF A LINE COMMON CORE ALGEBRA II
LESSON # - FORMS OF A LINE COMMON CORE ALGEBRA II Linear functions come in a variet of forms. The two shown below have been introduced in Common Core Algebra I and Common Core Geometr. TWO COMMON FORMS
More informationSample Solutions for Assignment 3.
AMath 383, Autumn Sample Solutions for Assignment 3. Reading: Chs. 4-5.. Eercise 7 of Chapter 3. If $X is invested toda at 3% interest compounded continuousl, then in ears it will be worth Xe (.3 ) = Xe.6.
More informationSection 8.5 Parametric Equations
504 Chapter 8 Section 8.5 Parametric Equations Man shapes, even ones as simple as circles, cannot be represented as an equation where is a function of. Consider, for eample, the path a moon follows as
More information2.2 Equations of Lines
660_ch0pp07668.qd 10/16/08 4:1 PM Page 96 96 CHAPTER Linear Functions and Equations. Equations of Lines Write the point-slope and slope-intercept forms Find the intercepts of a line Write equations for
More informationExperimental Uncertainty Review. Abstract. References. Measurement Uncertainties and Uncertainty Propagation
Experimental Uncertaint Review Abstract This is intended as a brief summar of the basic elements of uncertaint analsis, and a hand reference for laborator use. It provides some elementar "rules-of-thumb"
More informationINTRODUCTION TO DIFFERENTIAL EQUATIONS
INTRODUCTION TO DIFFERENTIAL EQUATIONS. Definitions and Terminolog. Initial-Value Problems.3 Differential Equations as Mathematical Models CHAPTER IN REVIEW The words differential and equations certainl
More informationMethods of Solving Ordinary Differential Equations (Online)
7in 0in Felder c0_online.te V3 - Januar, 05 0:5 A.M. Page CHAPTER 0 Methods of Solving Ordinar Differential Equations (Online) 0.3 Phase Portraits Just as a slope field (Section.4) gives us a wa to visualize
More information12.1 Systems of Linear equations: Substitution and Elimination
. Sstems of Linear equations: Substitution and Elimination Sstems of two linear equations in two variables A sstem of equations is a collection of two or more equations. A solution of a sstem in two variables
More informationSolutions to Two Interesting Problems
Solutions to Two Interesting Problems Save the Lemming On each square of an n n chessboard is an arrow pointing to one of its eight neighbors (or off the board, if it s an edge square). However, arrows
More informationWILEY. Differential Equations with MATLAB (Third Edition) Brian R. Hunt Ronald L. Lipsman John E. Osborn Jonathan M. Rosenberg
Differential Equations with MATLAB (Third Edition) Updated for MATLAB 2011b (7.13), Simulink 7.8, and Symbolic Math Toolbox 5.7 Brian R. Hunt Ronald L. Lipsman John E. Osborn Jonathan M. Rosenberg All
More informationApproximation of Continuous Functions
Approximation of Continuous Functions Francis J. Narcowich October 2014 1 Modulus of Continuit Recall that ever function continuous on a closed interval < a x b < is uniforml continuous: For ever ɛ > 0,
More informationMathematics 309 Conic sections and their applicationsn. Chapter 2. Quadric figures. ai,j x i x j + b i x i + c =0. 1. Coordinate changes
Mathematics 309 Conic sections and their applicationsn Chapter 2. Quadric figures In this chapter want to outline quickl how to decide what figure associated in 2D and 3D to quadratic equations look like.
More informationQUADRATIC AND CONVEX MINIMAX CLASSIFICATION PROBLEMS
Journal of the Operations Research Societ of Japan 008, Vol. 51, No., 191-01 QUADRATIC AND CONVEX MINIMAX CLASSIFICATION PROBLEMS Tomonari Kitahara Shinji Mizuno Kazuhide Nakata Toko Institute of Technolog
More informationQUADRATIC FUNCTION REVIEW
Name: Date: QUADRATIC FUNCTION REVIEW Linear and eponential functions are used throughout mathematics and science due to their simplicit and applicabilit. Quadratic functions comprise another ver important
More information6.4 graphs OF logarithmic FUnCTIOnS
SECTION 6. graphs of logarithmic functions 9 9 learning ObjeCTIveS In this section, ou will: Identif the domain of a logarithmic function. Graph logarithmic functions. 6. graphs OF logarithmic FUnCTIOnS
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for solving higher order ODEs and
More information14.1 Systems of Linear Equations in Two Variables
86 Chapter 1 Sstems of Equations and Matrices 1.1 Sstems of Linear Equations in Two Variables Use the method of substitution to solve sstems of equations in two variables. Use the method of elimination
More informationMath 107: Calculus II, Spring 2015: Midterm Exam II Monday, April Give your name, TA and section number:
Math 7: Calculus II, Spring 25: Midterm Exam II Monda, April 3 25 Give our name, TA and section number: Name: TA: Section number:. There are 5 questions for a total of points. The value of each part of
More informationExponential and Logarithmic Functions
Eponential and Logarithmic Functions 6 Figure Electron micrograph of E. Coli bacteria (credit: Mattosaurus, Wikimedia Commons) CHAPTER OUTLINE 6. Eponential Functions 6. Logarithmic Properties 6. Graphs
More informationLESSON #48 - INTEGER EXPONENTS COMMON CORE ALGEBRA II
LESSON #8 - INTEGER EXPONENTS COMMON CORE ALGEBRA II We just finished our review of linear functions. Linear functions are those that grow b equal differences for equal intervals. In this unit we will
More information8. BOOLEAN ALGEBRAS x x
8. BOOLEAN ALGEBRAS 8.1. Definition of a Boolean Algebra There are man sstems of interest to computing scientists that have a common underling structure. It makes sense to describe such a mathematical
More informationMethods for Advanced Mathematics (C3) Coursework Numerical Methods
Woodhouse College 0 Page Introduction... 3 Terminolog... 3 Activit... 4 Wh use numerical methods?... Change of sign... Activit... 6 Interval Bisection... 7 Decimal Search... 8 Coursework Requirements on
More informationPartial differentiation. Background mathematics review
Partial differentiation Background mathematics review David Miller Partial differentiation Partial derivatives Background mathematics review David Miller Partial derivative Suppose we have a function f(,)
More informationLESSON #11 - FORMS OF A LINE COMMON CORE ALGEBRA II
LESSON # - FORMS OF A LINE COMMON CORE ALGEBRA II Linear functions come in a variet of forms. The two shown below have been introduced in Common Core Algebra I and Common Core Geometr. TWO COMMON FORMS
More informationMidterm 2. MAT17C, Spring 2014, Walcott
Mierm 2 MAT7C, Spring 24, Walcott Note: There are SEVEN questions. Points for each question are shown in a table below and also in parentheses just after each question (and part). Be sure to budget our
More informationEquations of lines in
Roberto s Notes on Linear Algebra Chapter 6: Lines, planes an other straight objects Section 1 Equations of lines in What ou nee to know alrea: The ot prouct. The corresponence between equations an graphs.
More information8.1 Exponents and Roots
Section 8. Eponents and Roots 75 8. Eponents and Roots Before defining the net famil of functions, the eponential functions, we will need to discuss eponent notation in detail. As we shall see, eponents
More informationYou don't have to be a mathematician to have a feel for numbers. John Forbes Nash, Jr.
Course Title: Real Analsis Course Code: MTH3 Course instructor: Dr. Atiq ur Rehman Class: MSc-II Course URL: www.mathcit.org/atiq/fa5-mth3 You don't have to be a mathematician to have a feel for numbers.
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MATLAB sessions: Laboratory 4 MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for
More informationEigenvectors and Eigenvalues 1
Ma 2015 page 1 Eigenvectors and Eigenvalues 1 In this handout, we will eplore eigenvectors and eigenvalues. We will begin with an eploration, then provide some direct eplanation and worked eamples, and
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol. I - Numerical Methods for Ordinary Differential Equations and Dynamic Systems - E.A.
COMPUTATIOAL METHODS AD ALGORITHMS Vol. I umerical Methods for Ordinar Differential Equations and Dnamic Sstems E.A. oviov UMERICAL METHODS FOR ORDIARY DIFFERETIAL EQUATIOS AD DYAMIC SYSTEMS E.A. oviov
More information4452 Mathematical Modeling Lecture 13: Chaos and Fractals
Math Modeling Lecture 13: Chaos and Fractals Page 1 442 Mathematical Modeling Lecture 13: Chaos and Fractals Introduction In our tetbook, the discussion on chaos and fractals covers less than 2 pages.
More information4.7. Newton s Method. Procedure for Newton s Method HISTORICAL BIOGRAPHY
4. Newton s Method 99 4. Newton s Method HISTORICAL BIOGRAPHY Niels Henrik Abel (18 189) One of the basic problems of mathematics is solving equations. Using the quadratic root formula, we know how to
More information(Refer Slide Time: 00:32)
Nonlinear Dynamical Systems Prof. Madhu. N. Belur and Prof. Harish. K. Pillai Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture - 12 Scilab simulation of Lotka Volterra
More informationSpotlight on the Extended Method of Frobenius
113 Spotlight on the Extended Method of Frobenius See Sections 11.1 and 11.2 for the model of an aging spring. Reference: Section 11.4 and SPOTLIGHT ON BESSEL FUNCTIONS. Bessel functions of the first kind
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 75 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP. Use MATLAB solvers for solving higher order ODEs and systems
More informationChapter 8 Notes SN AA U2C8
Chapter 8 Notes SN AA U2C8 Name Period Section 8-: Eploring Eponential Models Section 8-2: Properties of Eponential Functions In Chapter 7, we used properties of eponents to determine roots and some of
More information3.7 InveRSe FUnCTIOnS
CHAPTER functions learning ObjeCTIveS In this section, ou will: Verif inverse functions. Determine the domain and range of an inverse function, and restrict the domain of a function to make it one-to-one.
More informationMath Notes on sections 7.8,9.1, and 9.3. Derivation of a solution in the repeated roots case: 3 4 A = 1 1. x =e t : + e t w 2.
Math 7 Notes on sections 7.8,9., and 9.3. Derivation of a solution in the repeated roots case We consider the eample = A where 3 4 A = The onl eigenvalue is = ; and there is onl one linearl independent
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationy = (1 y)cos t Answer: The equation is not autonomous because of the cos t term.
Math 211 Homework #4 Februar 9, 2001 2.9.2. = 1 2 + 2 Answer: Note that = 1 2 + 2 is autonomous, having form = f(). Solve the equation f()= 0 to find the equilibrium points. f()= 0 1 2 + 2 = 0 = 1. Thus,
More informationSymmetry Arguments and the Role They Play in Using Gauss Law
Smmetr Arguments and the Role The la in Using Gauss Law K. M. Westerberg (9/2005) Smmetr plas a ver important role in science in general, and phsics in particular. Arguments based on smmetr can often simplif
More informationLesson 9: Predator-Prey and ode45
Lesson 9: Predator-Prey and ode45 9.1 Applied Problem. In this lesson we will allow for more than one population where they depend on each other. One population could be the predator such as a fox, and
More informationThis document is stored in Documents/4C/nablaoperator.tex. Compile it with LaTex. VECTOR CALCULUS
This document is stored in Documents/4C/nablaoperator.tex. Compile it with LaTex. November 12, 2014 Hans P. Paar VECTOR CALCULUS 1 Introduction Vector calculus is a branch of mathematics that allows differentiation
More informationFeedforward Neural Networks
Chapter 4 Feedforward Neural Networks 4. Motivation Let s start with our logistic regression model from before: P(k d) = softma k =k ( λ(k ) + w d λ(k, w) ). (4.) Recall that this model gives us a lot
More information6. Linear transformations. Consider the function. f : R 2 R 2 which sends (x, y) (x, y)
Consider the function 6 Linear transformations f : R 2 R 2 which sends (x, ) (, x) This is an example of a linear transformation Before we get into the definition of a linear transformation, let s investigate
More informationLESSON #24 - POWER FUNCTIONS COMMON CORE ALGEBRA II
1 LESSON #4 - POWER FUNCTIONS COMMON CORE ALGEBRA II Before we start to analze polnomials of degree higher than two (quadratics), we first will look at ver simple functions known as power functions. The
More informationSection B. Ordinary Differential Equations & its Applications Maths II
Section B Ordinar Differential Equations & its Applications Maths II Basic Concepts and Ideas: A differential equation (D.E.) is an equation involving an unknown function (or dependent variable) of one
More informationImplicit Second Derivative Hybrid Linear Multistep Method with Nested Predictors for Ordinary Differential Equations
American Scientiic Research Journal or Engineering, Technolog, and Sciences (ASRJETS) ISSN (Print) -44, ISSN (Online) -44 Global Societ o Scientiic Research and Researchers http://asretsournal.org/ Implicit
More informationWeek 3 September 5-7.
MA322 Weekl topics and quiz preparations Week 3 September 5-7. Topics These are alread partl covered in lectures. We collect the details for convenience.. Solutions of homogeneous equations AX =. 2. Using
More informationENGI 3424 Mid Term Test Solutions Page 1 of 9
ENGI 344 Mid Term Test 07 0 Solutions Page of 9. The displacement st of the free end of a mass spring sstem beond the equilibrium position is governed b the ordinar differential equation d s ds 3t 6 5s
More informationPre-AP Algebra 2 Lesson 1-1 Basics of Functions
Lesson 1-1 Basics of Functions Objectives: The students will be able to represent functions verball, numericall, smbolicall, and graphicall. The students will be able to determine if a relation is a function
More informationWe have examined power functions like f (x) = x 2. Interchanging x
CHAPTER 5 Eponential and Logarithmic Functions We have eamined power functions like f =. Interchanging and ields a different function f =. This new function is radicall different from a power function
More information520 Chapter 9. Nonlinear Differential Equations and Stability. dt =
5 Chapter 9. Nonlinear Differential Equations and Stabilit dt L dθ. g cos θ cos α Wh was the negative square root chosen in the last equation? (b) If T is the natural period of oscillation, derive the
More informationThe Maze Generation Problem is NP-complete
The Mae Generation Problem is NP-complete Mario Alviano Department of Mathematics, Universit of Calabria, 87030 Rende (CS), Ital alviano@mat.unical.it Abstract. The Mae Generation problem has been presented
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationLESSON #28 - POWER FUNCTIONS COMMON CORE ALGEBRA II
1 LESSON #8 - POWER FUNCTIONS COMMON CORE ALGEBRA II Before we start to analze polnomials of degree higher than two (quadratics), we first will look at ver simple functions known as power functions. The
More informationSecond-Order Linear Differential Equations C 2
C8 APPENDIX C Additional Topics in Differential Equations APPENDIX C. Second-Order Homogeneous Linear Equations Second-Order Linear Differential Equations Higher-Order Linear Differential Equations Application
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationLinear Programming. Maximize the function. P = Ax + By + C. subject to the constraints. a 1 x + b 1 y < c 1 a 2 x + b 2 y < c 2
Linear Programming Man real world problems require the optimization of some function subject to a collection of constraints. Note: Think of optimizing as maimizing or minimizing for MATH1010. For eample,
More informationS#ff ODEs and Systems of ODEs
S#ff ODEs and Systems of ODEs Popula#on Growth Modeling Let the number of individuals in a given area at time t be. At time the number is so that is the number of individuals that have arrived in the area
More informationFunctions. Introduction
Functions,00 P,000 00 0 70 7 80 8 0 000 00 00 Figure Standard and Poor s Inde with dividends reinvested (credit "bull": modification of work b Praitno Hadinata; credit "graph": modification of work b MeasuringWorth)
More informationIdentification of Nonlinear Dynamic Systems with Multiple Inputs and Single Output using discrete-time Volterra Type Equations
Identification of Nonlinear Dnamic Sstems with Multiple Inputs and Single Output using discrete-time Volterra Tpe Equations Thomas Treichl, Stefan Hofmann, Dierk Schröder Institute for Electrical Drive
More informationFIRST- AND SECOND-ORDER IVPS The problem given in (1) is also called an nth-order initial-value problem. For example, Solve: Solve:
.2 INITIAL-VALUE PROBLEMS 3.2 INITIAL-VALUE PROBLEMS REVIEW MATERIAL Normal form of a DE Solution of a DE Famil of solutions INTRODUCTION We are often interested in problems in which we seek a solution
More informationSection 3.1. ; X = (0, 1]. (i) f : R R R, f (x, y) = x y
Paul J. Bruillard MATH 0.970 Problem Set 6 An Introduction to Abstract Mathematics R. Bond and W. Keane Section 3.1: 3b,c,e,i, 4bd, 6, 9, 15, 16, 18c,e, 19a, 0, 1b Section 3.: 1f,i, e, 6, 1e,f,h, 13e,
More informationSolving ODEs and PDEs in MATLAB. Sören Boettcher
16.02.2009 Introduction Quick introduction to syntax ODE in the form of Initial Value Problems (IVP) what equations can handle how to code into how to choose the right solver how to get the solver to do
More informationADDITIVELY SEPARABLE FUNCTIONS
ADDITIVELY SEPARABLE FUNCTIONS FUNCTIONS OF THE FORM F (, ) =f()+g() STEVEN F. BELLENOT Abstract. A sample project similar to the project for Spring 25. The picture below has a hortizontal space coordinate
More informationRemark on the Sensitivity of Simulated Solutions of the Nonlinear Dynamical System to the Used Numerical Method
International Journal of Mathematical Analysis Vol. 9, 2015, no. 55, 2749-2754 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2015.59236 Remark on the Sensitivity of Simulated Solutions of
More information1 GSW Gaussian Elimination
Gaussian elimination is probabl the simplest technique for solving a set of simultaneous linear equations, such as: = A x + A x + A x +... + A x,,,, n n = A x + A x + A x +... + A x,,,, n n... m = Am,x
More information7.7 LOTKA-VOLTERRA M ODELS
77 LOTKA-VOLTERRA M ODELS sstems, such as the undamped pendulum, enabled us to draw the phase plane for this sstem and view it globall In that case we were able to understand the motions for all initial
More informationn n. ( t) ( ) = = Ay ( ) a y
Sstems of ODE Example of sstem with ODE: = a+ a = a + a In general for a sstem of n ODE = a + a + + a n = a + a + + a n = a + a + + a n n n nn n n n Differentiation of matrix: ( t) ( t) t t = = t t ( t)
More informationMAT 127: Calculus C, Fall 2010 Solutions to Midterm I
MAT 7: Calculus C, Fall 00 Solutions to Midterm I Problem (0pts) Consider the four differential equations for = (): (a) = ( + ) (b) = ( + ) (c) = e + (d) = e. Each of the four diagrams below shows a solution
More information13.1 2X2 Systems of Equations
. X Sstems of Equations In this section we want to spend some time reviewing sstems of equations. Recall there are two basic techniques we use for solving a sstem of equations: Elimination and Substitution.
More informationGreen s Theorem Jeremy Orloff
Green s Theorem Jerem Orloff Line integrals and Green s theorem. Vector Fields Vector notation. In 8.4 we will mostl use the notation (v) = (a, b) for vectors. The other common notation (v) = ai + bj runs
More information