= D ( )( ) derivative
|
|
- Philippa Franklin
- 5 years ago
- Views:
Transcription
1 > restart; In this worksheet we are going to be considering the numeric solution of the following differential equation: > ode diff(y(x,x f(x,y(x; d ode dx y( x f ( x, y( x Here, f is taken to be a known function. Everything we do will follow from the Taylor series expansion of y(x: > series y(x+h series(y(x+h,h,4; series y ( x + h y( x + D( y ( x h + ( D ( ( y ( x h + ( D ( 3 ( y( x h 3 + O( h 4 The O(h^4 term represents all the terms in the Taylor series that are not explicitly written down. If we were to neglect all the missing terms, the RHS would give an approximation to y(x+h with an error that scales like h^4. We want to used the ode to remove derivatives of y in the above Taylor series. To do this, we need to express the ode using the D notation: > derivative[] convert(ode,d; derivative D( y ( x f ( x, y( x We also need expressions for the higher order derivatives of y in terms of f. These are obtained as follows: > for i from to 3 do: derivative[i] subs(seq(derivative[j],j..i-,convert(diff(derivative[i-],x,d: od; derivative ( D ( ( y ( x + f x, y( x D ( f ( x, y( x f ( x, y( x derivative 3 ( D ( 3 ( y ( x + ( f ( x, y( x f ( x, y( x + ( + ( f ( x, y( x f ( x, y( x f ( x, y( x + D ( D + D f ( x, y( x Notice how in the above, the expressions for lower order derivatives are used to simplify the expressions for the higher ones. Putting this into the original Taylor series expansion gives > series subs(convert(derivative,list,series; series y ( x + h y( x f ( x, y( x h + h, +, f x, y( x f ( x, y( x + ( + ( f ( x, y( x f ( x, y( x f ( x, y( x + f x, y( x ( D ( f ( x, y( x + D ( f ( x, y( x f ( x, y( x h 3 + O( h 4 This series expansion is entirely in terms of y(x and f. It will be convenient to re-arrange this so that all the terms of order h^ and above are on the right. First we execute the following: > series3 expand(series - convert(series(rhs(series,h,,polynom; series3 y ( x + h y( x f ( x, y( x h y( x f ( x, y( x h + f x, y( x h, f x, y( x + 3, +, + f x, y( x D ( f ( x, y( x + f ( x, y( x h 3 + O( h 4 y( x f ( x, y( x h We can make the RHS look nicer by expanding it in another Taylor series in h. To do this, we use the applyop command (lookup?applyop. Basically, applyop(series,,equation,h,n will perform a n^th order Taylor series expansion on the RHS of equation about h 0. > series4 applyop(series,,series3,h,4; series4 y ( x + h y( x f ( x, y( x h + h, + 3, f x, y( x f ( x, y( x +, f x, y( x f ( x, y( x + f x, y( x D + f x, y( x f ( x, y( x h 3 + O( h 4 Note that this expression holds for any x and h. So we can get another series expansion by putting h to -h and x to x+h:
2 > series5 subs(h-h,xx+h,series4; series5 y( x y ( x + h + f ( x + h, y ( x + h h + f x + h, y ( x + h f x + h, y ( x + h f ( x + h, y ( x + h ( h ( f ( x + h, y ( x + h + 3, f x + h, y ( x + h f ( x + h, y ( x + h + ( f ( x + h, y ( x + h f ( x + h, y ( x + h f x + h, y ( x + h D ( f ( x + h, y ( x + h f x + h, y ( x + h f ( x + h, y ( x + h ( h 3 + O (( h 4 Unlike series4, the RHS of this is not a simple polynomial in h due to the presence of terms like f(x+h,y(x+h. These can be removed by again re-expanding the RHS as a series: > series applyop(series,,series5,h,4; series y( x y ( x + h + f ( x + h, y ( x + h h + f x, y( x f ( x, y( x h f x, y( x ( D ( f ( x, y( x + D ( f ( x, y( x D( y ( x + ( ( f ( x, y( x + ( f ( x, y( x D( y ( x f ( x, y( x + 3, +, f x, y( x D( y ( x 3, ( f ( x, y( x f ( x, y( x f x, y( x D f x, y( x f ( x, y( x h 3 + O( h 4 series4 and series form the basis of our numerical scheme. We can use them to find an approximation to the value of y(x+h given knowledge of y(x. In this sense, the give us a "new" value of y from and "old" value. We re-label the various quantities accordingly: > Subs [y(x+hy_new,y(xy_old,x+hx_new,xx_old]: series7 subs(subs,series4; series8 subs(subs,series; series7 y_new y_old f ( x_old, y_old h + y_old y_old h y_old + 3, y_old + y_old f ( x_old, y_old y_old D y_old y_old f ( x_old, y_old h 3 + O( h 4 series8 y_old y_new + f ( x_new, y_new h + y_old y_old h y_old ( D y_old + D y_old D( y ( x_old ( ( f ( x_old, y_old + ( f ( x_old, y_old D( y ( x_old f ( x_old, y_old 3 +, y_old D( y ( x_old 3, y_old y_old y_old f ( x_old, y_old y_old D y_old y_old f ( x_old, y_old h 3 + O( h 4 The "stencil" for the so-called forward Euler scheme is obtained by discarding all the terms on the right of series7. This gives us an explicit formula for y_new in terms of y_old, accurate to order h: > ForwardEuler applyop(series,,isolate(series7,y_new,h,; ForwardEuler y_new y_old + f ( x_old, y_old h + O( h The same procedure on series8 yields another formula accurate to order h. But, in this formula one cannot explicitly solve for y_new because it appears as one of the arguments of f. Hence, this is an implicit stencil where the the new value of y must be obtained by solving an equation. > BackwardEuler applyop(series,,series8,h,; BackwardEuler y_old y_new + f ( x_new, y_new h O( h Another implicit scheme is obtained by subtracting series7 from series8. Notice that the h^ term in both is the same, so this subtraction removes the h^ error term, leaving a residual error of order h^3. But again, y_new is defined implicitly in terms of y_old. > Trapezoidal applyop(series,,(series7-series8/,h,3;
3 Trapezoidal y_new y_old f ( x_old, y_old h f ( x_new, y_new h O ( h3 Actually, the backward and trapezoidal stencils are only implicit for general f. If we specialize to a particular form of f as follows, they can all be made explicit (notice we are now dropping the big-o terms: > f (x,y -> lambda*y; ForwardEuler isolate(convert(forwardeuler,polynom,y_new; BackwardEuler isolate(convert(backwardeuler,polynom,y_new; Trapezoidal isolate(convert(trapezoidal,polynom,y_new; f ( x, y λy ForwardEuler y_new y_old + λy_oldh BackwardEuler y_new y_old + λh y_old + λy_oldh Trapezoidal y_new λh It is convenient to turn these into mappings using unapply. For example, _ForwardEuler is now a procedure that take y_old, lambda and h and returns y_new using the ForwardEuler stencil. > _ForwardEuler unapply(rhs(forwardeuler,y_old,lambda,h; _BackwardEuler unapply(rhs(backwardeuler,y_old,lambda,h; _Trapezoidal unapply(rhs(trapezoidal,y_old,lambda,h; _ForwardEuler ( y_old, λ, h y_old + λy_oldh _BackwardEuler ( y_old, λ, h y_old + λh y_old + λy_oldh _Trapezoidal ( y_old, λ, h λh We now write a procedure that calculates an approximate solution for y(x using one of the above stencils. We calcuate values of y at N + discrete points from x 0 to x xf. The spacing between the points is h xf/n. We take for initial data y(0 y0. The procedure returns a list of N+ ordered pairs corresponding to our numerical solution for (x,y: > EULER proc(lambda,xf,y0,n,stencil local h, x, y, i: h xf/n: x Array(0..N,[seq(evalf(i*h,i0..N]: y Array(0..N: y[0] y0: for i from to N do: y[i] evalf(stencil(y[i-],lambda,h: od: [seq([x[i],y[i]],i0..n]: end proc: Here is the output using the trapezoidal stencil and ten points (lambda : > EULER(,,,0,_Trapezoidal; [ [ 0., ], [ ,. ], [ , ], [ , ], [ , ], [ , ], [ , ], [ , ], [ , ], [ , ], [., ] ] Here is a plot of the numeric solution obtained from the three methods versus the exact solution y exp(lambda*x [in cyan]. Experiment with different values of N and lambda and you will see that trapezoidal method [blue] generally does better than the other two [green and blue]. > N 5; lambda /; plot([euler(lambda,5,,n,_forwardeuler,euler(lambda,5,,n,_backwardeuler,euler(lambda,5,,n,_trapezoidal,exp(lambda*x],x0..5,color[red,green,blue,cyan];
4 N 5 λ The procedure EULER_ERROR calculate the discrepancy between the numerical and analytical expressions for y at x xf for the various stencils. > EULER_ERROR proc(lambda,xf,y0,n,stencil: evalf(abs(euler(lambda,xf,y0,n,stencil[n+][]-exp(lambda*xf: end proc: > EULER_ERROR(/,0,,000,_BackwardEuler; It is interesting to plot the log(error versus log(n for the schemes: > lambda -: N_low 0: N_high 400: xf 0: forward_error evalf([seq([log(i,log(euler_error(lambda,xf,,i,_forwardeuler],in_low..n_high]: backward_error evalf([seq([log(i,log(euler_error(lambda,xf,,i,_backwardeuler],in_low..n_high]: trap_error
5 evalf([seq([log(i,log(euler_error(lambda,xf,,i,_trapezoidal],in_low..n_high]: plot([forward_error,backward_error,trap_error],color[red,green,blue],axesboxed; In the log-log plot, the error in the forward an backward schemes asymptote to a line of slope - for large N, while the error curve for the trapezoidal algorithm goes to a line of slope -. This means that the error in y(xf for the first two schemes goes like N^- while it goes like N^- for the trapezoidal scheme. Since N ~ /h, this gives errors of O(h and O(h^ for the forward/backward and trapezoidal scheme, respectively. Hence, if the error in one step of a numerical scheme is O(h^(p+, then one can expect the global error in the simulation to be O(h^p. By "global error" I mean the error in the value of y(xf obtained from y(0 by taking many small steps (roughly speaking.
FINITE DIFFERENCES. Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations.
FINITE DIFFERENCES Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations. 1. Introduction When a function is known explicitly, it is easy
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI097 TU) Thursday January 014, 18:0-1:0
More informationECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations
ECE257 Numerical Methods and Scientific Computing Ordinary Differential Equations Today s s class: Stiffness Multistep Methods Stiff Equations Stiffness occurs in a problem where two or more independent
More informationExercises, module A (ODEs, numerical integration etc)
FYTN HT18 Dept. of Astronomy and Theoretical Physics Lund University, Sweden Exercises, module A (ODEs, numerical integration etc) 1. Use Euler s method to solve y (x) = y(x), y() = 1. (a) Determine y
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More informationMultistep Methods for IVPs. t 0 < t < T
Multistep Methods for IVPs We are still considering the IVP dy dt = f(t,y) t 0 < t < T y(t 0 ) = y 0 So far we have looked at Euler s method, which was a first order method and Runge Kutta (RK) methods
More informationNotation Nodes are data points at which functional values are available or at which you wish to compute functional values At the nodes fx i
LECTURE 6 NUMERICAL DIFFERENTIATION To find discrete approximations to differentiation (since computers can only deal with functional values at discrete points) Uses of numerical differentiation To represent
More informationSection 15.5: The Chain Rule. 1 Objectives. 2 Assignments. 3 Maple Commands. 1. Apply all forms of the chain rule. (3,5,9,11,13,19,21)
Section 155: The Chain Rule 1 Objectives 1 Apply all forms of the chain rule (3,5,9,11,13,19,21) 2 Compute derivatives by implicit differentiation (25,27,29) 2 Assignments 1 Read Section 155 2 Problems:
More informationTHE TRAPEZOIDAL QUADRATURE RULE
y 1 Computing area under y=1/(1+x^2) 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Trapezoidal rule x THE TRAPEZOIDAL QUADRATURE RULE From Chapter 5, we have the quadrature
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... find the approximate solutions of derivative (first- and second-order) and antiderivative (definite integral only). Numerical Differentiation
More informationThen, Lagrange s interpolation formula is. 2. What is the Lagrange s interpolation formula to find, if three sets of values. are given.
SRI RAMAKRISHNA INSTITUTE OF TECHNOLOGY, COIMBATORE- 10 DEPARTMENT OF SCIENCE AND HUMANITIES SUBJECT NUMERICAL METHODS & LINEAR PROGRAMMING ( SEMESTER IV ) IV- INTERPOLATION, NUMERICAL DIFFERENTIATION
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 informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationSolving Ordinary Differential equations
Solving Ordinary Differential equations Taylor methods can be used to build explicit methods with higher order of convergence than Euler s method. The main difficult of these methods is the computation
More informationAQR Unit 4: Using Recursion in Models and Decision Making Sequence Notes. Name: Date: Sequences
Name: Date: Sequences A number sequence is a set of numbers, usually separated by commas, arranged in an order. The first term is referred to as t 1, the second term as t 2, the third term as t 3 and so
More informationSeries Solutions. 8.1 Taylor Polynomials
8 Series Solutions 8.1 Taylor Polynomials Polynomial functions, as we have seen, are well behaved. They are continuous everywhere, and have continuous derivatives of all orders everywhere. It also turns
More informationTHEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)
4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M
More informationNonconstant Coefficients
Chapter 7 Nonconstant Coefficients We return to second-order linear ODEs, but with nonconstant coefficients. That is, we consider (7.1) y + p(t)y + q(t)y = 0, with not both p(t) and q(t) constant. The
More informationIntegration of Ordinary Differential Equations
Integration of Ordinary Differential Equations Com S 477/577 Nov 7, 00 1 Introduction The solution of differential equations is an important problem that arises in a host of areas. Many differential equations
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 information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 20
2.29 Numerical Fluid Mechanics Fall 2011 Lecture 20 REVIEW Lecture 19: Finite Volume Methods Review: Basic elements of a FV scheme and steps to step-up a FV scheme One Dimensional examples d x j x j 1/2
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 2. Predictor-Corrector Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Matematics National Institute of Tecnology Durgapur Durgapur-7109 email: anita.buie@gmail.com 1 . Capter 8 Numerical Solution of Ordinary
More informationIntroduction to standard and non-standard Numerical Methods
Introduction to standard and non-standard Numerical Methods Dr. Mountaga LAM AMS : African Mathematic School 2018 May 23, 2018 One-step methods Runge-Kutta Methods Nonstandard Finite Difference Scheme
More informationComputational Techniques Prof. Dr. Niket Kaisare Department of Chemical Engineering Indian Institute of Technology, Madras
Computational Techniques Prof. Dr. Niket Kaisare Department of Chemical Engineering Indian Institute of Technology, Madras Module No. # 07 Lecture No. # 04 Ordinary Differential Equations (Initial Value
More informationAMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.
Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error
More informationMathematical Computing
IMT2b2β Department of Mathematics University of Ruhuna A.W.L. Pubudu Thilan Differential Equations Types of Differential Equations Differential equations can basically be classified as ordinary differential
More informationSecond Order ODEs. CSCC51H- Numerical Approx, Int and ODEs p.130/177
Second Order ODEs Often physical or biological systems are best described by second or higher-order ODEs. That is, second or higher order derivatives appear in the mathematical model of the system. For
More informationSolving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations
Solving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations Outline ODEs and initial conditions. Explicit and implicit Euler methods. Runge-Kutta methods. Multistep
More informationDifferential Equations
Differential Equations Definitions Finite Differences Taylor Series based Methods: Euler Method Runge-Kutta Methods Improved Euler, Midpoint methods Runge Kutta (2nd, 4th order) methods Predictor-Corrector
More information4 Stability analysis of finite-difference methods for ODEs
MATH 337, by T. Lakoba, University of Vermont 36 4 Stability analysis of finite-difference methods for ODEs 4.1 Consistency, stability, and convergence of a numerical method; Main Theorem In this Lecture
More informationLecture V: The game-engine loop & Time Integration
Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change
More informationMath 660 Lecture 4: FDM for evolutionary equations: ODE solvers
Math 660 Lecture 4: FDM for evolutionary equations: ODE solvers Consider the ODE u (t) = f(t, u(t)), u(0) = u 0, where u could be a vector valued function. Any ODE can be reduced to a first order system,
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 informationPh 22.1 Return of the ODEs: higher-order methods
Ph 22.1 Return of the ODEs: higher-order methods -v20130111- Introduction This week we are going to build on the experience that you gathered in the Ph20, and program more advanced (and accurate!) solvers
More informationWORD: EXAMPLE(S): COUNTEREXAMPLE(S): EXAMPLE(S): COUNTEREXAMPLE(S): WORD: EXAMPLE(S): COUNTEREXAMPLE(S): EXAMPLE(S): COUNTEREXAMPLE(S): WORD:
Bivariate Data DEFINITION: In statistics, data sets using two variables. Scatter Plot DEFINITION: a bivariate graph with points plotted to show a possible relationship between the two sets of data. Positive
More informationModule Two: Differential Calculus(continued) synopsis of results and problems (student copy)
Module Two: Differential Calculus(continued) synopsis of results and problems (student copy) Srikanth K S 1 Syllabus Taylor s and Maclaurin s theorems for function of one variable(statement only)- problems.
More informationLecture 5: Single Step Methods
Lecture 5: Single Step Methods J.K. Ryan@tudelft.nl WI3097TU Delft Institute of Applied Mathematics Delft University of Technology 1 October 2012 () Single Step Methods 1 October 2012 1 / 44 Outline 1
More informationThe Definition and Numerical Method of Final Value Problem and Arbitrary Value Problem Shixiong Wang 1*, Jianhua He 1, Chen Wang 2, Xitong Li 1
The Definition and Numerical Method of Final Value Problem and Arbitrary Value Problem Shixiong Wang 1*, Jianhua He 1, Chen Wang 2, Xitong Li 1 1 School of Electronics and Information, Northwestern Polytechnical
More informationMultistage Methods I: Runge-Kutta Methods
Multistage Methods I: Runge-Kutta Methods Varun Shankar January, 0 Introduction Previously, we saw that explicit multistep methods (AB methods) have shrinking stability regions as their orders are increased.
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 informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationMA4001 Engineering Mathematics 1 Lecture 15 Mean Value Theorem Increasing and Decreasing Functions Higher Order Derivatives Implicit Differentiation
MA4001 Engineering Mathematics 1 Lecture 15 Mean Value Theorem Increasing and Decreasing Functions Higher Order Derivatives Implicit Differentiation Dr. Sarah Mitchell Autumn 2014 Rolle s Theorem Theorem
More informationDefinition 5.1. A vector field v on a manifold M is map M T M such that for all x M, v(x) T x M.
5 Vector fields Last updated: March 12, 2012. 5.1 Definition and general properties We first need to define what a vector field is. Definition 5.1. A vector field v on a manifold M is map M T M such that
More informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More informationEuler s Method, cont d
Jim Lambers MAT 461/561 Spring Semester 009-10 Lecture 3 Notes These notes correspond to Sections 5. and 5.4 in the text. Euler s Method, cont d We conclude our discussion of Euler s method with an example
More informationAM205: Assignment 3 (due 5 PM, October 20)
AM25: Assignment 3 (due 5 PM, October 2) For this assignment, first complete problems 1, 2, 3, and 4, and then complete either problem 5 (on theory) or problem 6 (on an application). If you submit answers
More information(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ).
1 Interpolation: The method of constructing new data points within the range of a finite set of known data points That is if (x i, y i ), i = 1, N are known, with y i the dependent variable and x i [x
More informationDYNAMICAL SYSTEMS
0.42 DYNAMICAL SYSTEMS Week Lecture Notes. What is a dynamical system? Probably the best way to begin this discussion is with arguably a most general and yet least helpful statement: Definition. A dynamical
More information! 1.1 Definitions and Terminology
! 1.1 Definitions and Terminology 1. Introduction: At times, mathematics aims to describe a physical phenomenon (take the population of bacteria in a petri dish for example). We want to find a function
More informationOrdinary differential equations. Phys 420/580 Lecture 8
Ordinary differential equations Phys 420/580 Lecture 8 Most physical laws are expressed as differential equations These come in three flavours: initial-value problems boundary-value problems eigenvalue
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 informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
More informationMA1023-Methods of Mathematics-15S2 Tutorial 1
Tutorial 1 the week starting from 19/09/2016. Q1. Consider the function = 1. Write down the nth degree Taylor Polynomial near > 0. 2. Show that the remainder satisfies, < if > > 0 if > > 0 3. Show that
More information10. e tan 1 (y) 11. sin 3 x
MATH B FINAL REVIEW DISCLAIMER: WHAT FOLLOWS IS A LIST OF PROBLEMS, CONCEPTUAL QUESTIONS, TOPICS, AND SAMPLE PROBLEMS FROM THE TEXTBOOK WHICH COMPRISE A HEFTY BUT BY NO MEANS EXHAUSTIVE LIST OF MATERIAL
More informationAN OVERVIEW. Numerical Methods for ODE Initial Value Problems. 1. One-step methods (Taylor series, Runge-Kutta)
AN OVERVIEW Numerical Methods for ODE Initial Value Problems 1. One-step methods (Taylor series, Runge-Kutta) 2. Multistep methods (Predictor-Corrector, Adams methods) Both of these types of methods are
More informationGoals for This Lecture:
Goals for This Lecture: Learn the Newton-Raphson method for finding real roots of real functions Learn the Bisection method for finding real roots of a real function Look at efficient implementations of
More informationMathematics for chemical engineers. Numerical solution of ordinary differential equations
Mathematics for chemical engineers Drahoslava Janovská Numerical solution of ordinary differential equations Initial value problem Winter Semester 2015-2016 Outline 1 Introduction 2 One step methods Euler
More informationSolutions Definition 2: a solution
Solutions As was stated before, one of the goals in this course is to solve, or find solutions of differential equations. In the next definition we consider the concept of a solution of an ordinary differential
More information1 Implicit Differentiation
MATH 1010E University Mathematics Lecture Notes (week 5) Martin Li 1 Implicit Differentiation Sometimes a function is defined implicitly by an equation of the form f(x, y) = 0, which we think of as a relationship
More informationNUMERICAL SOLUTION OF ODE IVPs. Overview
NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental
More 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 informationLecture 6. Numerical Solution of Differential Equations B21/B1. Lecture 6. Lecture 6. Initial value problem. Initial value problem
Lecture 6 Numerical Solution of Differential Equations B2/B Initial value problem y = f(x,y), x [x,x M ] y(x ) = y Professor Endre Süli Numerical Solution of Differential Equations p./45 Numerical Solution
More information(1) Rate of change: A swimming pool is emptying at a constant rate of 90 gal/min.
CHAPTER 1 Introduction 1. Bacground Models of physical situations from Calculus (1) Rate of change: A swimming pool is emptying at a constant rate of 90 gal/min. With V = volume in gallons and t = time
More informationBindel, Spring 2012 Intro to Scientific Computing (CS 3220) Week 12: Monday, Apr 16. f(x) dx,
Panel integration Week 12: Monday, Apr 16 Suppose we want to compute the integral b a f(x) dx In estimating a derivative, it makes sense to use a locally accurate approximation to the function around the
More informationNumerical integration - I. M. Peressi - UniTS - Laurea Magistrale in Physics Laboratory of Computational Physics - Unit V
Numerical integration - I M. Peressi - UniTS - Laurea Magistrale in Physics Laboratory of Computational Physics - Unit V deterministic methods in 1D equispaced points (trapezoidal, Simpson...), others...
More informationEngg. Math. II (Unit-IV) Numerical Analysis
Dr. Satish Shukla of 33 Engg. Math. II (Unit-IV) Numerical Analysis Syllabus. Interpolation and Curve Fitting: Introduction to Interpolation; Calculus of Finite Differences; Finite Difference and Divided
More informationChapter 11 ORDINARY DIFFERENTIAL EQUATIONS
Chapter 11 ORDINARY DIFFERENTIAL EQUATIONS The general form of a first order differential equations is = f(x, y) with initial condition y(a) = y a We seek the solution y = y(x) for x > a This is shown
More informationDIFFERENTIATION RULES
3 DIFFERENTIATION RULES DIFFERENTIATION RULES The functions that we have met so far can be described by expressing one variable explicitly in terms of another variable. y For example,, or y = x sin x,
More informationI. Numerical Computing
I. Numerical Computing A. Lectures 1-3: Foundations of Numerical Computing Lecture 1 Intro to numerical computing Understand difference and pros/cons of analytical versus numerical solutions Lecture 2
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 5 Chapter 21 Numerical Differentiation PowerPoints organized by Dr. Michael R. Gustafson II, Duke University 1 All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationOn the General Solution of Initial Value Problems of Ordinary Differential Equations Using the Method of Iterated Integrals. 1
On the General Solution of Initial Value Problems of Ordinary Differential Equations Using the Method of Iterated Integrals. 1 Ahsan Amin ahsanamin2999@gmail.com This First version January 2017 First Preliminary
More informationMAPLE Worksheet Number 7 Derivatives in Calculus
MAPLE Worksheet Number 7 Derivatives in Calculus The MAPLE command for computing the derivative of a function depends on which of the two ways we used to define the function, as a symbol or as an operation.
More informationComputational Physics
Interpolation, Extrapolation & Polynomial Approximation Lectures based on course notes by Pablo Laguna and Kostas Kokkotas revamped by Deirdre Shoemaker Spring 2014 Introduction In many cases, a function
More informationDepartment of Applied Mathematics Preliminary Examination in Numerical Analysis August, 2013
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August, 013 August 8, 013 Solutions: 1 Root Finding (a) Let the root be x = α We subtract α from both sides of x n+1 = x
More informationTopics Covered in Calculus BC
Topics Covered in Calculus BC Calculus BC Correlation 5 A Functions, Graphs, and Limits 1. Analysis of graphs 2. Limits or functions (including one sides limits) a. An intuitive understanding of the limiting
More informationCS205B / CME306 Homework 3. R n+1 = R n + tω R n. (b) Show that the updated rotation matrix computed from this update is not orthogonal.
CS205B / CME306 Homework 3 Rotation Matrices 1. The ODE that describes rigid body evolution is given by R = ω R. (a) Write down the forward Euler update for this equation. R n+1 = R n + tω R n (b) Show
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS ( WI3097 TU AESB0 ) Thursday April 6
More informationMA 3021: Numerical Analysis I Numerical Differentiation and Integration
MA 3021: Numerical Analysis I Numerical Differentiation and Integration Suh-Yuh Yang ( 楊肅煜 ) Department of Mathematics, National Central University Jhongli District, Taoyuan City 32001, Taiwan syyang@math.ncu.edu.tw
More informationEuler s Method and Logistic Growth (BC Only)
Euler s Method Students should be able to: Approximate numerical solutions of differential equations using Euler s method without a calculator. Recognize the method as a recursion formula extension of
More information23. Implicit differentiation
23. 23.1. The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x and y. If a value of x is given, then a corresponding value of y is determined. For instance, if x = 1, then y
More informationReview: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function.
Taylor Series (Sect. 10.8) Review: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function. Review: Power series define functions Remarks:
More informationNUMERICAL METHODS. lor CHEMICAL ENGINEERS. Using Excel', VBA, and MATLAB* VICTOR J. LAW. CRC Press. Taylor & Francis Group
NUMERICAL METHODS lor CHEMICAL ENGINEERS Using Excel', VBA, and MATLAB* VICTOR J. LAW CRC Press Taylor & Francis Group Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup,
More informationMB4018 Differential equations
MB4018 Differential equations Part II http://www.staff.ul.ie/natalia/mb4018.html Prof. Natalia Kopteva Spring 2015 MB4018 (Spring 2015) Differential equations Part II 0 / 69 Section 1 Second-Order Linear
More informationTHE UNIVERSITY OF NEW SOUTH WALES SCHOOL OF PHYSICS PHYS2020 COMPUTATIONAL PHYSICS FINAL EXAM SESSION I Answer all questions
THE UNIVERSITY OF NEW SOUTH WALES SCHOOL OF PHYSICS PHYS2020 COMPUTATIONAL PHYSICS FINAL EXAM SESSION I 2007 Answer all questions Time allowed =2 hours Total number of questions =5 Marks =40 The questions
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Multistep Methods Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 1. Runge Kutta methods 2. Embedded RK methods
More informationBackward error analysis
Backward error analysis Brynjulf Owren July 28, 2015 Introduction. The main source for these notes is the monograph by Hairer, Lubich and Wanner [2]. Consider ODEs in R d of the form ẏ = f(y), y(0) = y
More informationReview Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester
HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE
More informationNovember 20, Interpolation, Extrapolation & Polynomial Approximation
Interpolation, Extrapolation & Polynomial Approximation November 20, 2016 Introduction In many cases we know the values of a function f (x) at a set of points x 1, x 2,..., x N, but we don t have the analytic
More informationHUDSONVILLE HIGH SCHOOL COURSE FRAMEWORK
HUDSONVILLE HIGH SCHOOL COURSE FRAMEWORK COURSE / SUBJECT A P C a l c u l u s ( B C ) KEY COURSE OBJECTIVES/ENDURING UNDERSTANDINGS OVERARCHING/ESSENTIAL SKILLS OR QUESTIONS Limits and Continuity Derivatives
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 information1.1 Motivation: Why study differential equations?
Chapter 1 Introduction Contents 1.1 Motivation: Why stu differential equations?....... 1 1.2 Basics............................... 2 1.3 Growth and decay........................ 3 1.4 Introduction to Ordinary
More informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon's method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
More informationMini project ODE, TANA22
Mini project ODE, TANA22 Filip Berglund (filbe882) Linh Nguyen (linng299) Amanda Åkesson (amaak531) October 2018 1 1 Introduction Carl David Tohmé Runge (1856 1927) was a German mathematician and a prominent
More informationPower Series Solutions We use power series to solve second order differential equations
Objectives Power Series Solutions We use power series to solve second order differential equations We use power series expansions to find solutions to second order, linear, variable coefficient equations
More informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
More informationNumerical Methods for Initial Value Problems; Harmonic Oscillators
1 Numerical Methods for Initial Value Problems; Harmonic Oscillators Lab Objective: Implement several basic numerical methods for initial value problems (IVPs), and use them to study harmonic oscillators.
More informationThis practice exam is intended to help you prepare for the final exam for MTH 142 Calculus II.
MTH 142 Practice Exam Chapters 9-11 Calculus II With Analytic Geometry Fall 2011 - University of Rhode Island This practice exam is intended to help you prepare for the final exam for MTH 142 Calculus
More information