Euler s Methods (a family of Runge- Ku9a methods)
|
|
- Winifred Shelton
- 6 years ago
- Views:
Transcription
1 Euler s Methods (a family of Runge- Ku9a methods)
2 ODE IVP An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable: The equation is coupled with an initial value/condition (i.e., value of y at x = 0). Hence, such equation is usually called ODE Initial Value Problem (IVP). Different notations for derivatives: Space/position: Time:
3 Recall Taylor Series (Lecture 1) If the y-axis and origin are moved a units to the left, the equation of the same curve relative to the new axis becomes y = f(a+x) and the function value at P is f(a). y Q y=f(a+x) P f(0) f(a) f(a+h) At point Q: 0 a h x
4 Euler s Method from Taylor Series The approximation used with Euler s method is to take only the first two terms of the Taylor series: In general form: new value = old value + step size slope If and as well as then (8.1)
5 Euler s (Forward) Method Alternatively, from step size and rearrange to we use the Taylor series to approximate the function around with step size Taking only the first derivative: where and is the differential equation evaluated at (8.2) This formula is referred to as Euler s forward method, or explicit Euler s method, or Euler-Cauchy method, or pointslope method.
6 Example 8.1 Obtain a numerical solution of the differential equation given the initial conditions that when the range with intervals of 0.2. for
7 Example 8.1 Obtain a numerical solution of the differential equation given the initial conditions that when the range with intervals of 0.2. for Solution: with By Euler s method:
8 Example 8.1 At where then If the step by step Euler s method is continued, we can present the results in a table: x i y i y i
9 Example 8.1 To compare the approximations from Euler s method with the exact solution, the ODE can be solved analytically using integrating factor method. Rearrange the equation into the form (8.3) so now where and Take the integrating factor and with The general solution of the equation (8.3) is: thus (8.4)
10 Example 8.1 from which The second term on the right hand side is solved using integration by parts. Hence When thus from which Therefore or:
11 Example 8.1
12 MATLAB ImplementaMons function [x,y] = EulerForward(f,xinit,xend,yinit,h) % Number of iterations N = (xend-xinit)/h; % Initialize arrays % The first elements take xinit and yinit, correspondingly, % the rest fill with 0s. x = zeros(1, N+1); y = zeros(1, N+1); x(1) = xinit; y(1) = yinit; for i=1:n x(i+1) = x(i)+h; y(i+1) = y(i) + h*feval(f,x(i),y(i)); end end
13 MATLAB ImplementaMons function dydx = EulerFunction(x,y) dydx = 3*(1+x)-y; end a = 1; b = 2; ya = 4; h = 0.2; N = (b-a)/h; t = a:h:b; [x,y] = EulerForward('EulerFunction', a, b, ya, h); ye = y; % Numerical solution from using Euler s forward method yi = 3*t+exp(1-t); % Exact solution hold on; plot(t,yi,'r','linewidth', 2); plot(t,ye,'b','linewidth', 2); hold off; box on;
14 Error Analysis The numerical solution of ODEs involves two types of error: (1) Truncation, or discretization errors caused by the nature of the techniques used to approximate values of y. (2) Round-off errors, caused by the limited numbers of significant digits that can be retained by a computer. The truncation errors are composed of two parts: (a) Local truncation error à from an application of the method in question over a single step. (b) Propagated truncation error à from approximations produced during the previous step The sum of (a) and (b) is the global truncation error.
15 Error Analysis The local truncation error is the difference between the numerical solution after one step y 1 and the exact solution at Using equation (8.2), the numerical solution is given by: (8.5) where For the exact solution, we use the Taylor expansion of the function y around x 0 : (8.6)
16 Error Analysis Subtracting equation (8.6) from (8.5) yields The difference between the exact and the numerical solutions is the true local truncation error E t : The general form of the true local truncation error is given by: E t = f 0 (x i,y i ) 2! h 2 + f 00 (x i,y i ) 3! (8.7) h O(h n+1 ) (8.8)
17 Error Analysis The expressions (8.7) and (8.8) show that for sufficiently small h, the errors usually decrease as the order increases, the result is represented as the approximate local truncation error, formulated as (8.9) or (8.10)
18 Error Analysis Because Euler s method uses straight-line segments to approximate the solution, hence the method is also referred to as a first-order method. From error estimates shown in Lecture 7, then the global truncation error In general, if the local truncation error is the global truncation error is
19 Example 8.2 Use eqn. (8.8) to estimate the error of the first step of the following equation when integrating from x = 0 to x = 4 dy dx = 2x3 + 12x 2 20x +8.5 with step size h = 0.5. Use the results to determine the error due to each higher-order term of the Taylor series expansion. Solution: For the above problem, eqn. (8.8) can be written in: E t = f 0 (x i,y i ) 2! h 2 + f 00 (x i,y i ) 3! h 3 + f 000 (x i,y i ) 4! h 4
20 Example 8.2 where f 0 (x i,y i )= 6x x 20 f 00 (x i,y i )= 12x + 24 f 000 (x i,y i )= 12 with derivative of order higher than 3 is equal to zero. The error due to truncation of the second-order term can be calculated as ( 6)(0) + (24)(0) 20 E t,2 = (0.5) 2 = 2.5 2
21 Example 8.2 The error due to truncation of the third-order term can be calculated as ( 12)(0) + 24 E t,3 = (0.5) 3 =0.5 6 and the error due to truncation of the fourth-order term can be calculated as E t,4 = (0.5)4 = The three results can be added to yield the total truncation error: E t = E t,2 + E t,3 + E t,4 =
22 Example 8.2 If the problem in question is equipped with additional information, i.e., initial condition where at x = 0, y = 1, the global error can be calculated using formula E t =true approximate We use eqn. (8.2): y i+1 = y i + f(x i,y i )h where y(0) = 1 and the slope at x = 0 is f(0, 1) = 2(0) (0) 2 20(0) =8.5 therefore y(0.5) = y(0) + f(0, 1)0.5 =5.25
23 Example 8.2 The true solution at x = 0.5 is y = 0.5x 4 +4x 3 10x x +1 = thus the global error is E t = =
24 Effect of Reduced Step on Euler s Method From the previous example Reducing the interval or step size to 0.1 within the range we get: x i y i y i
25 Effect of Reduced Step on Euler s Method
26 Improvements of Euler s Method (1) Heun s Method Recall in Euler s method: is used to extrapolate linearly to (8.11) (8.12) which is called a predictor equation. It provides an estimate of that allows the calculation of an estimated slope at the end of the interval: (8.13)
27 Improvements of Euler s Method Thus, the two slopes in (8.11) and (8.13) can be combined to obtain an average slope The average slope is used to extrapolate linearly from y i to y i+1 using Euler s method: Or Predictor: Corrector: (8.14) (8.15)
28 Improvements of Euler s Method (2) Midpoint Method To predict a value of y at the midpoint of the interval: This predicted value is used to calculate a slope at the midpoint:
29 Improvements of Euler s Method The slope is then used to extrapolate linearly from to get to (8.16)
30 Stability The (explicit) forward Euler s method is easy to implement. The drawback arises from the limitations on the step size to ensure numerical stability. To see this, let s examine a linear IVP given by (8.17) with The exact solution of (8.17) is and a very smooth solution with, which is a stable
31 Stability Applying the (explicit) forward Euler s method: The solution is decaying (stable) if To prevent the amplification of the errors in the iteration process, we require, for stability of the forward Euler s method: (8.18)
32 Backward Euler s Method The backward method computes the approximations using (8.19) which is an implicit method, in the sense that in order to find y i+1 the nonlinear equation (8.19) has to be solved. Using equation (8.17) gives us and the solution is decaying (stable) if. (8.20)
33 Example 8.3 Using the same example as in the forward method with h = 0.2 Solution: with By Euler s backward method: (8.21) where Substitute this into (8.21) and
34 Example 8.3 now with we solve for And we use it to get the slope
35 Example 8.3 The next iteration is where and the next point of iteration we solve for And continue the iteration until
36 Example 8.3
37 MATLAB ImplementaMon Write Matlab solver to approximate the following function using Heun s Method 1 x dy dx +4y =2 given the initial conditions x = 0 when y = 4 within the range x = 0 to x = 2 with intervals of 0.1.
38 MATLAB ImplementaMon function [x,y] = HeunsMethod(f,xinit,xend,yinit,h) % Number of iterations N = (xend-xinit)/h; % Initialize arrays % The first elements take xinit and yinit, correspondingly, % the rest fill with 0s. x = [xinit zeros(1, N)]; y = [yinit zeros(1, N)]; for i=1:n x(i+1) = x(i)+h; % Predictor ynew = y(i) + h*feval(f,x(i),y(i)); % Corrector y(i+1) = y(i) + h*feval(f,x(i),y(i))/2 + h*feval(f,x(i+1),ynew)/2; end end
39 MATLAB ImplementaMon function dydx = HeunFunction(x,y) dydx = 2*x 4*x*y; end clear all close all a = 0; b = 2; ya = 4; h = 0.1; t = a:h:b; [x,y] = HeunsMethod( HeunFunction', a, b, ya, h); ye = y; % Numerical solution from using Euler s backward method yi = 0.5*(1+7*exp(-2*x.^2)); % Exact solution, for comparison err = abs(yi ye)*100/yi; % relative error hold on; plot(t,yi,'r','linewidth', 2); plot(t,ye,'b','linewidth', 2); hold off; box on;
ODE IVP. An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable:
Euler s Methods ODE IVP An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable: The equation is coupled with an initial value/condition (i.e., value
More informationover The idea is to construct an algorithm to solve the IVP ODE (8.1)
Runge- Ku(a Methods Review of Heun s Method (Deriva:on from Integra:on) The idea is to construct an algorithm to solve the IVP ODE (8.1) over To obtain the solution point we can use the fundamental theorem
More informationover The idea is to construct an algorithm to solve the IVP ODE (9.1)
Runge- Ku(a Methods Review of Heun s Method (Deriva:on from Integra:on) The idea is to construct an algorithm to solve the IVP ODE (9.1) over To obtain the solution point we can use the fundamental theorem
More informationOrdinary Differential Equations
Next: Partial Differential Equations Up: Numerical Analysis for Chemical Previous: Numerical Differentiation and Integration Subsections Runge-Kutta Methods Euler's Method Improvement of Euler's Method
More informationModule 2: Single Step Methods Lecture 4: The Euler Method. The Lecture Contains: The Euler Method. Euler's Method (Analytical Interpretations)
The Lecture Contains: The Euler Method Euler's Method (Analytical Interpretations) An Analytical Example file:///g /Numerical_solutions/lecture4/4_1.htm[8/26/2011 11:14:40 AM] We shall now describe methods
More informationLecture 8: Euler s Methods
Lecture 8: Euler s Methods Forward Euler Method The formula for the forward Euler method is given by equation (8.2) in the lecture note for week 8, as y i+1 = y i + f(x i, y i )h. (1) where f(x i, y i
More informationReview Initial Value Problems Euler s Method Summary
THE EULER METHOD P.V. Johnson School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 INITIAL VALUE PROBLEMS The Problem Posing a Problem 3 EULER S METHOD Method Errors 4 SUMMARY OUTLINE 1 REVIEW 2 INITIAL
More informationMath 225 Scientific Computing II Outline of Lectures
Math 225 Scientific Computing II Outline of Lectures Spring Semester 2003 I. Interpolating polynomials Lagrange formulation of interpolating polynomial Uniqueness of interpolating polynomial of degree
More informationODEs occur quite often in physics and astrophysics: Wave Equation in 1-D stellar structure equations hydrostatic equation in atmospheres orbits
Solving ODEs General Stuff ODEs occur quite often in physics and astrophysics: Wave Equation in 1-D stellar structure equations hydrostatic equation in atmospheres orbits need workhorse solvers to deal
More informationQueens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 2018 16 Lecture 16 May 3, 2018 Numerical solution of systems
More informationx n x n stepnumber k order r error constant C r+1 1/2 5/12 3/8 251/720 abs. stab. interval (α,0) /11-3/10
MATH 573 LECTURE NOTES 77 13.8. Predictor-corrector methods. We consider the Adams methods, obtained from the formula xn+1 xn+1 y(x n+1 y(x n ) = y (x)dx = f(x,y(x))dx x n x n by replacing f by an interpolating
More informationChE 400: Applied Chemical Engineering Calculations Tutorial 6: Numerical Solution of ODE Using Excel and Matlab
ChE 400: Applied Chemical Engineering Calculations Tutorial 6: Numerical Solution of ODE Using Excel and Matlab Tutorial 6: Numerical Solution of ODE Gerardine G. Botte This handout contains information
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationSimulation in Computer Graphics. Particles. Matthias Teschner. Computer Science Department University of Freiburg
Simulation in Computer Graphics Particles Matthias Teschner Computer Science Department University of Freiburg Outline introduction particle motion finite differences system of first order ODEs second
More informationHomework Set #2-3, Math 475B
Homework Set #2-3, Math 475B Part I: Matlab In the last semester you learned a number of essential features of MATLAB. 1. In this instance, you will learn to make 3D plots and contour plots of z = f(x,
More informationFinite difference methods
Finite difference methods Siltanen/Railo/Kaarnioja Spring 8 Applications of matrix computations Applications of matrix computations Finite difference methods Spring 8 / Introduction Finite difference methods
More informationMathematics for chemical engineers
Mathematics for chemical engineers Drahoslava Janovská Department of mathematics Winter semester 2013-2014 Numerical solution of ordinary differential equations Initial value problem Outline 1 Introduction
More informationHomework 14 solutions
Section 9.1: Ex 1,5,10,15,16 Section 9.: Ex 1,8,9; AP 1-5,9 Section 9.3: AP 1-5 Section 9.1 Homework 14 solutions 1. (a) Show that y( is a solution to the differential equation by substitution. (b) Find
More informationDirection Fields; Euler s Method
Direction Fields; Euler s Method It frequently happens that we cannot solve first order systems dy (, ) dx = f xy or corresponding initial value problems in terms of formulas. Remarkably, however, this
More informationMathematical Methods and Modeling Laboratory class. Numerical Integration of Ordinary Differential Equations
Mathematical Methods and Modeling Laboratory class Numerical Integration of Ordinary Differential Equations Exact Solutions of ODEs Cauchy s Initial Value Problem in normal form: Recall: if f is locally
More informationComputer Simulations
Computer Simulations A practical approach to simulation Semra Gündüç gunduc@ankara.edu.tr Ankara University Faculty of Engineering, Department of Computer Engineering 2014-2015 Spring Term Ankara University
More informationSection 18-1: Graphical Representation of Linear Equations and Functions
Section 18-1: Graphical Representation of Linear Equations and Functions Prepare a table of solutions and locate the solutions on a coordinate system: f(x) = 2x 5 Learning Outcome 2 Write x + 3 = 5 as
More informationMATH2071: LAB 2: Explicit ODE methods
MATH2071: LAB 2: Explicit ODE methods 1 Introduction Introduction Exercise 1 Euler s method review Exercise 2 The Euler Halfstep (RK2) Method Exercise 3 Runge-Kutta Methods Exercise 4 The Midpoint Method
More informationThe listing says y(1) = How did the computer know this?
18.03 Class 2, Feb 3, 2010 Numerical Methods [1] How do you know the value of e? [2] Euler's method [3] Sources of error [4] Higher order methods [1] The study of differential equations has three parts:.
More informationdy = dx y dx Project 1 Consider the following differential equations: I. xy II. III.
Project 1 Consider the following differential equations: dy = dy = y dy x + sin( x) y = x I. xy II. III. 1. Graph the direction field in the neighborhood of the origin. Graph approximate integral curves
More informationOrdinary differential equations solving methods
Radim Hošek, A07237 radhost@students.zcu.cz Ordinary differential equations solving methods Problem: y = y2 (1) y = x y (2) y = sin ( + y 2 ) (3) Where it is possible we try to solve the equations analytically,
More informationAMSC/CMSC 460 Final Exam, Fall 2007
AMSC/CMSC 460 Final Exam, Fall 2007 Show all work. You may leave arithmetic expressions in any form that a calculator could evaluate. By putting your name on this paper, you agree to abide by the university
More informationModule1: Numerical Solution of Ordinary Differential Equations. Lecture 6. Higher order Runge Kutta Methods
Module1: Numerical Solution of Ordinary Differential Equations Lecture 6 Higher order Runge Kutta Methods Keywords: higher order methods, functional evaluations, accuracy Higher order Runge Kutta Methods
More informationTruncation Errors. Applied Numerical Methods with MATLAB for Engineers and Scientists, 2nd ed., Steven C. Chapra, McGraw Hill, 2008, Ch. 4.
Chapter 4: Roundoff and Truncation Errors Applied Numerical Methods with MATLAB for Engineers and Scientists, 2nd ed., Steven C. Chapra, McGraw Hill, 2008, Ch. 4. 1 Outline Errors Accuracy and Precision
More information2.1 Derivatives and Rates of Change
2.1 Derivatives and Rates of Change In this chapter we study a special type of limit, called a derivative, that occurs when we want to find a slope of a tangent line, or a velocity, or any instantaneous
More informationApplication 2.4 Implementing Euler's Method
Application 2.4 Implementing Euler's Method One's understanding of a numerical algorithm is sharpened by considering its implementation in the form of a calculator or computer program. Figure 2.4.13 in
More informationDifferentiation. The Derivative and the Tangent Line Problem 10/9/2014. Copyright Cengage Learning. All rights reserved.
Differentiation Copyright Cengage Learning. All rights reserved. The Derivative and the Tangent Line Problem Copyright Cengage Learning. All rights reserved. 1 Objectives Find the slope of the tangent
More informationFirst of all, we need to know what it means for a parameterize curve to be differentiable. FACT:
CALCULUS WITH PARAMETERIZED CURVES In calculus I we learned how to differentiate and integrate functions. In the chapter covering the applications of the integral, we learned how to find the length of
More informationA Study on Numerical Exact Solution of Euler, Improved Euler and Runge - Kutta Method
A Study on Numerical Exact Solution of Euler, Improved Euler and Runge - Kutta Method Mrs. P. Sampoornam P.K.R Arts College for Women, Gobichettipalayam Abstract: In this paper, I will discuss the Runge
More informationReals 1. Floating-point numbers and their properties. Pitfalls of numeric computation. Horner's method. Bisection. Newton's method.
Reals 1 13 Reals Floating-point numbers and their properties. Pitfalls of numeric computation. Horner's method. Bisection. Newton's method. 13.1 Floating-point numbers Real numbers, those declared to be
More informationAccuracy versus precision
Accuracy versus precision Accuracy is a consistent error from the true value, but not necessarily a good or precise error Precision is a consistent result within a small error, but not necessarily anywhere
More information8 Piecewise Polynomial Interpolation
Applied Math Notes by R. J. LeVeque 8 Piecewise Polynomial Interpolation 8. Pitfalls of high order interpolation Suppose we know the value of a function at several points on an interval and we wish to
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 6:
file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture6/6_1.htm 1 of 1 6/20/2012 12:24 PM The Lecture deals with: ADI Method file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture6/6_2.htm 1 of 2 6/20/2012
More informationCALCULUS MADE EASY - FUNCTIONALITY for the TiNspire CAS
CALCULUS MADE EASY - FUNCTIONALITY for the TiNspire CAS www.tinspireapps.com Functions READ: Linear Functions Find Slope Find y=mx+b All-in-one-Function Explorer Evaluate Function Find Domain of f(x) Find
More informationRasterization: Geometric Primitives
Rasterization: Geometric Primitives Outline Rasterizing lines Rasterizing polygons 1 Rasterization: What is it? How to go from real numbers of geometric primitives vertices to integer coordinates of pixels
More informationThe Euler and trapezoidal stencils to solve d d x y x = f x, y x
restart; Te Euler and trapezoidal stencils to solve d d x y x = y x Te purpose of tis workseet is to derive te tree simplest numerical stencils to solve te first order d equation y x d x = y x, and study
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 5 Chapter 19 Numerical Differentiation PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More information4 Visualization and. Approximation
4 Visualization and Approximation b A slope field for the differential equation y tan(x + y) tan(x) tan(y). It is not always possible to write down an explicit formula for the solution to a differential
More information(Part - 1) P. Sam Johnson. April 14, Numerical Solution of. Ordinary Differential Equations. (Part - 1) Sam Johnson. NIT Karnataka.
P. April 14, 2015 1/51 Overview We discuss the following important methods of solving ordinary differential equations of first / second order. Picard s method of successive approximations (Method of successive
More information7.1 First-order initial-value problems
7.1 First-order initial-value problems A first-order initial-value problem (IVP) is a first-order ordinary differential equation with a specified value: d dt y t y a y f t, y t That is, this says that
More informationMAT 275 Laboratory 3 Numerical Solutions by Euler and Improved Euler Methods (scalar equations)
MATLAB sessions: Laboratory 3 1 MAT 275 Laboratory 3 Numerical Solutions by Euler and Improved Euler Methods (scalar equations) In this session we look at basic numerical methods to help us understand
More informationDOWNLOAD PDF BIG IDEAS MATH VERTICAL SHRINK OF A PARABOLA
Chapter 1 : BioMath: Transformation of Graphs Use the results in part (a) to identify the vertex of the parabola. c. Find a vertical line on your graph paper so that when you fold the paper, the left portion
More informationQuasilinear First-Order PDEs
MODULE 2: FIRST-ORDER PARTIAL DIFFERENTIAL EQUATIONS 16 Lecture 3 Quasilinear First-Order PDEs A first order quasilinear PDE is of the form a(x, y, z) + b(x, y, z) x y = c(x, y, z). (1) Such equations
More informationEuler s Method for Approximating Solution Curves
Euler s Method for Approximating Solution Curves As you may have begun to suspect at this point, time constraints will allow us to learn only a few of the many known methods for solving differential equations.
More informationMath Homework 3
Math 0 - Homework 3 Due: Friday Feb. in class. Write on your paper the lab section you have registered for.. Staple the sheets together.. Solve exercise 8. of the textbook : Consider the following data:
More informationOutput Primitives Lecture: 3. Lecture 3. Output Primitives. Assuming we have a raster display, a picture is completely specified by:
Lecture 3 Output Primitives Assuming we have a raster display, a picture is completely specified by: - A set of intensities for the pixel positions in the display. - A set of complex objects, such as trees
More informationCALCULUS WITH PHYSICS APPLICATIONS - FUNCTIONALITY for the TiNspire CAS CX
CALCULUS WITH PHYSICS APPLICATIONS - FUNCTIONALITY for the TiNspire CAS CX www.tinspireapps.com Physics Apps Center of Mass (2D) 1. Moment of Mass about x and y-axis Mass of Lamina - f(x) Mass of Lamina
More information37 Self-Assessment. 38 Two-stage Runge-Kutta Methods. Chapter 10 Runge Kutta Methods
Chapter 10 Runge Kutta Methods In the previous lectures, we have concentrated on multi-step methods. However, another powerful set of methods are known as multi-stage methods. Perhaps the best known of
More informationDigital Differential Analyzer Bresenhams Line Drawing Algorithm
Bresenham s Line Generation The Bresenham algorithm is another incremental scan conversion algorithm. The big advantage of this algorithm is that, it uses only integer calculations. Difference Between
More informationOpenGL Graphics System. 2D Graphics Primitives. Drawing 2D Graphics Primitives. 2D Graphics Primitives. Mathematical 2D Primitives.
D Graphics Primitives Eye sees Displays - CRT/LCD Frame buffer - Addressable pixel array (D) Graphics processor s main function is to map application model (D) by projection on to D primitives: points,
More informationpractice: quadratic functions [102 marks]
practice: quadratic functions [102 marks] A quadratic function, f(x) = a x 2 + bx, is represented by the mapping diagram below. 1a. Use the mapping diagram to write down two equations in terms of a and
More informationAsymptotic Error Analysis
Asymptotic Error Analysis Brian Wetton Mathematics Department, UBC www.math.ubc.ca/ wetton PIMS YRC, June 3, 2014 Outline Overview Some History Romberg Integration Cubic Splines - Periodic Case More History:
More informationMAT 275 Laboratory 3 Numerical Solutions by Euler and Improved Euler Methods (scalar equations)
MAT 275 Laboratory 3 Numerical Solutions by Euler and Improved Euler Methods (scalar equations) In this session we look at basic numerical methods to help us understand the fundamentals of numerical approximations.
More informationLecture 22: Rules for Solving IVP
cs412: introduction to numerical analysis 12/7/10 Lecture 22: Rules for Solving IVP Instructor: Professor Amos Ron Scribes: Yunpeng Li, Pratap Ramamurthy, Mark Cowlishaw, Nathanael Fillmore 1 Review Recall
More information13.1. Functions of Several Variables. Introduction to Functions of Several Variables. Functions of Several Variables. Objectives. Example 1 Solution
13 Functions of Several Variables 13.1 Introduction to Functions of Several Variables Copyright Cengage Learning. All rights reserved. Copyright Cengage Learning. All rights reserved. Objectives Understand
More information2.9 Linear Approximations and Differentials
2.9 Linear Approximations and Differentials 2.9.1 Linear Approximation Consider the following graph, Recall that this is the tangent line at x = a. We had the following definition, f (a) = lim x a f(x)
More informationApproximate Solutions to Differential Equations Slope Fields (graphical) and Euler s Method (numeric) by Dave Slomer
Approximate Solutions to Differential Equations Slope Fields (graphical) and Euler s Method (numeric) by Dave Slomer Leonhard Euler was a great Swiss mathematician. The base of the natural logarithm function,
More informationRendering. A simple X program to illustrate rendering
Rendering A simple X program to illustrate rendering The programs in this directory provide a simple x based application for us to develop some graphics routines. Please notice the following: All points
More information16.1 Runge-Kutta Method
70 Chapter 6. Integration of Ordinary Differential Equations CITED REFERENCES AND FURTHER READING: Gear, C.W. 97, Numerical Initial Value Problems in Ordinary Differential Equations (Englewood Cliffs,
More informationCOMPUTER VISION > OPTICAL FLOW UTRECHT UNIVERSITY RONALD POPPE
COMPUTER VISION 2017-2018 > OPTICAL FLOW UTRECHT UNIVERSITY RONALD POPPE OUTLINE Optical flow Lucas-Kanade Horn-Schunck Applications of optical flow Optical flow tracking Histograms of oriented flow Assignment
More informationNumerical Methods Lecture 7 - Optimization
Numerical Methods Lecture 7 - Optimization Topics: numerical optimization - Newton again - Random search - Golden Section Search READING : text pgs. 331-349 Optimization - motivation What? Locating where
More informationChemical Reaction Engineering - Part 4 - Integration Richard K. Herz,
Chemical Reaction Engineering - Part 4 - Integration Richard K. Herz, rherz@ucsd.edu, www.reactorlab.net Integrating the component balance for a batch reactor For a single reaction in a batch reactor,
More informationMotion. 1 Introduction. 2 Optical Flow. Sohaib A Khan. 2.1 Brightness Constancy Equation
Motion Sohaib A Khan 1 Introduction So far, we have dealing with single images of a static scene taken by a fixed camera. Here we will deal with sequence of images taken at different time intervals. Motion
More informationPhysics Tutorial 2: Numerical Integration Methods
Physics Tutorial 2: Numerical Integration Methods Summary The concept of numerically solving differential equations is explained, and applied to real time gaming simulation. Some objects are moved in a
More informationExact equations are first order DEs of the form M(x, y) + N(x, y) y' = 0 for which we can find a function f(x, φ(x)) so that
Section 2.6 Exact Equations (ONLY) Key Terms: Exact equations are first order DEs of the form M(x, y) + N(x, y) y' = 0 for which we can find a function f(x, φ(x)) so that The construction of f(x, φ(x))
More informationObjectives. Materials
Activity 13 Objectives Understand what a slope field represents in terms of Create a slope field for a given differential equation Materials TI-84 Plus / TI-83 Plus Graph paper Introduction One of the
More informationComputational Methods. Sources of Errors
Computational Methods Sources of Errors Manfred Huber 2011 1 Numerical Analysis / Scientific Computing Many problems in Science and Engineering can not be solved analytically on a computer Numeric solutions
More informationFour equations are necessary to evaluate these coefficients. Eqn
1.2 Splines 11 A spline function is a piecewise defined function with certain smoothness conditions [Cheney]. A wide variety of functions is potentially possible; polynomial functions are almost exclusively
More information=.1( sin(2πx/24) y) y(0) = 70.
Differential Equations, Spring 2017 Computer Project 2, Due 11:30 am, Friday, March 10 The goal of this project is to implement the Euler Method and the Improved Euler Method, and to use these methods
More informationComputers in Engineering Root Finding Michael A. Hawker
Computers in Engineering COMP 208 Root Finding Michael A. Hawker Root Finding Many applications involve finding the roots of a function f(x). That is, we want to find a value or values for x such that
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationMath 104: Euler s Method and Applications of ODEs
Math 104: Euler s Method and Applications of ODEs Ryan Blair University of Pennsylvania Tuesday January 29, 2013 Ryan Blair (U Penn) Math 104:Euler s Method and Applications of ODEsTuesday January 29,
More informationPolynomial Approximation and Interpolation Chapter 4
4.4 LAGRANGE POLYNOMIALS The direct fit polynomial presented in Section 4.3, while quite straightforward in principle, has several disadvantages. It requires a considerable amount of effort to solve the
More informationMath 133A, September 10: Numerical Methods (Euler & Runge-Kutta) Section 1: Numerical Methods
Math 133A, September 10: Numerical Methods (Euler & Runge-Kutta) Section 1: Numerical Methods We have seen examples of first-order differential equations which can, by one method of another, be solved
More informationMath (Spring 2009): Lecture 5 Planes. Parametric equations of curves and lines
Math 18.02 (Spring 2009): Lecture 5 Planes. Parametric equations of curves and lines February 12 Reading Material: From Simmons: 17.1 and 17.2. Last time: Square Systems. Word problem. How many solutions?
More informationSuppose we want to find approximate values for the solution of the differential equation
A course in differential equations Chapter 2. Numerical methods of solution We shall discuss the simple numerical method suggested in the last chapter, and explain several more sophisticated ones. We shall
More informationThe Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
More information----- o Implicit Differentiation ID: A. dy r.---; d 2 Y 2. If- = '" 1-y- then - = dx 'dx 2. a c. -1 d. -2 e.
Name: Class: Date: ----- ID: A Implicit Differentiation Multiple Choice Identify the choice that best completes the statement or answers the question.. The slope of the line tangent to the curve y + (xy
More informationP1 REVISION EXERCISE: 1
P1 REVISION EXERCISE: 1 1. Solve the simultaneous equations: x + y = x +y = 11. For what values of p does the equation px +4x +(p 3) = 0 have equal roots? 3. Solve the equation 3 x 1 =7. Give your answer
More informationCS 450 Numerical Analysis. Chapter 7: Interpolation
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 informationNew formulations of the semi-lagrangian method for Vlasov-type equations
New formulations of the semi-lagrangian method for Vlasov-type equations Eric Sonnendrücker IRMA Université Louis Pasteur, Strasbourg projet CALVI INRIA Nancy Grand Est 17 September 2008 In collaboration
More informationRendering. A simple X program to illustrate rendering
Rendering A simple X program to illustrate rendering The programs in this directory provide a simple x based application for us to develop some graphics routines. Please notice the following: All points
More informationCS 130. Scan Conversion. Raster Graphics
CS 130 Scan Conversion Raster Graphics 2 1 Image Formation Computer graphics forms images, generally two dimensional, using processes analogous to physical imaging systems like: - Cameras - Human visual
More informationRasterization, or What is glbegin(gl_lines) really doing?
Rasterization, or What is glbegin(gl_lines) really doing? Course web page: http://goo.gl/eb3aa February 23, 2012 Lecture 4 Outline Rasterizing lines DDA/parametric algorithm Midpoint/Bresenham s algorithm
More informationLinear and Quadratic Least Squares
Linear and Quadratic Least Squares Prepared by Stephanie Quintal, graduate student Dept. of Mathematical Sciences, UMass Lowell in collaboration with Marvin Stick Dept. of Mathematical Sciences, UMass
More information[Anton, pp , pp ] & [Bourne, pp ]
hapter 3 Integral Theorems [Anton, pp. 1124 1130, pp. 1145 1160] & [Bourne, pp. 195 224] First of all some definitions which we will need in the following: Definition 3.1. (a) A domain (region) is an open
More information2.2 Graphs Of Functions. Copyright Cengage Learning. All rights reserved.
2.2 Graphs Of Functions Copyright Cengage Learning. All rights reserved. Objectives Graphing Functions by Plotting Points Graphing Functions with a Graphing Calculator Graphing Piecewise Defined Functions
More informationSOME PROPERTIES OF TRIGONOMETRIC FUNCTIONS. 5! x7 7! + = 6! + = 4! x6
SOME PROPERTIES OF TRIGONOMETRIC FUNCTIONS PO-LAM YUNG We defined earlier the sine cosine by the following series: sin x = x x3 3! + x5 5! x7 7! + = k=0 cos x = 1 x! + x4 4! x6 6! + = k=0 ( 1) k x k+1
More informationIndex. C m (Ω), 141 L 2 (Ω) space, 143 p-th order, 17
Bibliography [1] J. Adams, P. Swarztrauber, and R. Sweet. Fishpack: Efficient Fortran subprograms for the solution of separable elliptic partial differential equations. http://www.netlib.org/fishpack/.
More informationObjectives. Materials
S Activity 4 Objectives Materials Understand what a slope field represents in terms of dy Create a slope field for a given differential equation T-84 Plus / T-83 Plus Graph paper ntroduction ntroduction
More informationLecture 8. Divided Differences,Least-Squares Approximations. Ceng375 Numerical Computations at December 9, 2010
Lecture 8, Ceng375 Numerical Computations at December 9, 2010 Computer Engineering Department Çankaya University 8.1 Contents 1 2 3 8.2 : These provide a more efficient way to construct an interpolating
More informationNumerical Integration
1 Numerical Integration Jim Van Verth Insomniac Games jim@essentialmath.com Essential mathematics for games and interactive applications Talk Summary Going to talk about: Euler s method subject to errors
More informationMath 21a Tangent Lines and Planes Fall, What do we know about the gradient f? Tangent Lines to Curves in the Plane.
Math 21a Tangent Lines and Planes Fall, 2016 What do we know about the gradient f? Tangent Lines to Curves in the Plane. 1. For each of the following curves, find the tangent line to the curve at the point
More informationORDINARY DIFFERENTIAL EQUATIONS
Page 1 of 22 ORDINARY DIFFERENTIAL EQUATIONS Lecture 5 Visualization Tools for Solutions of First-Order ODEs (Revised 02 February, 2009 @ 08:05) Professor Stephen H Saperstone Department of Mathematical
More information1) Find. a) b) c) d) e) 2) The function g is defined by the formula. Find the slope of the tangent line at x = 1. a) b) c) e) 3) Find.
1 of 7 1) Find 2) The function g is defined by the formula Find the slope of the tangent line at x = 1. 3) Find 5 1 The limit does not exist. 4) The given function f has a removable discontinuity at x
More information