Lecture 21. Numerical differentiation. f ( x+h) f ( x) h h
|
|
- Charla Richardson
- 5 years ago
- Views:
Transcription
1 Lecture Numerical differentiation Introduction We can analytically calculate te derivative of any elementary function, so tere migt seem to be no motivation for calculating derivatives numerically. However we may need to estimate te derivative of a numerical function, or we may only ave a fixed set of sampled function values. In tese cases we need to estimate te derivative numerically. Finite difference approximation Te definition of te derivative of a function f ( x ) tat you will most often find in calculus textbooks is f ( x+) f (x ) df lim dx () Tis immediately suggests te approximation f ( x+) f ( x) df Δ f dx Δ x () were te step size is small but not zero. Tis is called a finite difference. Specifically it's a forward difference because we compare te function value at x wit its value at a point forward of tis along te x axis, x+. How small sould be? Because of round-off error, smaller is not always better. Let's use Scilab to estimate d x e dx x e e (3) for various values. Te absolute error in te estimate vs. is graped in Fig.. As decreases from down to 8 te error decreases also. However for 9 and smaller te error actually increases! Te culprit is round-off error in te form of te small difference of large numbers. Double precision aritmetic provides about 6 digits of precision. If 6 ten e e to 6 digits and te difference e e will be very inaccurate. Wen 8 te difference e e will be accurate to about 8 digits or about 8, te point at wic teoretical improvement in numerical accuracy is offset by iger round-off error. We typically ignore round-off error wen estimating numerical accuracy, but round-off error needs to be kept in mind wen implementing any algoritm. Let's investigate ow te numerical accuracy of our estimate varies wit step size. Assume we want to estimate te derivative of f (x ) at x. Write f as a power series f (x) f ()+ n EE Numerical Computing (n ) f () x n n! (4) 5-8-8
2 Lecture : Numerical differentiation / x Fig. : Error in forward-difference estimation of derivative of e at x vs. step size. 8 As decreases from to, numerical error decreases proportionally. As decreases furter, round-off error begins to exceed numerical error. Ten f () f () (n) f () n f ()+ f ()+ f () + n! 6 n (5) For small terefore f () f () f ()+ f () + (6) and we say tat te approximation is first-order accurate since te error (te second term on te rigt) varies as te first power of. Decreasing by a factor / will decrease te error by /. However, as we see in Fig. tis is only true up to te point tat round-off error begins to be significant. For double precision 8 is optimal. We can sift (6) along te x axis and rearrange to obtain te general forward-difference formula f (x ) f ( x+) f ( x) +O() (7) 3 Higer-order formulas Te forward-difference approximation (6) uses two samples of te function, namely f () and f (). Using tree samples we migt be able to get a better estimate of f (). Suppose we ave te samples f ( ), f (), f (). In terms of te power series representation of te function tese are EE Numerical Computing 5-8-8
3 Lecture : Numerical differentiation 3/ f ( ) f ( ) f () + f () f ( ) f () f ( )+ f ( ) + f () f (3) ( ) 3+ f (4) () 4 f (5 )( )5+ (8) f () + f (3) () 3+ f (4 ) ()4 + f (5) ()5+ Let's use a linear combination of tese values to estimate f () as f ()a f ( )+b f ( )+c f ( ) (9) were a, b, c are unknown coefficients tat we will coose to get te best possible estimate. Te sum a f ( )+b f ()+c f ( ) will include a term wit a factor of f (). We want tis to vanis. Tis requires (a+b+c) f ( ) a+b+c () wic is one equation in tree unknowns. Terms wit a factor of f () sould combine to give te derivative value f (). Tis requires ( a+c) f () f () a+c () We now ave two equations in tree unknowns. To get a tird equation we can require te next term, wic contains a factor of f (), to vanis. Tis gives us te equation (a+c) f ( ) a+c () Our tree equations in tree unknowns can be written as ( )( ) ( ) a b c Te solution is a, b, c and our approximation reads f () f ( ) f ( ) Using (8) we ave f () f ( ) f ()+ f (3) () + 6 so tis approximation is second-order accurate. Decreasing by a factor of / sould decrease te numerical error by a factor of /. Rearranging and writing tis for an arbitrary value of x we ave te formula EE Numerical Computing 5-8-8
4 Lecture : Numerical differentiation 4/ x Fig. : Error in central-difference estimate of derivative of e vs.. f (x ) f ( x+) f ( x ) +O( ) (3) Tis type of finite difference is called a central difference since it uses bot te forward sample f (x+) and te backward sample f (x ). Scilab code is given in te Appendix. Te error in te central-difference approximation d x e dx x e e is plotted in Fig.. Note ow te error reduces more rapidly wit decreasing. Tis allows te approximation to reac a greater accuracy before round-off error starts to become significant. Wit 5 te error is only about. We extend tis idea by using even more function samples. If we ave te five samples f ( ), f ( ), f (), f ( ), f ( ) we can form an estimate f ()a f ( )+b f ( )+c f ()+d f ()+e f ( ) (4) Tis as five unknowns, so we need to form five equations. In terms of te Taylor series representation of f ( x ) our five samples ave te form EE Numerical Computing 5-8-8
5 Lecture : Numerical differentiation f ( ) 8 (3) 6 (4 ) 3 (5) f () 3+ f ()4 f () 5+ f ( ) f () + f () f (3) () 3+ f (4) () 4 f (5) ( ) 5+ (5) f ( ) f ( )+ f () + f () + f (3) () 3+ f (4) () 4 + f (5) () (4) 3 (5) f ( )+ f () + f ( ) + f (3) ()3+ f () 4+ f () 5+ f ( ) f () f ( ) f ( ) f ()+ f () f () 5/ To get te f () terms in (4) to vanis requires (a+b+c+d +e) f () a+b+c+d +e (6) To get te f () terms to produce te value f () requires ( a b+d + e ) f () f ( ) a+b+d + e (7) Te remaining tree equations are obtained by requiring te f (), f (3) () and f (4 )( ) terms to vanis: ( ) b d a+ + + e f () 4 a+b+d +4 e ( 8 a b+d +8 e ) f (3) ()3 8 a b+d +8 e ( 6 a+b+d +6 e ) f (4) () 4 6 a+b+d +6 e ( 6+b+d +6 e ) f (4) ( ) 4 6 a+b+d +6 e Our five equations in five unknowns form te system ( )( ) ( ) a b c d e (8) wic as te solution a 8 8, b, c, d, e (9) Our approximation is terefore f () EE Numerical Computing f ( ) 8 f ( )+8 f () f ( ) () 5-8-8
6 Lecture : Numerical differentiation 6/ x Fig. 3: Fourt-order accurate central-difference approximation to te derivative of e. Using (5) te f (5) () terms are f (5) 5 () ( ) f (5) () 4 3 () so f ( ) 8 f ( )+8 f () f ( ) (5) 4 f () f () + 3 () We see tat tis is a fourt-order accurate approximation. Decreasing by a factor of / sould decrease te numerical error by a factor of /,. Tis is illustrated in Fig. 3 for x 3 4 f (x )e. For te error is only about. For arbitrary x value our fourt-order central difference approximation is f ( x) f (x ) 8 f ( x )+8 f ( x+) f ( x+ ) 4 +O( ) (3) Scilab code is given in te Appendix. 4 Te numderivative function (Scilab) Scilab as a numerical derivative function named numderivative. To calculate te derivative of te function f ( x ) at xx we execute fp numderivative(f,x); If f is a function several variables f (x), ten numderivative will return te gradient of f. It is also possible to specify te step size and te order of te approximation (, or 4). fp numderivative(f,x,,order); EE Numerical Computing 5-8-8
7 Lecture : Numerical differentiation 7/ Te default is second order (central difference) and Scilab cooses an optimal value of. Some examples: -->deff('yf(x)','yexp(-x)'); -->numderivative(f,) //default central difference wit optimal ans //exact value is -exp(-) >(f(.)-f())/. //forward difference. ans >numderivative(f,,.,) //forward difference. ans >(f(.)-f(.9))/. //central difference. ans >numderivative(f,,.,) //central difference. ans In most cases te default suffices. 5 Second derivative Wit reference to (8), if we want to approximate te second derivative f () as f ()a f ( )+b f ( )+c f ( ) (4) (a+b+c) f ( ) a+b+c (5) we would require ( a+c) f ( ) a+c (6) and (a+c) f () f () a+c (7) Te solution to tese tree equations is ac, b (8) and we find f () f ()+ f ( ) f ()+ (4) f ( ) + (9) so te approximation is second-order accurate. For arbitrary x value we ave EE Numerical Computing 5-8-8
8 Lecture : Numerical differentiation 8/ f (x) f ( x+) f ( x)+ f ( x ) +O ( ) (3) 6 Partial derivatives For a function f (x, y ) te partial derivative can be defined as x f ( x+, y ) f (x, y) lim x tat is, we old y fixed and compute te derivative as if f was only a function of x. A centraldifference approximation is f ( x+, y) f ( x, y) x (3) f ( x, y+) f ( x, y ) y (3) f f ( x+, y) f ( x, y )+ f ( x, y ) x (33) Likewise and f Mixed partial derivative approximations suc as can be developed in steps suc as y x f y x [ ] [ ] x y+ x y f (x+, y+) f ( x, y+) f (x+, y ) f ( x, y ) f ( x+, y+) f (x, y+) f (x+, y )+ f ( x, y ) 4 (34) 7 Differential equations 7. Ordinary differential equations An ordinary differential equation (ODE) relates a single independent variable, e.g., x, to a function f (x ) and its derivatives f (x ), f (x),. Most pysical laws are expressed in terms of differential equations, ence teir great importance. Certain classes of ODEs can be solved analytically but many cannot. In eiter case our derivative formulas can be used to develop numerical solutions. Suppose a pysical problem is described by a differential equation of te form EE Numerical Computing 5-8-8
9 Lecture : Numerical differentiation 9/ f + f +7 f (35) x (36) One can verify tat f ( x )e cos(4 x) solves (35) by taking derivatives and substituting into te equation. A numerical approximation to (35) is given by (using (3) and (3)) f ( x+) f ( x)+ f ( x ) f (x+) f ( x ) + +7 f ( x) (37) Solving tis for f (x+) we obtain f (x+) ( 7 ) f (x ) ( ) f ( x ) + (38) Let's use tis to calculate f (x ) for x,,, 3,. To get started we need te first two values f ( x ), f ( x )e cos (4 ) (39) Ten we can apply (38) to get f (x 3x +), f ( x 4 x3 +) and so on as long as we wis. In Scilab tis looks someting like.; x ::5; f() ; f() exp(-)*cos(4*); for i:n- y(i+) ((-7*^)*y(i)-(-)*y(i-))/(+); end Te resulting numerical solution and te exact solution are sown in Fig. 4. Te agreement is excellent. Function odecentdiff in te Appendix uses tis idea to numerically solve a second-order equation of te form y + p( x ) y +q ( x) yr ( x) Given and initial x value x, a step size and te two function values y ( x ) and y (x +). Fig. 4 compares te numerical solution of (35) using odecentdiff wit te exact solution y f ( x )e x cos (4 x ) (4) for a step size.. 7. Partial differential equations A partial differential equation (PDE) relates two or more independent variables, e.g., x,y, to a f f f f,,,,. One of te most function f ( x, y ) and its partial derivatives x y x x y important PDEs is Laplace's equation EE Numerical Computing 5-8-8
10 Lecture : Numerical differentiation / Fig. 4: Solid line: numerical solution of (35); circles: exact solution. f f + x y (4) Numerically we can write f f + x y f (x+, y ) f ( x, y )+ f ( x, y ) f ( x, y+) f (x, y)+ f ( x, y ) (4) + f ( x+, y)+ f ( x, y)+ f (x, y+)+ f ( x, y ) 4 f ( x, y) Te last expression is zero wen f (x, y ) [ f ( x+, y )+ f (x, y )+ f ( x, y +)+ f ( x, y ) ] 4 (43) wic relates te value f (x, y ) to its neigboring values. Specifically f (x, y ) is equal to te average of its neigbor's values. EE Numerical Computing 5-8-8
11 Lecture : Numerical differentiation / 8 Appendix Scilab code 8. nd order central difference // derivsecondorder.sci // 4--5,, for pedagogic purposes // Numerical estimation of derivative of f(x) using nd-order // accurate central difference and "optimum" step size. function ypderivsecondorder(f, x) e-5*(+abs(x)); //step size scales wit x, no less tan e-5 yp (f(x+)-f(x-))/(*); endfunction 8. 4t order central difference // derivfourtorder.sci // 4--5,, for pedagogic purposes // Numerical estimation of derivative of f(x) using 4t-order // accurate central difference and "optimum" step size. function ypderivfourtorder(f, x) e-3*(+abs(x)); //step size scales wit x, no less tan e-3 yp (f(x-*)-8*f(x-)+8*f(x+)-f(x+*))/(*); endfunction 8.3 Differential equation solver using nd order central difference // odecentdiff.sci // 4--5,, for pedagogic purposes // Uses nd-order accurate central difference approximation to // derivatives to solve ode y''+p(x)y'+q(x)yr(x) // approximations are // y' (y(x+)-y(x-))/() and y'' (y(x+)-y(x)+y(x-))/^ // p,q,r are functions, x is te initial x value, is step size, // n is number of points to solve for, yy(x), yy(x+). function [x, y]odecentdiff(p, q, r, x,, n, y, y) x zeros(n,); y zeros(n,); x() x; x() x()+; y() y; y() y; *; for i:n- p *p(x(i)); x(i+) x(i)+; y(i+) (**r(x(i))+(4-**q(x(i)))*y(i)+(p-)*y(i-))/(+p); end endfunction EE Numerical Computing 5-8-8
NUMERICAL DIFFERENTIATION
NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationThe total error in numerical differentiation
AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and
More informationHigher Derivatives. Differentiable Functions
Calculus 1 Lia Vas Higer Derivatives. Differentiable Functions Te second derivative. Te derivative itself can be considered as a function. Te instantaneous rate of cange of tis function is te second derivative.
More informationSection 15.6 Directional Derivatives and the Gradient Vector
Section 15.6 Directional Derivatives and te Gradient Vector Finding rates of cange in different directions Recall tat wen we first started considering derivatives of functions of more tan one variable,
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
More informationMath 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0
3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
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 informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationSection 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is
Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit
More informationDerivatives of Exponentials
mat 0 more on derivatives: day 0 Derivatives of Eponentials Recall tat DEFINITION... An eponential function as te form f () =a, were te base is a real number a > 0. Te domain of an eponential function
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationf a h f a h h lim lim
Te Derivative Te derivative of a function f at a (denoted f a) is f a if tis it exists. An alternative way of defining f a is f a x a fa fa fx fa x a Note tat te tangent line to te grap of f at te point
More informationThe Derivative as a Function
Section 2.2 Te Derivative as a Function 200 Kiryl Tsiscanka Te Derivative as a Function DEFINITION: Te derivative of a function f at a number a, denoted by f (a), is if tis limit exists. f (a) f(a + )
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More information2.11 That s So Derivative
2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point
More informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationDefinition of the Derivative
Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationQuaternion Dynamics, Part 1 Functions, Derivatives, and Integrals. Gary D. Simpson. rev 01 Aug 08, 2016.
Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals Gary D. Simpson gsim1887@aol.com rev 1 Aug 8, 216 Summary Definitions are presented for "quaternion functions" of a quaternion. Polynomial
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationExam 1 Review Solutions
Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),
More informationMath 1241 Calculus Test 1
February 4, 2004 Name Te first nine problems count 6 points eac and te final seven count as marked. Tere are 120 points available on tis test. Multiple coice section. Circle te correct coice(s). You do
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationFirst we will go over the following derivative rule. Theorem
Tuesday, Feb 1 Tese slides will cover te following 1 d [cos(x)] = sin(x) iger-order derivatives 3 tangent line problems 4 basic differential equations First we will go over te following derivative rule
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationTest 2 Review. 1. Find the determinant of the matrix below using (a) cofactor expansion and (b) row reduction. A = 3 2 =
Test Review Find te determinant of te matrix below using (a cofactor expansion and (b row reduction Answer: (a det + = (b Observe R R R R R R R R R Ten det B = (((det Hence det Use Cramer s rule to solve:
More informationIntegral Calculus, dealing with areas and volumes, and approximate areas under and between curves.
Calculus can be divided into two ke areas: Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and minima problems Integral
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Matematics and Computer Science. ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI3097 TU) Tuesday January 9 008, 9:00-:00
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More information1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist
Mat 1120 Calculus Test 2. October 18, 2001 Your name Te multiple coice problems count 4 points eac. In te multiple coice section, circle te correct coice (or coices). You must sow your work on te oter
More information1 Introduction to Optimization
Unconstrained Convex Optimization 2 1 Introduction to Optimization Given a general optimization problem of te form min x f(x) (1.1) were f : R n R. Sometimes te problem as constraints (we are only interested
More informationSolution. Solution. f (x) = (cos x)2 cos(2x) 2 sin(2x) 2 cos x ( sin x) (cos x) 4. f (π/4) = ( 2/2) ( 2/2) ( 2/2) ( 2/2) 4.
December 09, 20 Calculus PracticeTest s Name: (4 points) Find te absolute extrema of f(x) = x 3 0 on te interval [0, 4] Te derivative of f(x) is f (x) = 3x 2, wic is zero only at x = 0 Tus we only need
More informationMath 161 (33) - Final exam
Name: Id #: Mat 161 (33) - Final exam Fall Quarter 2015 Wednesday December 9, 2015-10:30am to 12:30am Instructions: Prob. Points Score possible 1 25 2 25 3 25 4 25 TOTAL 75 (BEST 3) Read eac problem carefully.
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More informationPolynomial Functions. Linear Functions. Precalculus: Linear and Quadratic Functions
Concepts: definition of polynomial functions, linear functions tree representations), transformation of y = x to get y = mx + b, quadratic functions axis of symmetry, vertex, x-intercepts), transformations
More informationReflection Symmetries of q-bernoulli Polynomials
Journal of Nonlinear Matematical Pysics Volume 1, Supplement 1 005, 41 4 Birtday Issue Reflection Symmetries of q-bernoulli Polynomials Boris A KUPERSHMIDT Te University of Tennessee Space Institute Tullaoma,
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationAverage Rate of Change
Te Derivative Tis can be tougt of as an attempt to draw a parallel (pysically and metaporically) between a line and a curve, applying te concept of slope to someting tat isn't actually straigt. Te slope
More informationA.P. CALCULUS (AB) Outline Chapter 3 (Derivatives)
A.P. CALCULUS (AB) Outline Capter 3 (Derivatives) NAME Date Previously in Capter 2 we determined te slope of a tangent line to a curve at a point as te limit of te slopes of secant lines using tat point
More informationSection 3: The Derivative Definition of the Derivative
Capter 2 Te Derivative Business Calculus 85 Section 3: Te Derivative Definition of te Derivative Returning to te tangent slope problem from te first section, let's look at te problem of finding te slope
More informationMath 312 Lecture Notes Modeling
Mat 3 Lecture Notes Modeling Warren Weckesser Department of Matematics Colgate University 5 7 January 006 Classifying Matematical Models An Example We consider te following scenario. During a storm, a
More informationTHE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225
THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:
More informationFinding and Using Derivative The shortcuts
Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex
More informationMATH 155A FALL 13 PRACTICE MIDTERM 1 SOLUTIONS. needs to be non-zero, thus x 1. Also 1 +
MATH 55A FALL 3 PRACTICE MIDTERM SOLUTIONS Question Find te domain of te following functions (a) f(x) = x3 5 x +x 6 (b) g(x) = x+ + x+ (c) f(x) = 5 x + x 0 (a) We need x + x 6 = (x + 3)(x ) 0 Hence Dom(f)
More information1watt=1W=1kg m 2 /s 3
Appendix A Matematics Appendix A.1 Units To measure a pysical quantity, you need a standard. Eac pysical quantity as certain units. A unit is just a standard we use to compare, e.g. a ruler. In tis laboratory
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationThe Derivative The rate of change
Calculus Lia Vas Te Derivative Te rate of cange Knowing and understanding te concept of derivative will enable you to answer te following questions. Let us consider a quantity wose size is described by
More informationChapter XI. Solution of Ordinary Differential Equations
dy d 7 Capter XI Solution of Ordinary Differential Equations Several approaces are described for te solution of ODE (ordinary differential equations). Tese include: Te Taylor series Metod Te Euler Metod
More informationINTRODUCTION TO CALCULUS LIMITS
Calculus can be divided into two ke areas: INTRODUCTION TO CALCULUS Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and
More informationFractional Derivatives as Binomial Limits
Fractional Derivatives as Binomial Limits Researc Question: Can te limit form of te iger-order derivative be extended to fractional orders? (atematics) Word Count: 669 words Contents - IRODUCIO... Error!
More informationExponentials and Logarithms Review Part 2: Exponentials
Eponentials and Logaritms Review Part : Eponentials Notice te difference etween te functions: g( ) and f ( ) In te function g( ), te variale is te ase and te eponent is a constant. Tis is called a power
More informationLab 6 Derivatives and Mutant Bacteria
Lab 6 Derivatives and Mutant Bacteria Date: September 27, 20 Assignment Due Date: October 4, 20 Goal: In tis lab you will furter explore te concept of a derivative using R. You will use your knowledge
More informationName: Answer Key No calculators. Show your work! 1. (21 points) All answers should either be,, a (finite) real number, or DNE ( does not exist ).
Mat - Final Exam August 3 rd, Name: Answer Key No calculators. Sow your work!. points) All answers sould eiter be,, a finite) real number, or DNE does not exist ). a) Use te grap of te function to evaluate
More informationPhysics 121, April 1, Equilibrium. Physics 121. April 1, Physics 121. April 1, Course Information. Discussion of Exam # 2
Pysics 121, April 1, 2008. Pysics 121. April 1, 2008. Course Information Discussion of Exam # 2 Topics to be discussed today: Requirements for Equilibrium Gravitational Equilibrium Sample problems Pysics
More informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
More information4.2 - Richardson Extrapolation
. - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence
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 informationPreface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.
Preface Here are my online notes for my course tat I teac ere at Lamar University. Despite te fact tat tese are my class notes, tey sould be accessible to anyone wanting to learn or needing a refreser
More informationDifferential Calculus (The basics) Prepared by Mr. C. Hull
Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit
More information3.1 Extreme Values of a Function
.1 Etreme Values of a Function Section.1 Notes Page 1 One application of te derivative is finding minimum and maimum values off a grap. In precalculus we were only able to do tis wit quadratics by find
More informationSECTION 1.10: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES
(Section.0: Difference Quotients).0. SECTION.0: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES Define average rate of cange (and average velocity) algebraically and grapically. Be able to identify, construct,
More informationExercises for numerical differentiation. Øyvind Ryan
Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can
More informationFall 2014 MAT 375 Numerical Methods. Numerical Differentiation (Chapter 9)
Fall 2014 MAT 375 Numerical Metods (Capter 9) Idea: Definition of te derivative at x Obviuos approximation: f (x) = lim 0 f (x + ) f (x) f (x) f (x + ) f (x) forward-difference formula? ow good is tis
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More informationSection 2: The Derivative Definition of the Derivative
Capter 2 Te Derivative Applied Calculus 80 Section 2: Te Derivative Definition of te Derivative Suppose we drop a tomato from te top of a 00 foot building and time its fall. Time (sec) Heigt (ft) 0.0 00
More information1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow
1.7, Groundwater Hydrology Prof. Carles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow Simulation: Te prediction of quantities of interest (dependent variables) based upon an equation
More information(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a?
Solutions to Test 1 Fall 016 1pt 1. Te grap of a function f(x) is sown at rigt below. Part I. State te value of eac limit. If a limit is infinite, state weter it is or. If a limit does not exist (but is
More information(4.2) -Richardson Extrapolation
(.) -Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Suppose tat lim G 0 and lim F L. Te function F is said to converge to L as
More informationThe Priestley-Chao Estimator
Te Priestley-Cao Estimator In tis section we will consider te Pristley-Cao estimator of te unknown regression function. It is assumed tat we ave a sample of observations (Y i, x i ), i = 1,..., n wic are
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationWYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationlecture 35: Linear Multistep Mehods: Truncation Error
88 lecture 5: Linear Multistep Meods: Truncation Error 5.5 Linear ultistep etods One-step etods construct an approxiate solution x k+ x(t k+ ) using only one previous approxiation, x k. Tis approac enoys
More informationMATH1901 Differential Calculus (Advanced)
MATH1901 Dierential Calculus (Advanced) Capter 3: Functions Deinitions : A B A and B are sets assigns to eac element in A eactl one element in B A is te domain o te unction B is te codomain o te unction
More information, meant to remind us of the definition of f (x) as the limit of difference quotients: = lim
Mat 132 Differentiation Formulas Stewart 2.3 So far, we ave seen ow various real-world problems rate of cange and geometric problems tangent lines lead to derivatives. In tis section, we will see ow to
More information1 2 x Solution. The function f x is only defined when x 0, so we will assume that x 0 for the remainder of the solution. f x. f x h f x.
Problem. Let f x x. Using te definition of te derivative prove tat f x x Solution. Te function f x is only defined wen x 0, so we will assume tat x 0 for te remainder of te solution. By te definition of
More informationThese errors are made from replacing an infinite process by finite one.
Introduction :- Tis course examines problems tat can be solved by metods of approximation, tecniques we call numerical metods. We begin by considering some of te matematical and computational topics tat
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationMA455 Manifolds Solutions 1 May 2008
MA455 Manifolds Solutions 1 May 2008 1. (i) Given real numbers a < b, find a diffeomorpism (a, b) R. Solution: For example first map (a, b) to (0, π/2) and ten map (0, π/2) diffeomorpically to R using
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
More informationSolutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014
Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.
More informationAnalytic Functions. Differentiable Functions of a Complex Variable
Analytic Functions Differentiable Functions of a Complex Variable In tis capter, we sall generalize te ideas for polynomials power series of a complex variable we developed in te previous capter to general
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
More informationMATH1151 Calculus Test S1 v2a
MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums
More information