1 The Euler Forward scheme (schéma d'euler explicite)
|
|
- Roxanne Fox
- 6 years ago
- Views:
Transcription
1 TP : Finite dierence metod for a European option M2 Modélisation aléatoire - Université Denis-Diderot Cours EDP en Finance et Métodes Numériques. December Te Euler Forward sceme (scéma d'euler explicite) We look for a numerical approximation of te European put function v = v(t, s), t [0, T ], s [0, S max ]. It satises in rst approximation te Black and Scoles backward PDE on te truncated domain Ω = [S min, S max ]: s2 2 s 2 v rs s v + rv = 0, t (0, T ), s (S min, S max ) t v σ2 2 v(t, S min ) = v l (t) Ke rt S min, t (0, T ) v(t, S max ) = v r (t) 0, t (0, T ) v(0, s) = ϕ(s) := (K s) +, s (S min, S max ). We will consider te following parameters: K = 100, S min = 0, S max = 200, T = 1, σ = 0.2, r = 0.1 In particular, we aim at computing v(t, s) at nal time t = T. We rst introduce a discrete mes as follows. Let := Smax S min I+1 be (spatial) mes step, and t := T N be te time step. Ten s := S min +, = 0,..., I + 1 are te mes points, and t n = n t, n = 0,..., N te time mes. We are looking for U n, an approximation of v(t n, s ). For any function v C 2 (or v C 3 for (4)), we recall te following approximations, as 0, v (s ) = v(s ) v(s 1 ) v (s ) = v(s +1) v(s ) v (s ) = v(s +1) v(s 1 ) 2 (1) + O() (2) + O() (3) + O( 2 ). (4) We terefore obtain several possible approximations by nite dierences for te rst order derivative: S v(t n, s ) U n U n 1 S v(t n, s ) U n +1 U n S v(t n, s ) U n +1 U n 1 2 (backward dierence approximation) (5) (forward dierence approximation) (6) (centered dierence) (7) Te rst two approximations are said to be consistent of order 1 (in space), wile te second one is consistent of order 2. 1
2 We also recall te approximation 2 SSv(t n, s ) U n 1 + 2U n U n +1 2, (8) wic is consitent of order 2 in space. Hence we obtain te so-called "Euler Forward sceme" (or Explicit Euler sceme), abreviated "EE" using te centered approximation, as follows: U n+1 U n t + σ2 2 s2 U n 1 + 2U n U n +1 2 rs U n +1 U n ru n = 0 n = 0,..., N 1, = 1,..., I U n 0 = v l (t n ) Ke rtn S min, n = 0,..., N (9) U n I+1 = v r (t n ) 0, U 0 = ϕ(s ) (K s ) +, n = 0,..., N = 1,..., I Let us remark tat we ave taken = 1 and = I as extremal indices in. For = 1, te sceme utilizes te known value U0 n := v l(t n ) (left boundary value). For = I, te sceme utilizes te known value UI+1 n := v r(t n ) (rigt boundary value). 2 Programming Euler Forward 2.1 Preliminaries. We coose to work wit te unkown te vector corresponding to (v(t n, s )) =1,...,I : U n = U n 1. U n I. We would like to write (9) under te vector form as follows: U n+1 U n + AU n + q(t n ) = 0, (10) t were A is a square matrix of dimension I and q(t) is a column vector of size I. Let us denote α := σ2 s 2 2 2, β := r s 2. We look for A and q(t) suc tat ( α i + β i )U n i 1 + (2α i + r)u n i + ( α i β i )U n i+1 (AU + q(t n )) i. By identication we see tat A is a tridiagonal matrix 2α 1 + r α 1 β 1 0 α 2 + β 2 2α 2 + r α 2 β A := α i + β i 2α i + r α i β i α I + β I 2α I + r 2
3 and q(t) will contain te known boundary values U 0 and U I+1 : ( α 1 + β 1 )v l (t) q(t) := ( α I β I )v r (t) 2.2 Getting into te program. a) Download te working program 1 (tp1.m if using Matlab, or tp1.sci if using Scilab). Tis program as some lines to be completed in order to work properly. In te case of te matlab program, you will also ave to download oter function les BS.m and ploot.m (you can also take te solution le tp1sol.m) b) Ceck te program for te payo function u0, ul (for v l ) and ur (for v r ) and program tem correctly. (notice tat te instruction y=max(k-s,0) works for s vector, K scalar, and return a vector of same size as s.) c) In te section MESH, complete te value of te mes step and program te vector s containing te (s ) 1 I values. Typically it sould look like: dt=t/n; =(Smax-Smin)/(I+1); s=smin+ (1:I)'*; // column vector (s 1,..., s I ) T d) Program te matrix A and test te program. Under te matlab command window, type tp1. (Under Scilab, type exec tp1.sci or exec('tp1.sci',-1)) e) Program te function q(t). For instance in matlab: q [(-alpa(1) + bet(1))* ul(t); zeros(i-2,1); (-alpa(end) - bet(end))* ur(t)]; In Scilab tis may look like function y=ul(t); y=k*exp(-r*t); endfunction function y=ur(t); y=0; endfunction; function y=q(t) y=zeros(s); // vector of zeros wit same size as te s vector y(1)= (-a(1)+b(1))*ul(t); u($)= (-a($)-b($))*ur(t); endfunction 2.2 Euler Forward sceme (or "Euler Explicit" sceme) a) Program te explicit form of te vector U n+1, in terms of U n, in te main loop, using te matrix A and te function q. In te end, it sould look like (in matlab) case 'EE' P = (Id - dt*a)*p - dt*q(t); 1 See ttps://lll.mat.upmc.fr/bokanowski/enseignement/ens2016.tml 3
4 Note tat te grapic function ploot also plots te exact Black and Scoles formula (see te given function BS in te le BS.m). b) Correct te lign errli=0 in te main loop, in order to compute correctly te maximum norm. between te sceme values and te Black and Scoles values: U V BS = max 1 i I U i V BS (s i ) 2.3 Solution. A solution le is given in tp1sol.m 3 First numerical tests a) Test te Euler forward sceme (EE). First x N = 10 and take I = 10, 20, 50,.... Ten take te following N = I values : 10, 20, 50, 100. Observe tat: - te sceme is not always numerically stable - it does not always give a positive solution (i.e. we do not always ave U n 0 n, ). b) Try to understand te origin of te oscillations wen tey occur. One sould look at te "amplication" matrix dened as B := I d ta (Te coecients of B are tey positive? Do tey ave a modulus smaller tan 1?) c) Fill in te CFL number dened ere as µ := t 2 /( 1 σ 2 S max 2 ) and print it. Ceck tat tere is no stability problem wen ν is suciently small. d) In te case wen σ = and for instance wit N = 10 and I = 50: observe tat te sceme is stable but tat te solution given by te sceme is not positive. e) Order of te sceme. We rst consider te following values of I and of N: I = 10, 20, 40, 80, 160, and N = I 2 /10 (tat is, N = 10, 40,...). Tis is in order to ave t 2. Fill in te following error table and compute te corresponding (spatial) order of te metod. I N e k order α k More precisely, if te error is e k for a given parameter I = I k, we will compute te order at step k by te formula α k := log(e k 1/e k ) log( k 1 / k ) log(e k 1/e k ). log(2) Te idea is to try to detect a beavior of te form e k = C α k, were C est a constant and were k is te spatial mes size corresponding to I k. 4
5 Te numerical (spatial) order sould be close to two. But ere setting N const I 2 is costly in terms of number of operations. More precisely, since t = const 2 one can see tat error also beaves as O( t), tat is, a rst order error beavior wit respect to te time discretisation. Tis motivates te use of implicit scemes in order to avoid te time-step condition (or "CFL" condition). Exercice: backward and forward dierence approximations a) Program in te same way te oter approximations (backward dierences / décentrage droit : 'DROIT' and te forward dierences - décentrage gauce 'GAUCHE'). Te parameter CENTRAGE is dened at te begining of te program (section NUMERICAL DATA) and determines te type of nite dierence approximation tat is used. b) Using te forward dierences, x temporarily σ = and N = 10, I = 50, and ceck tat now te solution keeps positive. c) Using te backward dierences, x temporarily σ = and N = 10, I = 50, and ceck tat now te numerical solution is unstable (understand te problem by looking at te coecients of te amplication matrix). 4 Implicit Euler sceme (scéma d'euler implicite) Te CFL constraint imposes some restriction on te time step t. Implicit scemes may allow us to get rid of tis restriction. For instance, te implicit euler sceme, ereafter abrieviated "EI" (for Frenc "Euler Implicit"), wit centered dierence approximation for te rst spatial derivatives reads U n+1 U n t + σ2 2 s2 U n+1 1 n+1 + 2U U+1 n+1 U+1 n+1 2 rs U 1 n ru n+1 = 0 n = 0,..., N 1, = 1,..., I U0 n+1 = v l (t n+1 ) Ke rt n+1 S min, n = 0,..., N 1 (11) U n+1 I+1 = v r(t n+1 ) 0, n = 0,..., N 1 U 0 = (K s ), = 1,..., I a) Write te sceme in vector form, in analogy wit (10). b) Program EI: set te parameter SCHEMA='EI' at te begining of te program, and complete accordingly, in te main loop, te part case 'EI'... Note tat in order to solve a linear system of te form Ax = b one may use te linear solver in matlab (resp. Scilab) as follows x=a\b; c) Ceck tat wit te EI sceme tere is no more stability problems (wit for instance N = 10 and I = 50). 5 Crank-Nicolson sceme a) Program te Crank-Nicolson sceme (CN) (tis is te θ-sceme wit θ = 1 2 ). 5
6 b) Numerically study te convergence of te centered CN sceme. To do so, one can compute te numerical error obtained wit te following values N = I + 1 {10, 20, 40, 80, 160}. Ten do te same study as for te Euler implicit sceme and compare by making an error table (and wit an order estimation) One sould observe a second order beavior for CN, wile it sould be rst order for EI. c) Ceck tat tese results are coerent wit te teory. 6 Exercices Exercice. 1 ("Call option") a) Using te put-call parity formula, do a function tat also computes te black and scoles formula for te call. b) Propose adequate left and rigt boundary conditions (tat is, at s = S min and s = S max ) for te call option, in te form (v(t, S min ) = v l (t) and v(t, S max ) = v r (t), were v l, v r are function to determine. 2 c) Propose an apropriate PDE for te call option, using tese boundary conditions. d) Write an Euler explicit sceme for te call option, and program it. Exercice. 2 (Improved eciency using sparse matrices) Te matrix A as only a few non zero elements (about 3I non zero elements). It is possible to code only te non-zero elements of A by using te matrix type sparse. Type elp sparse or doc sparse for documentation in matlab. use A=spzeros(I,I) to initialize a sparse 0-matrix of size I I. type A or full(a) to ceck te values, sparse(a) corresponds to te sparse matrix of a full matrix A. use speye(i,i) to get te sparse identity matrix. Program te EI or CN sceme by using only sparse matrices, and compare te speed of te new code wit respect to te old one for large I values. Execution time is in general reduced. It can be evaluated by using te commands tic; toc; (resp. time() and time() in Scilab) as follows: tic; % INSTRUCTIONS; %... t=toc; printf('cpu time : t=%5.2f, t); (ten te variable t contains te time elapsed between te calls of tic and toc.) 2 v l (t, s) 0, v r(t, s) s Ke rt. 6
1 The Euler Forward scheme (schéma d Euler explicite)
TP : Fiite differece metod for a Europea optio M2 Modélisatio aléatoire - Uiversité Deis-Diderot Cours EDP e Fiace et Métodes Numériques December 8, 2017 1 Te Euler Forward sceme (scéma d Euler explicite)
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 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 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 21. Numerical differentiation. f ( x+h) f ( x) h h
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
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 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 information3.4 Worksheet: Proof of the Chain Rule NAME
Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are
More informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
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 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 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 informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
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 informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
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 informationSection 2.1 The Definition of the Derivative. We are interested in finding the slope of the tangent line at a specific point.
Popper 6: Review of skills: Find tis difference quotient. f ( x ) f ( x) if f ( x) x Answer coices given in audio on te video. Section.1 Te Definition of te Derivative We are interested in finding te slope
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 information1 Upwind scheme for advection equation with variable. 2 Modified equations: numerical dissipation and dispersion
1 Upwind sceme for advection equation wit variable coefficient Consider te equation u t + a(x)u x Applying te upwind sceme, we ave u n 1 = a (un u n 1), a 0 u n 1 = a (un +1 u n ) a < 0. CFL condition
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 informationCS522 - Partial Di erential Equations
CS5 - Partial Di erential Equations Tibor Jánosi April 5, 5 Numerical Di erentiation In principle, di erentiation is a simple operation. Indeed, given a function speci ed as a closed-form formula, its
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 informationParabolic PDEs: time approximation Implicit Euler
Part IX, Capter 53 Parabolic PDEs: time approximation We are concerned in tis capter wit bot te time and te space approximation of te model problem (52.4). We adopt te metod of line introduced in 52.2.
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 informationdu(l) 5 dl = U(0) = 1 and 1.) Substitute for U an unspecified trial function into governing equation, i.e. dx + = =
Consider an ODE of te form: Finite Element Metod du fu g d + wit te following Boundary Conditions: U(0) and du(l) 5 dl.) Substitute for U an unspecified trial function into governing equation, i.e. ^ U
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 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 information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More informationRobotic manipulation project
Robotic manipulation project Bin Nguyen December 5, 2006 Abstract Tis is te draft report for Robotic Manipulation s class project. Te cosen project aims to understand and implement Kevin Egan s non-convex
More informationNUMERICAL 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 informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationChapter 1 Functions and Graphs. Section 1.5 = = = 4. Check Point Exercises The slope of the line y = 3x+ 1 is 3.
Capter Functions and Graps Section. Ceck Point Exercises. Te slope of te line y x+ is. y y m( x x y ( x ( y ( x+ point-slope y x+ 6 y x+ slope-intercept. a. Write te equation in slope-intercept form: x+
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
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 informationMAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016
MAT244 - Ordinary Di erential Equations - Summer 206 Assignment 2 Due: July 20, 206 Full Name: Student #: Last First Indicate wic Tutorial Section you attend by filling in te appropriate circle: Tut 0
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 informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
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 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 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 information7 Semiparametric Methods and Partially Linear Regression
7 Semiparametric Metods and Partially Linear Regression 7. Overview A model is called semiparametric if it is described by and were is nite-dimensional (e.g. parametric) and is in nite-dimensional (nonparametric).
More informationPOLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY
APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky
More informationMath Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim
Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z
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 informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationMath 34A Practice Final Solutions Fall 2007
Mat 34A Practice Final Solutions Fall 007 Problem Find te derivatives of te following functions:. f(x) = 3x + e 3x. f(x) = x + x 3. f(x) = (x + a) 4. Is te function 3t 4t t 3 increasing or decreasing wen
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 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 informationChapter 1D - Rational Expressions
- Capter 1D Capter 1D - Rational Expressions Definition of a Rational Expression A rational expression is te quotient of two polynomials. (Recall: A function px is a polynomial in x of degree n, if tere
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 informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
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 informationTreatment of Discontinuities
Treatment of Discontinuities Today, we sall look at te problem of dealing wit discontinuities in models. Models from engineering often exibit discontinuities tat describe situations suc as switcing, limiters,
More informationSTABILITY RESULTS FOR SECOND ORDER BACKWARD DIFFERENTIATION SCHEMES FOR PARABOLIC HAMILTON-JACOBI-BELLMAN EQUATIONS
STABILITY RESULTS FOR SECOND ORDER BACKWARD DIFFERENTIATION SCHEMES FOR PARABOLIC HAMILTON-JACOBI-BELLMAN EQUATIONS OLIVIER BOKANOWSKI, ATHENA PICARELLI, AND CHRISTOPH REISINGER Abstract. We propose a
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationPractice Problem Solutions: Exam 1
Practice Problem Solutions: Exam 1 1. (a) Algebraic Solution: Te largest term in te numerator is 3x 2, wile te largest term in te denominator is 5x 2 3x 2 + 5. Tus lim x 5x 2 2x 3x 2 x 5x 2 = 3 5 Numerical
More informationSharp Korn inequalities in thin domains: The rst and a half Korn inequality
Sarp Korn inequalities in tin domains: Te rst and a alf Korn inequality Davit Harutyunyan (University of Uta) joint wit Yury Grabovsky (Temple University) SIAM, Analysis of Partial Dierential Equations,
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 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 information2.3 Product and Quotient Rules
.3. PRODUCT AND QUOTIENT RULES 75.3 Product and Quotient Rules.3.1 Product rule Suppose tat f and g are two di erentiable functions. Ten ( g (x)) 0 = f 0 (x) g (x) + g 0 (x) See.3.5 on page 77 for a proof.
More informationDifferential equations. Differential equations
Differential equations A differential equation (DE) describes ow a quantity canges (as a function of time, position, ) d - A ball dropped from a building: t gt () dt d S qx - Uniformly loaded beam: wx
More informationInf sup testing of upwind methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Met. Engng 000; 48:745 760 Inf sup testing of upwind metods Klaus-Jurgen Bate 1; ;, Dena Hendriana 1, Franco Brezzi and Giancarlo
More information. Compute the following limits.
Today: Tangent Lines and te Derivative at a Point Warmup:. Let f(x) =x. Compute te following limits. f( + ) f() (a) lim f( +) f( ) (b) lim. Let g(x) = x. Compute te following limits. g(3 + ) g(3) (a) lim
More informationBlanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS
Opuscula Matematica Vol. 26 No. 3 26 Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Abstract. In tis work a new numerical metod is constructed for time-integrating
More informationChapter 9 - Solved Problems
Capter 9 - Solved Problems Solved Problem 9.. Consider an internally stable feedback loop wit S o (s) = s(s + ) s + 4s + Determine weter Lemma 9. or Lemma 9. of te book applies to tis system. Solutions
More informationQuantum Numbers and Rules
OpenStax-CNX module: m42614 1 Quantum Numbers and Rules OpenStax College Tis work is produced by OpenStax-CNX and licensed under te Creative Commons Attribution License 3.0 Abstract Dene quantum number.
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationA SHORT INTRODUCTION TO BANACH LATTICES AND
CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,
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 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 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 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 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 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 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 informationMATH1131/1141 Calculus Test S1 v8a
MATH/ Calculus Test 8 S v8a October, 7 Tese solutions were written by Joann Blanco, typed by Brendan Trin and edited by Mattew Yan and Henderson Ko Please be etical wit tis resource It is for te use of
More informationUNIVERSITY OF MANITOBA DEPARTMENT OF MATHEMATICS MATH 1510 Applied Calculus I FIRST TERM EXAMINATION - Version A October 12, :30 am
DEPARTMENT OF MATHEMATICS MATH 1510 Applied Calculus I October 12, 2016 8:30 am LAST NAME: FIRST NAME: STUDENT NUMBER: SIGNATURE: (I understand tat ceating is a serious offense DO NOT WRITE IN THIS TABLE
More informationCopyright 2012 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future
Copyrigt 212 IEEE. Personal use of tis material is permitted. Permission from IEEE must be obtained for all oter uses, in any current or future media, including reprinting/republising tis material for
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 informationCalculus I Practice Exam 1A
Calculus I Practice Exam A Calculus I Practice Exam A Tis practice exam empasizes conceptual connections and understanding to a greater degree tan te exams tat are usually administered in introductory
More informationNumerical Solution to Parabolic PDE Using Implicit Finite Difference Approach
Numerical Solution to arabolic DE Using Implicit Finite Difference Approac Jon Amoa-Mensa, Francis Oene Boateng, Kwame Bonsu Department of Matematics and Statistics, Sunyani Tecnical University, Sunyani,
More information2016 PRELIM 2 PAPER 2 MARK SCHEME
06 River Valley Hig Scool Prelim Paper Mark Sceme 06 PRELIM PAPER MARK SCHEME (a) V 5.00 X 85. 9V 3 I.7 0 X V I X V I X 0.03 0. 85.9 5.00.7 X 48.3 00 X X 900 00 [A0] Anomalous data can be identified. Systematic
More information3.4 Algebraic Limits. Ex 1) lim. Ex 2)
Calculus Maimus.4 Algebraic Limits At tis point, you sould be very comfortable finding its bot grapically and numerically wit te elp of your graping calculator. Now it s time to practice finding its witout
More informationMonoidal Structures on Higher Categories
Monoidal Structures on Higer Categories Paul Ziegler Monoidal Structures on Simplicial Categories Let C be a simplicial category, tat is a category enriced over simplicial sets. Suc categories are a model
More informationChapter 2 Limits and Continuity
4 Section. Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 6) Quick Review.. f () ( ) () 4 0. f () 4( ) 4. f () sin sin 0 4. f (). 4 4 4 6. c c c 7. 8. c d d c d d c d c 9. 8 ( )(
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 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 information1 + t5 dt with respect to x. du = 2. dg du = f(u). du dx. dg dx = dg. du du. dg du. dx = 4x3. - page 1 -
Eercise. Find te derivative of g( 3 + t5 dt wit respect to. Solution: Te integrand is f(t + t 5. By FTC, f( + 5. Eercise. Find te derivative of e t2 dt wit respect to. Solution: Te integrand is f(t e t2.
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 informationControllability of a one-dimensional fractional heat equation: theoretical and numerical aspects
Controllability of a one-dimensional fractional eat equation: teoretical and numerical aspects Umberto Biccari, Víctor Hernández-Santamaría To cite tis version: Umberto Biccari, Víctor Hernández-Santamaría.
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 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 informationLecture 10: Carnot theorem
ecture 0: Carnot teorem Feb 7, 005 Equivalence of Kelvin and Clausius formulations ast time we learned tat te Second aw can be formulated in two ways. e Kelvin formulation: No process is possible wose
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 informationA First-Order System Approach for Diffusion Equation. I. Second-Order Residual-Distribution Schemes
A First-Order System Approac for Diffusion Equation. I. Second-Order Residual-Distribution Scemes Hiroaki Nisikawa W. M. Keck Foundation Laboratory for Computational Fluid Dynamics, Department of Aerospace
More information