1 The Euler Forward scheme (schéma d Euler explicite)
|
|
- Anissa Gardner
- 6 years ago
- Views:
Transcription
1 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, Te Euler Forward sceme (scéma d Euler explicite) We look for a umerical approximatio of te Europea put fuctio v = v(t, s), t [0, T ], s [0, S max ] It satisfies i first approximatio te Black ad Scoles backward PDE o te trucated domai Ω = [S mi, S max ]: s2 2 s 2 v rs s v + rv = 0, t (0, T ), s (S mi, S max ) t v σ2 2 v(t, S mi ) = v l (t) Ke rt S mi, t (0, T ) v(t, S max ) = v r (t) 0, t (0, T ) v(0, s) = ϕ(s) := (K s) +, s (S mi, S max ) We will cosider te followig parameters: K = 100, S mi = 0, S max = 200, T = 1, σ = 02, r = 01 I particular, we aim at computig v(t, s) at fial time t = T We first itroduce a discrete mes as follows Let := Smax S mi I+1 be spatial mes step, ad τ := T N te time step Te s := S mi +, = 0,, I + 1 are te mes poits, ad t = τ, = 0,, N te time mes We are lookig for U, a approximatio of v(t, s ) For ay fuctio v C 2 (or v C 3 for (4)), we recall te followig approximatios, as 0, v (s ) = v(s ) v(s 1 ) v (s ) = v(s +1) v(s ) v (s ) = v(s +1) v(s 1 ) (1) + O() (2) + O() (3) + O( 2 ) (4) We terefore obtai several possible approximatios by fiite differeces for te first order derivative: s v(t, s ) U U 1 (backward differece approximatio) (5) s v(t, s ) U +1 U (forward differece approximatio) (6) s v(t, s ) U +1 U 1 (cetered approximatio) (7) Te first two approximatios are said to be cosistet of order 1 (i space), wile te secod oe is cosistet of order 2 1
2 We also recall te approximatio wic is secod-order cosistet i space 2 s 2 v(t, s ) U 1 + 2U U +1 2, (8) Hece we obtai te so-called "Euler Forward sceme" (or Explicit Euler sceme), abreviated "EE" usig te cetered approximatio, as follows: U +1 τ U + σ2 2 s2 U 1 + 2U U +1 2 rs U +1 U 1 + ru = 0 = 0,, N 1, = 1,, I U 0 = v l (t ) Ke rt S mi, = 0,, N (9) U I+1 = v r (t ) 0, U 0 = ϕ(s ) (K s ) +, = 0,, N = 1,, I Let us remark tat we ave take = 1 ad = I as extremal idices i For = 1, te sceme utilizes te kow value U0 := v l(t ) (left boudary value) For = I, te sceme utilizes te kow value UI+1 := v r(t ) (rigt boudary value) 2 Programmig Euler Forward 21 Prelimiaries We coose to work wit te ukow te vector correspodig to (v(t, s )) =1,,I : U = U 1 U I We would like to write (9) uder te vector form as follows: U +1 U + AU + q(t ) = 0, (10) τ were A is a square matrix of dimesio I ad q(t) is a colum vector of size I Let us deote α := σ2 s 2 2 2, β := r s We look for A ad q(t) suc tat α i ( U i 1 + 2U i U i+1) β i (U i+1 U i 1) + ru i = ( α i + β i )Ui 1 + (2α i + r)ui + ( α i β i )Ui+1 (AU + q(t )) i 2
3 By idetificatio we see tat A is a tridiagoal matrix 2α 1 + r α 1 β 1 0 α 2 + β 2 2α 2 + r α 2 β 2 A := α i + β i 2α i + r α i β i 0 α I + β I 2α I + r ad q(t ) cotais te boudary values U0 = v l(t ) ad UI+1 = v r(t ) as follows: ( α 1 + β 1 )U0 ( α 1 + β 1 )v l (t ) 0 0 q(t ) := 0 ( α I β I )U I+1 0 ( α I β I )v r (t ) 22 Gettig ito te program a) Dowload te workig program 1 (tp1m if usig Matlab, or tp1sci if usig Scilab) Tis program as some lies to be completed i order to work properly I te case of te matlab program, you sould also dowload te oter files (BSm, plootm, orderm ad te solutio file tp1solm) b) Ceck te program for te payoff fuctio u0 (for ϕ), ul (for v l ) ad ur (for v r ) ad program tem correctly (otice tat te istructio y=max(k-s,0) works for s vector, K scalar, ad retur a vector of same size as s) c) I te sectio MESH, complete te value of te mes step ad program te vector s cotaiig te (s ) 1 I values Typically it sould look like: dt=t/n; =(Smax-Smi)/(I+1); s=smi+ (1:I) *; // colum vector (s 1,, s I ) T d) Program te matrix A ad test te program Uder te matlab commad widow, type tp1 (Scilab: type exec tp1sci or exec( tp1sci,-1)) e) Program te fuctio q(t) For istace i matlab: q [(-alpa(1) + bet(1))* ul(t); zeros(i-2,1); (-alpa(ed) - bet(ed))* ur(t)]; I Scilab tis may look like fuctio y=ul(t); y=k*exp(-r*t); edfuctio fuctio y=ur(t); y=0; edfuctio; fuctio 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); edfuctio 1 See ttps://lllmatupmcfr/bokaowski/eseigemet/2017/m2mo/m2motml 3
4 22 Euler Forward sceme (or "Euler Explicit" sceme) a) Program te explicit form of te vector U +1, i terms of U, i te mai loop, usig te matrix A ad te fuctio q I te ed, it sould look like (i matlab) case EE P = (Id - dt*a)*p - dt*q(t); Note tat te grapic fuctio ploot also plots te exact Black ad Scoles formula (see te give fuctio BS i te file BSm) b) Correct te lig errli=0 i te mai loop, i order to compute correctly te maximum orm betwee te sceme values ad te Black ad Scoles values: U V BS = max 1 i I U i V BS (s i ) (use orm(x, if ) to get te l orm of a give vector x) 23 Solutio A solutio file is give i tp1solm 3 First umerical tests a) Test te Euler forward sceme (EE) First fix N = 10 ad take I = 10, 20, 50, Te take te followig N = I values : 10, 20, 50, 100 Observe tat: - te sceme is ot always umerically stable - it does ot always give a positive solutio (ie we do ot always ave U 0, ) b) Uderstad te origi of te oscillatios we tey occur For istace, fix N = 10 ad I = 50, ad look at te "amplificatio" matrix defied as B := I d τa (ceck tat te coefficiets of B are ot all positive ad tat tey may ave a modulus greater ta 1) c) Fill i te CFL umber defied ere as µ := τ 2 σ2 S 2 max ad prit it Ceck tat tere is o stability problem we µ is sufficietly small d) I te case we σ = 0002 ad for istace wit N = 10 ad I = 50: observe tat te sceme is stable but tat te solutio give by te sceme is ot positive e) Order of te sceme We first cosider te followig values of I ad of N: I = 10, 20, 40, 80, 160, ad N = I 2 /10 (tat is, N = 10, 40, ) Tis is i order to ave τ 2 Fill i te followig error tables ad compute te correspodig (spatial) order of te metod (you may use te file orderm) More precisely, if te error is e k for a give 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 ) were k is te spatial mes step at step k Te idea is to try to detect a beavior of te form e k = C α k, were C est a costat ad were k is te spatial mes size correspodig to I k 4
5 I N = I e k order α k I N I 2 e k order α k So if k 0 /2 k, as i te left table, te we get e k /e k 1 1/2 α ad te previous formula gives α k = log(e k 1/e k ) log(2) wic sould give a estimate of α Te umerical (spatial) order sould be close to two But ere settig N cost I 2 is costly i terms of umber of operatios More precisely, sice τ = cost 2 oe ca see tat error also beaves as O(τ), tat is, a first order error beavior wit respect to te time approximatio Tis motivates te use of implicit scemes i order to avoid te time-step coditio (or "CFL" coditio) Exercice: backward ad forward differece approximatios a) Program i te same way te oter approximatios (backward differeces / décetrage droit : DROIT ad te forward differeces - décetrage gauce GAUCHE ) Te parameter CENTRAGE is defied at te begiig of te program (sectio NUMERICAL DATA) ad determies te type of fiite differece approximatio tat is used b) Usig te forward differeces, fix temporarily σ = 0002 ad N = 10, I = 50, ad ceck tat ow te solutio keeps positive c) Usig te backward differeces, fix temporarily σ = 0002 ad N = 10, I = 50, ad ceck tat ow te umerical solutio is ustable (uderstad te problem by lookig at te coefficiets of te amplificatio matrix) 4 Implicit Euler sceme (scéma d Euler implicite) Te CFL costrait imposes some restrictio o te time step τ Implicit scemes may allow us to get rid of tis restrictio For istace, te implicit euler sceme, ereafter abrieviated "EI" (for Frec "Euler Implicit"), wit cetered differece approximatio for te first spatial derivatives reads U +1 τ U + σ2 2 s2 U U U+1 +1 U rs U ru +1 = 0 = 0,, N 1, = 1,, I U0 +1 = v l (t +1 ) Ke rt +1 S mi, = 0,, N 1 (11) U +1 I+1 = v r(t +1 ) 0, = 0,, N 1 U 0 = (K s ), = 1,, I 5
6 a) Sow tat te sceme, i vector form, i aalogy wit (10), is (U +1 U )/τ + AU +1 + q(t +1 ) = 0 b) Program EI: set te parameter SCHEMA= EI at te begiig of te program, ad complete accordigly, i te mai loop, te part case EI Note tat i order to solve a liear system of te form Ax = b oe may use te liear solver i matlab (resp Scilab) as follows x=a\b; c) Ceck tat wit te EI sceme tere is o more stability problems (wit for istace N = 10 ad I = 50) 5 Crak-Nicolso sceme a) Program te Crak-Nicolso sceme (CN) (θ-sceme wit θ = 1 2 ) I vector form te sceme reads (U +1 U )/τ (AU +1 + q(t +1 ) (AU + q(t )) = 0 b) Numerically study te covergece of te cetered CN sceme To do so, oe ca compute te umerical error obtaied wit te followig values N = I + 1 {10, 20, 40, 80, 160} Te do te same study as for te Euler implicit sceme ad compare by makig a error table (ad wit a order estimatio) Oe sould observe a secod order beavior for CN, wile it sould be first order for EI c) Ceck tat tese results are coeret wit te teory 6 Exercices Exercice 1 ("Call optio") a) Usig te put-call parity formula, do a fuctio tat also computes te black ad scoles formula for te call b) Propose adequate left ad rigt boudary coditios (tat is, at s = S mi ad s = S max ) for te call optio, i te form (v(t, S mi ) = v l (t) ad v(t, S max ) = v r (t), were v l, v r are fuctio to determie 2 c) Propose a apropriate PDE for te call optio, usig tese boudary coditios d) Write a Euler explicit sceme for te call optio, ad program it Exercice 2 (Improved efficiecy usig sparse matrices) Te matrix A as oly a few o zero elemets (about 3I o zero elemets) It is possible to code oly te o-zero elemets of A by usig te matrix type sparse Type elp sparse or doc sparse for documetatio i matlab A=spzeros(I,I) full(a) sparse(a) speye(i,i), speye(i) iitialize a sparse 0-matrix of size I I create correspodig full matrix of a sparse matrix create te correspodig sparse matrix of A create a sparse idetity matrix of size I Program te EI or CN sceme by usig oly sparse matrices, ad compare te speed of te ew code wit respect to te old oe for large I values Executio time is i geeral reduced It ca be evaluated by usig te commads tic; toc; (resp time() ad time() i Scilab) as follows: 2 v l (t, s) 0, v r(t, s) s Ke rt 6
7 tic; % INSTRUCTIONS; % t=toc; pritf( cpu time : t=%52f, t); (te te variable t cotais te time elapsed betwee te calls of tic ad toc) 7
1 The Euler Forward scheme (schéma d'euler explicite)
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 2016 1 Te Euler Forward sceme (scéma d'euler explicite)
More informationx x x 2x x N ( ) p NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS By Newton-Raphson formula
NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS. If g( is cotiuous i [a,b], te uder wat coditio te iterative (or iteratio metod = g( as a uique solutio i [a,b]? '( i [a,b].. Wat
More informationThe Advection-Diffusion equation!
ttp://www.d.edu/~gtryggva/cf-course/! Te Advectio-iffusio equatio! Grétar Tryggvaso! Sprig 3! Navier-Stokes equatios! Summary! u t + u u x + v u y = P ρ x + µ u + u ρ y Hyperbolic part! u x + v y = Elliptic
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationStability analysis of numerical methods for stochastic systems with additive noise
Stability aalysis of umerical metods for stoctic systems wit additive oise Yosiiro SAITO Abstract Stoctic differetial equatios (SDEs) represet pysical peomea domiated by stoctic processes As for determiistic
More informationME 501A Seminar in Engineering Analysis Page 1
Accurac, Stabilit ad Sstems of Equatios November 0, 07 Numerical Solutios of Ordiar Differetial Equatios Accurac, Stabilit ad Sstems of Equatios Larr Caretto Mecaical Egieerig 0AB Semiar i Egieerig Aalsis
More informationPartial Differential Equations
EE 84 Matematical Metods i Egieerig Partial Differetial Eqatios Followig are some classical partial differetial eqatios were is assmed to be a fctio of two or more variables t (time) ad y (spatial coordiates).
More informationMath 257: Finite difference methods
Math 257: Fiite differece methods 1 Fiite Differeces Remember the defiitio of a derivative f f(x + ) f(x) (x) = lim 0 Also recall Taylor s formula: (1) f(x + ) = f(x) + f (x) + 2 f (x) + 3 f (3) (x) +...
More informationA widely used display of protein shapes is based on the coordinates of the alpha carbons - - C α
Nice plottig of proteis: I A widely used display of protei shapes is based o the coordiates of the alpha carbos - - C α -s. The coordiates of the C α -s are coected by a cotiuous curve that roughly follows
More informationLECTURE 2 LEAST SQUARES CROSS-VALIDATION FOR KERNEL DENSITY ESTIMATION
Jauary 3 07 LECTURE LEAST SQUARES CROSS-VALIDATION FOR ERNEL DENSITY ESTIMATION Noparametric kerel estimatio is extremely sesitive to te coice of badwidt as larger values of result i averagig over more
More informationReview Problems 1. ICME and MS&E Refresher Course September 19, 2011 B = C = AB = A = A 2 = A 3... C 2 = C 3 = =
Review Problems ICME ad MS&E Refresher Course September 9, 0 Warm-up problems. For the followig matrices A = 0 B = C = AB = 0 fid all powers A,A 3,(which is A times A),... ad B,B 3,... ad C,C 3,... Solutio:
More informationA sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece,, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet as
More informationDifferentiation Techniques 1: Power, Constant Multiple, Sum and Difference Rules
Differetiatio Teciques : Power, Costat Multiple, Sum ad Differece Rules 97 Differetiatio Teciques : Power, Costat Multiple, Sum ad Differece Rules Model : Fidig te Equatio of f '() from a Grap of f ()
More informationSome examples of vector spaces
Roberto s Notes o Liear Algebra Chapter 11: Vector spaces Sectio 2 Some examples of vector spaces What you eed to kow already: The te axioms eeded to idetify a vector space. What you ca lear here: Some
More informationInfinite Sequences and Series
Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet
More informationMore Elementary Aspects of Numerical Solutions of PDEs!
ttp://www.d.edu/~gtryggva/cfd-course/ Outlie More Elemetary Aspects o Numerical Solutios o PDEs I tis lecture we cotiue to examie te elemetary aspects o umerical solutios o partial dieretial equatios.
More informationFourier Techniques lecture by Håkon Hoel
CSC Nada DN55 Sprig - Differetial Equatios II JOp p (7) Fourier Teciques lecture by Håo Hoel Fourier series... Applicatios to liear PDE... 3 Numerical metods...3 3. Vo Neuma aalysis...4 3. Spectral metods...6
More information( ) (( ) ) ANSWERS TO EXERCISES IN APPENDIX B. Section B.1 VECTORS AND SETS. Exercise B.1-1: Convex sets. are convex, , hence. and. (a) Let.
Joh Riley 8 Jue 03 ANSWERS TO EXERCISES IN APPENDIX B Sectio B VECTORS AND SETS Exercise B-: Covex sets (a) Let 0 x, x X, X, hece 0 x, x X ad 0 x, x X Sice X ad X are covex, x X ad x X The x X X, which
More informationMost text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t
Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said
More informationECE 901 Lecture 12: Complexity Regularization and the Squared Loss
ECE 90 Lecture : Complexity Regularizatio ad the Squared Loss R. Nowak 5/7/009 I the previous lectures we made use of the Cheroff/Hoeffdig bouds for our aalysis of classifier errors. Hoeffdig s iequality
More informationFinite Difference Derivations for Spreadsheet Modeling John C. Walton Modified: November 15, 2007 jcw
Fiite Differece Derivatios for Spreadsheet Modelig Joh C. Walto Modified: November 15, 2007 jcw Figure 1. Suset with 11 swas o Little Platte Lake, Michiga. Page 1 Modificatio Date: November 15, 2007 Review
More informationTMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods
TMA4205 Numerical Liear Algebra The Poisso problem i R 2 : diagoalizatio methods September 3, 2007 c Eiar M Røquist Departmet of Mathematical Scieces NTNU, N-749 Trodheim, Norway All rights reserved A
More informationThe Random Walk For Dummies
The Radom Walk For Dummies Richard A Mote Abstract We look at the priciples goverig the oe-dimesioal discrete radom walk First we review five basic cocepts of probability theory The we cosider the Beroulli
More informationOn Exact Finite-Difference Scheme for Numerical Solution of Initial Value Problems in Ordinary Differential Equations.
O Exact Fiite-Differece Sceme for Numerical Solutio of Iitial Value Problems i Ordiar Differetial Equatios. Josua Suda, M.Sc. Departmet of Matematical Scieces, Adamawa State Uiversit, Mubi, Nigeria. E-mail:
More informationProblem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t =
Mathematics Summer Wilso Fial Exam August 8, ANSWERS Problem 1 (a) Fid the solutio to y +x y = e x x that satisfies y() = 5 : This is already i the form we used for a first order liear differetial equatio,
More informationMa 530 Introduction to Power Series
Ma 530 Itroductio to Power Series Please ote that there is material o power series at Visual Calculus. Some of this material was used as part of the presetatio of the topics that follow. What is a Power
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationPrinciple Of Superposition
ecture 5: PREIMINRY CONCEP O RUCUR NYI Priciple Of uperpositio Mathematically, the priciple of superpositio is stated as ( a ) G( a ) G( ) G a a or for a liear structural system, the respose at a give
More informationLecture 7 Testing Nonlinear Inequality Restrictions 1
Eco 75 Lecture 7 Testig Noliear Iequality Restrictios I Lecture 6, we discussed te testig problems were te ull ypotesis is de ed by oliear equality restrictios: H : ( ) = versus H : ( ) 6= : () We sowed
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More informationThe Method of Least Squares. To understand least squares fitting of data.
The Method of Least Squares KEY WORDS Curve fittig, least square GOAL To uderstad least squares fittig of data To uderstad the least squares solutio of icosistet systems of liear equatios 1 Motivatio Curve
More informationFinite Difference Method for the Estimation of a Heat Source Dependent on Time Variable ABSTRACT
Malaysia Joural of Matematical Scieces 6(S): 39-5 () Special Editio of Iteratioal Worsop o Matematical Aalysis (IWOMA) Fiite Differece Metod for te Estimatio of a Heat Source Depedet o Time Variable, Allabere
More informationSlide Set 13 Linear Model with Endogenous Regressors and the GMM estimator
Slide Set 13 Liear Model with Edogeous Regressors ad the GMM estimator Pietro Coretto pcoretto@uisa.it Ecoometrics Master i Ecoomics ad Fiace (MEF) Uiversità degli Studi di Napoli Federico II Versio: Friday
More informationSequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet
More informationb i u x i U a i j u x i u x j
M ath 5 2 7 Fall 2 0 0 9 L ecture 1 9 N ov. 1 6, 2 0 0 9 ) S ecod- Order Elliptic Equatios: Weak S olutios 1. Defiitios. I this ad the followig two lectures we will study the boudary value problem Here
More informationSOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.
SOLUTIONS TO EXAM 3 Problem Fid the sum of the followig series 2 + ( ) 5 5 2 5 3 25 2 2 This series diverges Solutio: Note that this defies two coverget geometric series with respective radii r 2/5 < ad
More informationFurther Methods for Advanced Mathematics (FP2) WEDNESDAY 9 JANUARY 2008
ADVANCED GCE 7/ MATHEMATICS (MEI) Furter Metods for Advaced Matematics (F) WEDNESDAY 9 JANUARY 8 Additioal materials: Aswer Booklet (8 pages) Grap paper MEI Eamiatio Formulae ad Tables (MF) Afteroo Time:
More informationRecurrence Relations
Recurrece Relatios Aalysis of recursive algorithms, such as: it factorial (it ) { if (==0) retur ; else retur ( * factorial(-)); } Let t be the umber of multiplicatios eeded to calculate factorial(). The
More informationInteger Linear Programming
Iteger Liear Programmig Itroductio Iteger L P problem (P) Mi = s. t. a = b i =,, m = i i 0, iteger =,, c Eemple Mi z = 5 s. t. + 0 0, 0, iteger F(P) = feasible domai of P Itroductio Iteger L P problem
More informationMATH 31B: MIDTERM 2 REVIEW
MATH 3B: MIDTERM REVIEW JOE HUGHES. Evaluate x (x ) (x 3).. Partial Fractios Solutio: The umerator has degree less tha the deomiator, so we ca use partial fractios. Write x (x ) (x 3) = A x + A (x ) +
More informationd y f f dy Numerical Solution of Ordinary Differential Equations Consider the 1 st order ordinary differential equation (ODE) . dx
umerical Solutio o Ordiar Dieretial Equatios Cosider te st order ordiar dieretial equatio ODE d. d Te iitial coditio ca be tae as. Te we could use a Talor series about ad obtai te complete solutio or......!!!
More informationLIMITS AND DERIVATIVES
Capter LIMITS AND DERIVATIVES. Overview.. Limits of a fuctio Let f be a fuctio defied i a domai wic we take to be a iterval, say, I. We sall study te cocept of it of f at a poit a i I. We say f ( ) is
More information9.3 Power Series: Taylor & Maclaurin Series
9.3 Power Series: Taylor & Maclauri Series If is a variable, the a ifiite series of the form 0 is called a power series (cetered at 0 ). a a a a a 0 1 0 is a power series cetered at a c a a c a c a c 0
More informationLecture 8: October 20, Applications of SVD: least squares approximation
Mathematical Toolkit Autum 2016 Lecturer: Madhur Tulsiai Lecture 8: October 20, 2016 1 Applicatios of SVD: least squares approximatio We discuss aother applicatio of sigular value decompositio (SVD) of
More informationMath 116 Practice for Exam 3
Math 6 Practice for Eam 3 Geerated April 4, 26 Name: SOLUTIONS Istructor: Sectio Number:. This eam has questios. Note that the problems are ot of equal difficulty, so you may wat to skip over ad retur
More informationNEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE
UPB Sci Bull, Series A, Vol 79, Iss, 207 ISSN 22-7027 NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE Gabriel Bercu We itroduce two ew sequeces of Euler-Mascheroi type which have fast covergece
More informationLecture 2: Finite Difference Methods in Heat Transfer
Lecture 2: Fiite Differece Methods i Heat Trasfer V.Vuorie Aalto Uiversity School of Egieerig Heat ad Mass Trasfer Course, Autum 2016 November 1 st 2017, Otaiemi ville.vuorie@aalto.fi Overview Part 1 (
More informationChimica Inorganica 3
himica Iorgaica Irreducible Represetatios ad haracter Tables Rather tha usig geometrical operatios, it is ofte much more coveiet to employ a ew set of group elemets which are matrices ad to make the rule
More informationEssential Microeconomics EXISTENCE OF EQUILIBRIUM Core ideas: continuity of excess demand functions, Fixed point theorems
Essetial Microecoomics -- 5.3 EXISTENCE OF EQUILIBRIUM Core ideas: cotiuity of excess demad fuctios, Fixed oit teorems Two commodity excage ecoomy 2 Excage ecoomy wit may commodities 5 Discotiuous demad
More informationCHAPTER 5. Theory and Solution Using Matrix Techniques
A SERIES OF CLASS NOTES FOR 2005-2006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 3 A COLLECTION OF HANDOUTS ON SYSTEMS OF ORDINARY DIFFERENTIAL
More informationLecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods
Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace
More informationDEGENERACY AND ALL THAT
DEGENERACY AND ALL THAT Te Nature of Termodyamics, Statistical Mecaics ad Classical Mecaics Termodyamics Te study of te equilibrium bulk properties of matter witi te cotext of four laws or facts of experiece
More informationLIMITS AND DERIVATIVES NCERT
. Overview.. Limits of a fuctio Let f be a fuctio defied i a domai wic we take to be a iterval, say, I. We sall study te cocept of it of f at a poit a i I. We say f ( ) is te epected value of f at a give
More informationChapter 4 : Laplace Transform
4. Itroductio Laplace trasform is a alterative to solve the differetial equatio by the complex frequecy domai ( s = σ + jω), istead of the usual time domai. The DE ca be easily trasformed ito a algebraic
More informationThe multiplicative structure of finite field and a construction of LRC
IERG6120 Codig for Distributed Storage Systems Lecture 8-06/10/2016 The multiplicative structure of fiite field ad a costructio of LRC Lecturer: Keeth Shum Scribe: Zhouyi Hu Notatios: We use the otatio
More informationEfficient GMM LECTURE 12 GMM II
DECEMBER 1 010 LECTURE 1 II Efficiet The estimator depeds o the choice of the weight matrix A. The efficiet estimator is the oe that has the smallest asymptotic variace amog all estimators defied by differet
More informationIt is often useful to approximate complicated functions using simpler ones. We consider the task of approximating a function by a polynomial.
Taylor Polyomials ad Taylor Series It is ofte useful to approximate complicated fuctios usig simpler oes We cosider the task of approximatig a fuctio by a polyomial If f is at least -times differetiable
More informationIterative Techniques for Solving Ax b -(3.8). Assume that the system has a unique solution. Let x be the solution. Then x A 1 b.
Iterative Techiques for Solvig Ax b -(8) Cosider solvig liear systems of them form: Ax b where A a ij, x x i, b b i Assume that the system has a uique solutio Let x be the solutio The x A b Jacobi ad Gauss-Seidel
More informationALLOCATING SAMPLE TO STRATA PROPORTIONAL TO AGGREGATE MEASURE OF SIZE WITH BOTH UPPER AND LOWER BOUNDS ON THE NUMBER OF UNITS IN EACH STRATUM
ALLOCATING SAPLE TO STRATA PROPORTIONAL TO AGGREGATE EASURE OF SIZE WIT BOT UPPER AND LOWER BOUNDS ON TE NUBER OF UNITS IN EAC STRATU Lawrece R. Erst ad Cristoper J. Guciardo Erst_L@bls.gov, Guciardo_C@bls.gov
More information, then cv V. Differential Equations Elements of Lineaer Algebra Name: Consider the differential equation. and y2 cos( kx)
Cosider the differetial equatio y '' k y 0 has particular solutios y1 si( kx) ad y cos( kx) I geeral, ay liear combiatio of y1 ad y, cy 1 1 cy where c1, c is also a solutio to the equatio above The reaso
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationThe z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j
The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.
More informationREGRESSION (Physics 1210 Notes, Partial Modified Appendix A)
REGRESSION (Physics 0 Notes, Partial Modified Appedix A) HOW TO PERFORM A LINEAR REGRESSION Cosider the followig data poits ad their graph (Table I ad Figure ): X Y 0 3 5 3 7 4 9 5 Table : Example Data
More informationSummary: CORRELATION & LINEAR REGRESSION. GC. Students are advised to refer to lecture notes for the GC operations to obtain scatter diagram.
Key Cocepts: 1) Sketchig of scatter diagram The scatter diagram of bivariate (i.e. cotaiig two variables) data ca be easily obtaied usig GC. Studets are advised to refer to lecture otes for the GC operatios
More informationOn the convergence, consistence and stability of a standard finite difference scheme
AMERICAN JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH 2, Sciece Huβ, ttp://www.sciub.org/ajsir ISSN: 253-649X, doi:.525/ajsir.2.2.2.74.78 O te covergece, cosistece ad stabilit of a stadard fiite differece
More information5.1 Review of Singular Value Decomposition (SVD)
MGMT 69000: Topics i High-dimesioal Data Aalysis Falll 06 Lecture 5: Spectral Clusterig: Overview (cotd) ad Aalysis Lecturer: Jiamig Xu Scribe: Adarsh Barik, Taotao He, September 3, 06 Outlie Review of
More informationM A T H F A L L CORRECTION. Algebra I 1 4 / 1 0 / U N I V E R S I T Y O F T O R O N T O
M A T H 2 4 0 F A L L 2 0 1 4 HOMEWORK ASSIGNMENT #4 CORRECTION Algebra I 1 4 / 1 0 / 2 0 1 4 U N I V E R S I T Y O F T O R O N T O P r o f e s s o r : D r o r B a r - N a t a Correctio Homework Assigmet
More informationOptimal Estimator for a Sample Set with Response Error. Ed Stanek
Optial Estiator for a Saple Set wit Respose Error Ed Staek Itroductio We develop a optial estiator siilar to te FP estiator wit respose error tat was cosidered i c08ed63doc Te first 6 pages of tis docuet
More informationTHE ENERGY BALANCE ERROR FOR CIRCUIT TRANSIENT ANALYSIS
THE ENERGY BALANCE ERROR FOR CIRCUIT TRANSIENT ANALYSIS FLORIN CONSTANTINESCU, ALEXANDRU GABRIEL GHEORGHE, MIRUNA NIŢESCU Key words: Trasiet aalysis, Eergy balace error, Time step coice. Two algoritms
More informationLøsningsførslag i 4M
Norges tekisk aturviteskapelige uiversitet Istitutt for matematiske fag Side 1 av 6 Løsigsførslag i 4M Oppgave 1 a) A sketch of the graph of the give f o the iterval [ 3, 3) is as follows: The Fourier
More informationSection 13.3 Area and the Definite Integral
Sectio 3.3 Area ad the Defiite Itegral We ca easily fid areas of certai geometric figures usig well-kow formulas: However, it is t easy to fid the area of a regio with curved sides: METHOD: To evaluate
More informationLet us give one more example of MLE. Example 3. The uniform distribution U[0, θ] on the interval [0, θ] has p.d.f.
Lecture 5 Let us give oe more example of MLE. Example 3. The uiform distributio U[0, ] o the iterval [0, ] has p.d.f. { 1 f(x =, 0 x, 0, otherwise The likelihood fuctio ϕ( = f(x i = 1 I(X 1,..., X [0,
More informationSection 11.8: Power Series
Sectio 11.8: Power Series 1. Power Series I this sectio, we cosider geeralizig the cocept of a series. Recall that a series is a ifiite sum of umbers a. We ca talk about whether or ot it coverges ad i
More informationProblem Set 4 Due Oct, 12
EE226: Radom Processes i Systems Lecturer: Jea C. Walrad Problem Set 4 Due Oct, 12 Fall 06 GSI: Assae Gueye This problem set essetially reviews detectio theory ad hypothesis testig ad some basic otios
More informationYou may work in pairs or purely individually for this assignment.
CS 04 Problem Solvig i Computer Sciece OOC Assigmet 6: Recurreces You may work i pairs or purely idividually for this assigmet. Prepare your aswers to the followig questios i a plai ASCII text file or
More informationEXPLICIT FORMULAE FOR THE BRENTI S POLYNOMIALS γ a1,...,a r
EXPLICIT FORMULAE FOR THE BRENTI S POLYNOMIALS γ,,a r Jia-yi Si Departmet of Matematics, East Cia Normal Uiversity, Sagai, 0006, P R Cia Astract We itroduce two kids of geeralized Catala umers a a, Z ad
More informationSummary and Discussion on Simultaneous Analysis of Lasso and Dantzig Selector
Summary ad Discussio o Simultaeous Aalysis of Lasso ad Datzig Selector STAT732, Sprig 28 Duzhe Wag May 4, 28 Abstract This is a discussio o the work i Bickel, Ritov ad Tsybakov (29). We begi with a short
More informationw (1) ˆx w (1) x (1) /ρ and w (2) ˆx w (2) x (2) /ρ.
2 5. Weighted umber of late jobs 5.1. Release dates ad due dates: maximimizig the weight of o-time jobs Oce we add release dates, miimizig the umber of late jobs becomes a sigificatly harder problem. For
More informationCALCULUS BASIC SUMMER REVIEW
CALCULUS BASIC SUMMER REVIEW NAME rise y y y Slope of a o vertical lie: m ru Poit Slope Equatio: y y m( ) The slope is m ad a poit o your lie is, ). ( y Slope-Itercept Equatio: y m b slope= m y-itercept=
More informationAlgebra of Least Squares
October 19, 2018 Algebra of Least Squares Geometry of Least Squares Recall that out data is like a table [Y X] where Y collects observatios o the depedet variable Y ad X collects observatios o the k-dimesioal
More informationImprovement of Generic Attacks on the Rank Syndrome Decoding Problem
Improvemet of Geeric Attacks o the Rak Sydrome Decodig Problem Nicolas Arago, Philippe Gaborit, Adrie Hauteville, Jea-Pierre Tillich To cite this versio: Nicolas Arago, Philippe Gaborit, Adrie Hauteville,
More information(b) What is the probability that a particle reaches the upper boundary n before the lower boundary m?
MATH 529 The Boudary Problem The drukard s walk (or boudary problem) is oe of the most famous problems i the theory of radom walks. Oe versio of the problem is described as follows: Suppose a particle
More informationAn Insight into Differentiation and Integration
Differetiatio A Isigt ito Differetiatio a Itegratio Differetiatio is basically a task to fi out ow oe variable is cagig i relatio to aoter variable, te latter is usually take as a cause of te cage. For
More informationZeros of Polynomials
Math 160 www.timetodare.com 4.5 4.6 Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered with fidig the solutios of polyomial equatios of ay degree
More informationLinear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy
Liear Differetial Equatios of Higher Order Basic Theory: Iitial-Value Problems d y d y dy Solve: a( ) + a ( )... a ( ) a0( ) y g( ) + + + = d d d ( ) Subject to: y( 0) = y0, y ( 0) = y,..., y ( 0) = y
More informationLinear Elliptic PDE s Elliptic partial differential equations frequently arise out of conservation statements of the form
Liear Elliptic PDE s Elliptic partial differetial equatios frequetly arise out of coservatio statemets of the form B F d B Sdx B cotaied i bouded ope set U R. Here F, S deote respectively, the flux desity
More informationINFINITE SEQUENCES AND SERIES
11 INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES 11.4 The Compariso Tests I this sectio, we will lear: How to fid the value of a series by comparig it with a kow series. COMPARISON TESTS
More informationCS 270 Algorithms. Oliver Kullmann. Growth of Functions. Divide-and- Conquer Min-Max- Problem. Tutorial. Reading from CLRS for week 2
Geeral remarks Week 2 1 Divide ad First we cosider a importat tool for the aalysis of algorithms: Big-Oh. The we itroduce a importat algorithmic paradigm:. We coclude by presetig ad aalysig two examples.
More informationCHAPTER 5 SOME MINIMAX AND SADDLE POINT THEOREMS
CHAPTR 5 SOM MINIMA AND SADDL POINT THORMS 5. INTRODUCTION Fied poit theorems provide importat tools i game theory which are used to prove the equilibrium ad eistece theorems. For istace, the fied poit
More informationMachine Learning Brett Bernstein
Machie Learig Brett Berstei Week 2 Lecture: Cocept Check Exercises Starred problems are optioal. Excess Risk Decompositio 1. Let X = Y = {1, 2,..., 10}, A = {1,..., 10, 11} ad suppose the data distributio
More informationNICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =
AN INTRODUCTION TO SCHRÖDER AND UNKNOWN NUMBERS NICK DUFRESNE Abstract. I this article we will itroduce two types of lattice paths, Schröder paths ad Ukow paths. We will examie differet properties of each,
More informationΩ ). Then the following inequality takes place:
Lecture 8 Lemma 5. Let f : R R be a cotiuously differetiable covex fuctio. Choose a costat δ > ad cosider the subset Ωδ = { R f δ } R. Let Ωδ ad assume that f < δ, i.e., is ot o the boudary of f = δ, i.e.,
More informationDefinition 4.2. (a) A sequence {x n } in a Banach space X is a basis for X if. unique scalars a n (x) such that x = n. a n (x) x n. (4.
4. BASES I BAACH SPACES 39 4. BASES I BAACH SPACES Sice a Baach space X is a vector space, it must possess a Hamel, or vector space, basis, i.e., a subset {x γ } γ Γ whose fiite liear spa is all of X ad
More informationThe Jordan Normal Form: A General Approach to Solving Homogeneous Linear Systems. Mike Raugh. March 20, 2005
The Jorda Normal Form: A Geeral Approach to Solvig Homogeeous Liear Sstems Mike Raugh March 2, 25 What are we doig here? I this ote, we describe the Jorda ormal form of a matrix ad show how it ma be used
More informationA Pseudo Spline Methods for Solving an Initial Value Problem of Ordinary Differential Equation
Joural of Matematics ad Statistics 4 (: 7-, 008 ISSN 549-3644 008 Sciece Publicatios A Pseudo Splie Metods for Solvig a Iitial Value Problem of Ordiary Differetial Equatio B.S. Ogudare ad G.E. Okeca Departmet
More informationFIR Filter Design: Part II
EEL335: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we cosider how we might go about desigig FIR filters with arbitrary frequecy resposes, through compositio of multiple sigle-peak
More information1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations
. Liearizatio of a oliear system give i the form of a system of ordiary differetial equatios We ow show how to determie a liear model which approximates the behavior of a time-ivariat oliear system i a
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationA TYPE OF PRIMITIVE ALGEBRA*
A TYPE OF PRIMITIVE ALGEBRA* BT J. H. M. WEDDERBURN I a recet paper,t L. E. Dickso has discussed the liear associative algebra, A, defied by the relatios xy = yo(x), y = g, where 8 ( x ) is a polyomial
More information1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y
Questio (a) A square matrix A= A is called positive defiite if the quadratic form waw > 0 for every o-zero vector w [Note: Here (.) deotes the traspose of a matrix or a vector]. Let 0 A = 0 = show that:
More information