arxiv: v2 [math.pr] 6 Mar 2013

Transcription

1 Strong Convergence of uler Approximations of Stochastic Differential quations with Delay under Local Lipschitz Condition arxiv: v math.pr 6 Mar 013 Chaman Kumar and Sotirios Sabanis School of Mathematics, University of dinburgh, dinburgh, H9 3JZ, U.K. October 16, 018 Abstract The strong convergence of uler approximations of stochastic delay differential equations is proved under general conditions. The assumptions on drift and diffusion coefficients have been relaxed to include polynomial growth and only continuity in the arguments corresponding to delays. Furthermore, the rate of convergence is obtained under one-sided and polynomial Lipschitz conditions. Finally, our findings are demonstrated with the help of numerical simulations. Keywords: stochastic delay differential equations, uler approximations, strong convergence, rate of convergence, local Lipschitz condition. 1 Introduction In modeling many real world phenomena, the future states of the system depend not only on the present state but also on its past state(s). The models based on stochastic delay differential equations (SDDs) could be used in such situations. They have found numerous applications in various fields, for example, in communications, physics, biology, ecology, economics and finance. One could refer to 1-10 and references therein. The authors would like to thank Lukas Szpruch for his useful discussions. 1

2 In many applications, one finds SDDs for which solutions can not be obtained explicitly and hence this necessitates the development of numerical schemes. Over the past few years, many authors have shown their interest in studying strong convergence properties of numerical schemes ofsdds, tonameafew, Recently therehasbeenagrowinginterest inweak convergence of numerical schemes as well, see for example, 17, 18 and 19. We also mention a paper by 0 on pathwise approximations and its relation to strong convergence. For the purpose of this article, we study strong convergence of uler approximations of SDDs due to its wide applicability. We exploit the fact that SDDs can be regarded as special cases of stochastic differential equations (SDs) with random coefficients. To this end, first we establish strong convergence results of uler schemes of SDs with random coefficients under mild conditions. This, in turn, allows us to prove the convergence in the case of SDDs under more relaxed conditions than those existing in the literature. More precisely, to the best of our knowledge, strong convergence has been proved so far by assuming linear growth and local Lipschitz conditions in both variables corresponding to delay and non-delay arguments, see for example 16. Whereas, we relax these conditions by assuming polynomial growth and continuity in the argument corresponding to delays. Moreover no smoothness condition on the initial data is assumed. Therefore, the setting in this article is more general than those present in the literature, see for example, 11, 14, 15, 16 and 0. In addition, it is shown here that the rate of convergence is one-fourth when the drift coefficient satisfies one-sided Lipschitz and polynomial Lipschitz in the non-delay and delay arguments respectively. This result is in agreement with the findings of 1. We conclude this section with the introduction of some basic notation. The norm of a vector x R d and the Hilbert-Schmidt norm of a matrix A R d m are respectively denoted by x and A. The transpose of a vector x R d is denoted by x T and the inner product of two vectors x,y R d is denoted by x,y := x T y. The integer part of a real number x is denoted by x. We denote the indicator function of set A by ½ A. Moreover, L p = L p (Ω,F,P) denotes the space of random variables X with a norm X p := ( X p) 1/p < for p > 0. Finally, P denotes the predictable σ-algebra on R + Ω and B(V), the σ-algebra of Borel sets of topological spaces V.

3 Main Result Let (Ω,{F t } {t 0},F,P) be a filtered probability space satisfying the usual conditions, i.e. the filtration is increasing and right continuous. Let {W(t)} {t 0} be an m-dimensional Wiener martingale. Furthermore, it is assumed that β(t,y 1,...,y k,x) and α(t,y 1,...,y k,x) are B(R + ) B(R d k ) B(R d )-measurable functions and take values in R d and R d m respectively. For a fixed T > 0, let the stochastic delay differential equation (SDD) on (Ω,{F t } {t 0},F,P) be defined as follows, dx(t) = β(t,y(t),x(t))dt +α(t,y(t),x(t))dw(t), t 0,T, X(t) = ξ(t), t H,0, (1) where{ξ(t) : H t 0} CF b 0 ( H,0;R d )forsomeh > 0andY(t) := (X(δ 1 (t)),...,x(δ k (t))). The delay parameters δ 1 (t),...,δ k (t) are increasing functions of t and satisfy H δ j (t) t τ τ for some τ > 0 and j = 1,...,k. Further, for every n 1 and t 0,T, the uler scheme of SDD (1) is given by dx n (t) = β(t,y n (t),x n (κ n (t)))dt+α(t,y n (t),x n (κ n (t)))dw(t), t 0,T, X n (t) = ξ(t), t H,0, () where Y n (t) := (X n (δ 1 (t)),...,x n (δ k (t))) and κ n is defined by κ n (t) := n(t ) n + (3) with the observation that for equation () one takes = 0 in (3). We note that two popular cases of delay viz. δ i (t) = t τ and δ i (t) = t τ τ can be addressed by our findings. This type of delay parameters have been extensively discussed and found wide applications in the literature. The reader could consult the following, for example, 1-6, 10 and references therein. Remark 1. We note that the uler scheme () defines approximations to SDDs in an explicit way without a discretization of the delay terms. From our main theorem and corollaries on convergence of scheme (), one could easily obtain results on convergence of this scheme with discretized delay terms. Therefore for matters of notational simplicity, we choose the former approach. Let y := (y 1,...,y k ). We make the following assumptions for our result. 3

4 C-1. There exist constants G > 0 and l > 0 such that for any t 0,T, β(t,y,x) + α(t,y,x) G(1+ y l + x ) for all x R d and y R d k. C-. For every R > 0, there exists a constant K R > 0 such that for any t 0,T, x z,β(t,y,x) β(t,y,z) α(t,y,x) α(t,y,z) K R x z whenever x, y, z < R. C-3. The functions β(t,y,x) and α(t,y,x) are continuous in y uniformly in x from compacts, i.e. for every R > 0 and t 0,T, sup { β(t,y,x) β(t,y,x) + α(t,y,x) α(t,y,x) } 0 as y y. x R The conditions C-1 and C- are sufficient for existence and uniqueness of solution of SDD (1) and the uler scheme ()(see 1, ). We state the main result of this paper in the following theorem. Theorem 1. Suppose C-1 to C-3 hold, then the uler scheme () converges to the true solution of SDD (1) in L p -sense, i.e. lim n sup X(t) X n (t) p = 0 0 t T for all p > 0. Remark. We note that we can assume, without loss of generality, that T is a multiple of τ. If not, then SDD (1) and its M scheme () are defined for T > T with T = Nτ, where N is a positive integer. The results given in this report are then recovered for the original SDD (1) by choosing drift and diffusion parameters as β1 {t T} and α1 {t T}. First, we develop requisite theory for SDs with random coefficients in the following section and then prove Theorem 1 in Section 4. 4

5 3 SDs with Random Coefficients Let b(t,x) and σ(t,x) be P B(R d )-measurable functions which take values in R d and R d m respectively. Suppose 0 < t 1 T, then let us consider an SD with random coefficients given by dx(t) = b(t,x(t))dt+σ(t,x(t))dw(t), t,t 1, (4) with initial value X( ) which is an almost surely finite F t0 -measurable random variable. For every n 1, let b n (t,x) and σ n (t,x) be P B(R d )-measurable functions which take values in R d and R d m respectively. Then for any t,t 1, we define the following uler scheme corresponding to the SD (4) as, dx n (t) = b n (t,x n (κ n (t)))dt+σ n (t,x n (κ n (t)))dw(t) (5) with initial value X n ( ), which is an F t0 -measurable, almost surely finite, random variable and κ n as defined in (3). We make the following assumptions. A-1. There exists a constant c > 0 and non-negative random variables {M n } n 1 and M with M p sup n 1 Mn p < N for every p > 0 and some N := N(p) > 0 such that, almost surely, b(t,x) + σ(t,x) c(m + x ) b n (t,x) + σ n (t,x) c(m n + x ) for any t 0,T and x R d. A-. For every R > 0, there exists an F t0 -measurable random variable C R such that, almost surely, for any t,t 1, x y,b(t,x) b(t,y) σ(t,x) σ(t,y) C R x y for all x, y R. Moreover, there exists a function f : R + R + such that 5

6 lim P(C R > f(r)) = 0. R A-3. For every R > 0 and p > 0, lim n t1 sup { b n (t,x) b(t,x) p + σ n (t,x) σ(t,x) p }dt x R A-4. The initial values of SD (4) and its uler scheme (5) satisfy = 0. lim n X( ) X n ( ) p = 0 for every p > 0. A-5. The initial values of SD (4) and its uler scheme (5) have bounded p-th moments, i.e. there exists L := L(p) > 0 such that X( ) p sup X n ( ) p < L n 1 for every p > 0. We note that assumptions A-1 and A- are sufficient for existence and uniqueness of solution of SD (4) (see 3 and 4). Before proving the main result (Theorem ) of this section, we establish some lemmas. Lemma 1. Suppose that A-1 and A-5 hold, then for some K := K(p,T,L) > 0, and sup X(t) p < K t t 1 sup sup X n (t) p < K (6) n 1 t t 1 for every p > 0. Proof. First one chooses any p >. Then for the true solution {X(t)} t t0,t 1, by using Hölder s 6

7 inequality, one obtains, { t X(t) p 3 p 1 X( ) p + t p 1 b(s,x(s)) p ds+ p } σ(s,x(s))dw(s). Thuson taking supremumover,u for someu,t 1 andby applyingburkholder-davis-gundy inequality with constant c := c(p) > 0, one observes that u sup X(t) p 3 { X(t p 1 0 ) p +T p 1 b(s,x(s)) p ds t u Now A-1 and A-5 give t u } + c σ(s,x(s)) p ds. sup X(t) p t u { X(t p 1 0 ) p + T p 1 + c u } c p { } pds M + X(s) 3 { X(t p 1 0 ) p + p 1 T p 1 + c c p TN + p 1 T p 1 + c u } c p sup X(r) ds p r s which on the application of Gronwall s inequality yields sup X(t) p t t 1 < 3 p 1{ X( ) p + p 1 T p 1 + c c p TN } e 6p 1 (T p 1 + c)c pt. For any p >, we adopt similar arguments for {X n (t)} {t t0,t 1 } in order to prove (6) with the remark that A-5 guarantees that K does not depend on n. Then one observes that the application of Hölder s inequality completes the proof. For every R > 0 and n 1, let us define the stopping times, τ R := inf{t : X(t) R}, σ nr := inf{t : X n (t) R} and ν nr := τ R σ nr, (7) where inf =. We prove below a very useful lemma of this article as it is used in the penultimate section in order to recover the rate of convergence. 7

8 Lemma. Let us assume that A-1 and A-5 hold. Then for every p > 0. lim n t1 X n (s ν nr ) X n (κ n (s ν nr )) p ds = 0 Proof. In order to prove the result, one chooses p. Then one immediately writes, t1 t1 X n (s ν nr ) X n (κ n (s ν nr )) p ds = s νnr p + σ n (r,x n (κ n (r)))dw(r) ds κ n(s ν nr ) s νnr κ n(s ν nr ) b n (r,x n (κ n (r)))dr which implies on applying Hölder s inequality t1 X n (s ν nr ) X n (κ n (s ν nr )) p ds t1 ( s νnr ) p 1 s ν nr κ n (s ν nr ) p 1 b n (r,x n (κ n (r))) p dr ds κ n(s ν nr ) t1 s νnr + p 1 σ n (r,x n (κ n (r)))dw(r) p ds ( p 1 t1 n) κ n(s ν nr ) s νnr + p 1 t1 b n (r,x n (κ n (r))) p dr κ n(s ν nr ) s νnr κ n(s ν nr ) ds σ n (r,x n (κ n (r)))dw(r) p ds. (8) One finds bounds for the integrand of the first term of (8) as follows, s νnr s νnr b n (r,x n (κ n (r))) p dr c p ( M n + X n (κ n (r)) ) p dr κ n(s ν nr ) s νnr κ n(s ν nr ) c p p 1 ( M n p + X n (κ n (r)) p )dr κ n(s ν nr ) ( c p p 1 M n p + sup cp p 1 ( M n p + n X n (r) p) s ν nr κ n (s ν nr ) r t 1 ) sup X n (r) p r t 1 cp p 1 (K +N). (using Lemma 1) (9) n 8

9 For the second term of (8), one writes, s νnr p s σ n (r,x n (κ n (r)))dw(r) = κ n(s ν nr ) ( ) p(p 1) p/ s s κ n (s) (p )/ κ n(s) κ n(s) p ½ κn(s ν nr ),s ν nr σ n (r,x n (κ n (r)))dw(r) ½ κn(s ν nr ),s ν nr σ n (r,x n (κ n (r))) p dr which on the application of A-1 and Lemma 1 yields s νnr κ n(s ν nr ) σ n (r,x n (κ n (r)))dw(r) p ( p(p 1) ) p/ 1 n p/cp p 1 (K +N). (10) Substituting (9) and (10) in (8) yields t1 X n (s ν nr ) X n (κ n (s ν nr )) p ds T(K +N)4p 1 c p { 1 n p/ n p/ + ( ) p(p 1) p/ } O(n p/ ). (11) Then by letting n on both sides, one obtains lim n t1 X n (s ν nr ) X n (κ n (s ν nr )) p ds for all p. Then one observes that the application of Hölder s and Jensen s inequalities completes the proof. = 0 In order to prove the following theorem, a similar approach to the one adopted by 5 for SDs is followed here, but in a more general context. Theorem. Let A-1 to A-5 hold and suppose there exists a unique solution {X n (t)} {t0 t t 1 } of the uler scheme (5). Then, the uler scheme (5) converges to the SD (4) in L p -sense, i.e. lim n sup X(t) X n (t) p t t 1 = 0 (1) for all p > 0. Proof. For every R > 0 and n 1, we consider stopping times as defined in (7). First we fix p > 4. 9

10 Then sup X(t) X n (t) p t t 1 + sup X(t) X n (t) p ½ {τr t 1 or σ nr t 1 or C R >f(r)} t t 1 sup t t 1 X(t ν nr ) X n (t ν nr) p ½ {CR f(r)}. (13) To estimate the first term of (13), one can use Young s inequality for q > p (1/p + 1/q = 1) and η > 0 to obtain the following. sup t t 1 X(t) X n (t) p ½ {τr t 1 or σ nr t 1 or C R >f(r)} { ηp q q 1 sup X(t) q t t 1 { ηp q q K + q p qη p/(q p) ηp q sup X(t) X n (t) q t t 1 + q p qη p/(q p)p(τ R t 1 orσ nr t 1 or C R > f(r)) } + sup X n (t) q t t 1 + q p { } qη p/(q p) P(τ R t 1 )+P(σ nr t 1 )+P(C R > f(r)) } X(τ R ) p X n (σ nr ) p + +P(C R > f(r)), R p R p which becomes sup X(t) X n (t) p ½ {τr t 1 or σ nr t 1 or C R >f(r)} t t 1 ηp q q K + q p qη p/(q p) { K R p +P(C R > f(r)) }. (14) Further for estimating the second term of (13), let e n (s) := X(s) X n (s) for s,t 1, then one could write de n (s) = b n (s)ds+ σ n (s)dw(s) where b n (s) := b(s,x(s)) b n (s,x n (κ n (s))) and σ n (s) := σ(s,x(s)) σ n (s,x n (κ n (s))). Thus, on the application of Itô s formula, one obtains m m d e n (s) = e n (s), b n (s) ds+ σ n i (s) ds+ e n (s), σ n i (s) dwi (s) i=1 i=1 10

11 which on integrating over s,t ν nr for some t,t 1 yields e n (t ν nr ) = e n ( ) + + m i=1 t νnr t νnr e n (s), b n (s) ds+ m i=1 e n (s), σ i n (s) dwi (s) t νnr σ i n(s) ds where σ i n and W i denote the i-th column of (d m)-matrix σ and the i-th element of Wiener (m 1)-vector W respectively for all i {1,..., m}. Further, after raising the power to p/ for some p > 4 and applying Hölder s and Burkholder-Davis-Gundy inequalities, one obtains the following, sup e n (t ν nr ) p ½ {CR f(r)} t u +(8T) p +(4mT) p ½ {CR f(r)} m i=1 p e n ( ) p u νnr ½ {CR f(r)} + 3p 4 m p T p 4 4 c e n (s), b n (s) p ds m i=1 u νnr ½ {CR f(r)} σ i n(s) p ds u νnr e n (s), σ n i ds (s) p = I +II +III +IV (15) where c := c(p) is the constant of the Burkholder-Davis-Gundy inequality. In order to estimate II, one observes that the application of A- and Cauchy-Schwarz inequality yields, e n (s), b n (s) = X(s) X n (κ n (s)),b(s,x(s)) b(s,x n (κ n (s))) + X(s) X n (κ n (s)),b(s,x n (κ n (s))) b n (s,x n (κ n (s))) + X n (κ n (s)) X n (s),b(s,x(s)) b n (s,x n (κ n (s))) (16) C R X(s) X n (κ n (s)) + X(s) X n (κ n (s)) b(s,x n (κ n (s))) b n (s,x n (κ n (s))) + X n (s) X n (κ n (s)) b(s,x(s)) b n (s,x n (κ n (s))) 11

12 and on further application of Young s inequality and A-1, this becomes e n (s), b n (s) (C R +1) e n (s) +(C R +1) X n (s) X n (κ n (s)) + 1 b(s,x n(κ n (s))) b n (s,x n (κ n (s))) + X n (s) X n (κ n (s)) { c( X(s) + X(κ n (s)) +M +M n ) }. Therefore raising power p/ on both sides, we get e n (s), b n (s) p p 3{ (C R +1) p en (s) p +(C R +1) p Xn (s) X n (κ n (s)) p + b(s,x n (κ n (s))) b n (s,x n (κ n (s))) p } + X n (s) X n (κ n (s)) p { c( X(s) + X(κn (s)) +M +M n ) }p for every s,u ν nr. Hence one obtains the following estimate on using Hölder s inequality, u νnr II := (3T) p ½ {CR f(r)} (3T) p { (f(r)+1) p u e n (s), b n (s) p ds sup e n (r ν nr ) p ½ CR f(r) ds r s t1 +(f(r)+1) p X n (s ν nr ) X n (κ n (s ν nr )) p ds t1 + b(s ν nr,x n (κ n (s ν nr ))) b n (s ν nr,x n (κ n (s ν nr ))) p ds t1 + p+1/ Tc p X n (s ν nr ) X n (κ n (s ν nr )) p ds } K +N. (17) Now in order to estimate III of (15), one observes that for every i {1,...,m}, the assumption A- yields σ i n(s) p = σ i (s,x(s)) σ i n(s,x n (κ n (s))) p 3 p 1{ C p R e n (s) p +C p R X n (s) X n (κ n (s)) p + σ i (s,x n (κ n (s))) σ i n(s,x n (κ n (s))) p} 1

13 which implies III := (4mT) p + m i=1 3 p 1 (4mT) p { ½ {CR f(r)} mf(r) p u u νnr sup r s σ i n(s) p ds e n (r ν nr ) p ½ {CR f(r)} t1 +mf(r) p X n (s ν nr ) X n (κ n (s ν nr )) p ds m } t1 σ i (s ν nr,x n (κ n (s ν nr ))) σn i (s ν nr,x n (κ n (s ν nr ))) p ds. (18) i=1 Finally, to estimate IV, observes that for i {1,...,m}, an application of assumption A- yields ds e n (s), σ i n(s) p e n (s) p σ i n (s) p = en (s) p σ i (s,x(s)) σ i n(s,x n (κ n (s))) p which by using Young s inequality becomes e n (s), σ i n(s) p { 3 p C p 4 R e n (s) p +C p 4 R e n (s) p Xn (s) X n (κ n (s)) p 3p + e n (s) p σ i (s,x n (κ n (s))) σ i n (s,x n(κ n (s))) p {( ) 3C p 4 R +1 e n (s) p +C p 4 R X n (s) X n (κ n (s)) p + σ i (s,x n (κ n (s))) σ i n(s,x n (κ n (s))) p}. } Thus IV := 3p 4 m p T p 4 4 c 3(p ) 3 p m p T p 4 4 c + m u νnr ½ {CR f(r)} e n (s), σ n i ds (s) p i=1 { ) m (3f(R) p u 4 +1 sup e n (r ν nr ) p ½ {CR f(r)} r s t1 +mf(r) p 4 X n (s ν nr ) X n (κ n (s ν nr )) p ds m } t1 σ i (s ν nr,x n (κ n (s ν nr ))) σn(s ν i nr,x n (κ n (s ν nr ))) p ds.(19) i=1 ds 13

14 Hence using (17), (18) and (19) in (15) and then after applying Gronwall s inequality, one obtains { sup e n (t ν nr ) p ½ {CR f(r)} e K 1 T p e n ( ) p t t 1 + K t1 X n (s ν nr ) X n (κ n (s ν nr )) p ds + K t1 3 b(s ν nr,x n (κ n (s ν nr ))) b n (s ν nr,x n (κ n (s ν nr ))) p ds m + K } t1 5 σ i (s ν nr,x n (κ n (s ν nr ))) σn(s ν i nr,x n (κ n (s ν nr ))) p ds i=1 t1 + K 4 X n (s ν nr ) X n (κ n (s ν nr )) p ds (0) where constants K 1, K, K 3, K 4 and K 5 are appropriately defined and depend explicitly on f(r) but not on R and they are independent of n. Thus one can choose η sufficiently small and R sufficiently large such that for ǫ > 0 (however small), ηp q q K < ǫ 3, q p K qη p/(q p) R p < ǫ 3 and q p qη p/(q p)p(c R > f(r)) < ǫ 3. Therefore, substituting equations (14) and (0) in (13) and then using assumption A-3, A-4 and Lemma, oneobtains(1)forp > 4. TheapplicationofHölder sinequalitycompletestheproof. 4 Proof of Main Result In this section, we shall prove Theorem 1 by considering SDD (1) as a special case of an SD with random coefficients. In particular, we set b(t,x) := β(t,y(t),x),σ(t,x) := α(t,y(t),x),b n (t,x) := β(t,y n (t),x),σ n (t,x) := α(t,y n (t),x)(1) for any t 0,T and x R d. Proof of theorem 1. We shall prove the result using an induction method so as to show that, lim n sup X(t) X n (t) p = 0 () (i 1)τ t iτ 14

15 for every i {1,...,N} and for all p > 0. Case: t 0,τ. When t 0,τ, then SDD (1) and uler scheme () become SD (4) and scheme (5) with = 0, t 1 = τ, X(0) = X n (0) = ξ(0) and coefficients given by (1). Also one observes that A-1 to A-5 hold due to C-1 to C-3. In particular, A-1 is a consequence of C-1 with M = M n = 1 + Ψ l, where Ψ := sup t 0,τ (ξ(δ 1 (t)),...,ξ(δ k (t))), which is uniformly bounded due to the fact that ξ CF b 0. Further, A- is a consequence of C- since one observes that C R is deterministic and in fact C R = f(r), where f(r) = K Ψ ½ {0<R<Ψ} + K R ½ {R Ψ} and therefore P(C R > f(r)) = 0 for any R > 0. Finally, A-3, A-4 and A-5 hold trivially. Therefore, for i = 1, equation () holds due to Theorem and Lemma 1 for = 0 and t 1 = τ. For inductive arguments, we assume that when i = r, i.e. t (r 1)τ, rτ for some r {1,...,N 1}, equation () is satisfied and Lemma 1 holds for = (r 1)τ and t 1 = rτ. Then we claim that when i = r +1, i.e. t rτ,(r +1)τ, equation () holds and Lemma 1 is true for = rτ and t 1 = (r +1)τ. Case: t rτ,(r+1)τ. When t rτ,(r+1)τ, then SDD (1) and its uler scheme () become SD (4) and scheme (5) with = rτ, t 1 = (r + 1)τ, X( ) = X(rτ), X n ( ) = X n (rτ) and coefficients given by (1). Verify A-1. Consider for t rτ,(r +1)τ and x R d, then due to C-1, b(t,x) + σ(t,x) = β(t,y(t),x) + α(t,y(t),x) G(M + x ) and b n (t,x) + σ n (t,x) = β(t,y n (t),x) + α(t,y n (t),x) G(M n + x ) where M := 1+sup rτ t (r+1)τ Y(t) l and M n := 1+sup rτ t (r+1)τ Y n (t) l which are bounded in L p for every p > 0 because of Lemma 1 and the inductive assumptions. Verify A-. For every R > 0, x, z R and t rτ,(r +1)τ, from C-, one obtains, x z,b(t,x) b(t,z) = x z,β(t,y(t),x) β(t,y(t),z) C R x z σ(t,x) σ(t,z) = α(t,y(t),x) α(t,y(t),z) C R x z 15

16 where the random variable C R is given by C R := K R 1 ΩR + K j+1 1 {Ωj+1 \Ω j } j=r with Ω j := {ω Ω : sup t rτ,(r+1)τ Y(t) j}. Oneobserves that C R is an F rτ -measurable random variable. Further one takes f(r) := K R for every R > 0 so that P(C R > f(r)) = P(C R > K R ) = 1 P(C R K R ) ( ) 1 P(Ω R ) = P Y(t) > R 0 as R. (3) sup rτ t<(r+1)τ Verify A-3. For R > 0 and x R, one observes that Y(t) Y n (t) for every p > 0 when n. Therefore from C-3, P 0 as Y(t) Y n (t) p 0 sup β(t,y n (t),x) β(t,y(t),x) p P 0 and sup α(t,y n (t),x) α(t,y(t),x) p P 0. x R x R Furthermore, one observes that sequences { } { } sup β(t,y n (t),x) β(t,y(t),x) p and sup α(t,y n (t),x) α(t,y(t),x) p x R n 1 x R n 1 are uniformly integrable since they are bounded in L q for any q > 1, sup β(t,y n (t),x) β(t,y(t),x) pq = sup β(t,y n (t),x) β(t,y(t),x) pq x R x R pq 1 sup β(t,y n (t),x) + pq pq 1 sup β(t,y(t),x) pq x R x R 6 pq 1 G pq{ 1+ Y n (t) lpq + R pq} +6 pq 1 G pq{ 1+ Y(t) lpq + R pq} 6 pq 1 G pq { 1+K + R pq} and similarly, for the sequence {sup x R α(t,y n (t),x) α(t,y(t),x) p } n 1. Therefore, sup x R { β(t,y n (t),x) β(t,y(t),x) p + α(t,y n (t),x) α(t,y(t),x) p} 0 16

17 as n due to Dominated Convergence Theorem which also yields lim n Thus A-3 is satisfied. (r+1)τ rτ sup { b(s,x) b n (s,x) p + σ(s,x) σ n (s,x) p }ds x R Verify A-4. This holds due to the inductive assumptions. Verify A-5. This follows due to the inductive assumptions. FinallyLemma1holdsfor = rτ andt 1 = (r+1)τ duetothefactthata-1holdsfort rτ,(r+1)τ and x R d and A-5 is true for = rτ. This completes the proof. = 0. 5 Rate of Convergence In this section, we shall recover the rates of convergence of the uler scheme () under different set of assumptions. First we state below the relevant assumptions. C-4. There exist constants C > 0 and l 1 > 0 such that for any t 0,T, x x,β(t,y,x) β(t,y,x ) α(t,y,x) α(t,y,x ) C x x β(t,y,x) β(t,y,x) + α(t,y,x) α(t,y,x) C(1+ y l 1 + y l 1 ) y y for all x,x R d and y,y R d k. C-5. There exist constants C > 0 and l > 0 such that for any t 0,T, β(t,y,x) β(t,y,x ) C(1+ x l + x l ) x x for all x,x R d and y R d k. In Corollary 1, we prove that under C-1 and C-4, the rate of convergence of uler scheme () is one-fourth. On the other hand, in Corollary we prove that if one makes assumption C-5 in addition to C-1 and C-4, then the classical rate of convergence (one-half) can be recovered (see also 11, 14, 15, 16 and 0). Corollary 1. Suppose C-1 and C-4 hold, then for all p > 0, the uler scheme () converges to the true solution (1) in L p -sense with rate 1/4, i.e. 17

18 sup X(t) X n (t) p Cn p 4 0 t T where the constant C > 0 does not depend on n. Proof. In equation (14), we choose η = n p, R = n q p (q p), q > p 4, then ηp q q K = pkq n p, q q p K (q p)k q p qη p/(q p) = n p Rp and q qη p/(q p)p(c R > f(r)) 0 (4) where we choose C R C and f(r) C. Then one can show that, sup X(t) X n (t) p Cn p 4 (5) (k 1)τ t kτ for each k {1,...,N}. As seen before, SDD (1) corresponds to the ordinary SD (4) with X(0) = X n (0) = ξ(0), Y(t) = Y n (t) = ξ(δ(t)) for all t 0,τ. Therefore, one obtains the rate as in (5) for k = 1, which is 1/4 instead of 1/, due to the last term of (17), that is also the last term in (0), the estimate in (11) and the fact that τ b(s ν nr,x n (κ n (s ν nr ))) b n (s ν nr,x n (κ n (s ν nr ))) p ds 0 0 and τ σ i (s ν nr,x n (κ n (s ν nr ))) σn(s ν i nr,x n (κ n (s ν nr ))) p ds 0 0 for every i {1,...,m}. One then assumes that (5) holds for k = r, i.e. t (r 1)τ,rτ for some r {1,...,N 1}, so as to show that it also holds for t rτ,(r+1)τ. Due to inductive assumption, the initial data of SD (4) satisfy X(rτ) X n (rτ) p Cn p 4 (6) 18

19 and due to (1), C-4, Lemma 1 and inductive assumptions, one obtains the following estimates, (r+1)τ b(s ν nr,x n (κ n (s ν nr ))) b n (s ν nr,x n (κ n (s ν nr ))) p ds rτ (r+1)τ = β(s,y(s),x n (κ n (s))) β(s,y n (s),x n (κ n (s))) p ½ {s νnr }ds rτ (r+1)τ rτ 3 p 1 Cn p. ( 1+ Y(s) l 1 + Y n (s) l ) p Y(s) Yn (s) p ds (1+K) (r+1)τ rτ Y(s) Y n (s) p ds (7) One similarly estimates (r+1)τ σ i (s ν nr,x n (κ n (s ν nr ))) σn i (s ν nr,x n (κ n (s ν nr ))) p ds Cn p (8) rτ for every i {1,...,m}. Therefore, (5) holds due to (0) by taking into account (11), (4), (6), (7) and (8) for t rτ,(r+1)τ and p 4. Then the application of Hölder s inequality completes the proof. Corollary. Suppose C-1, C-4 and C-5 hold, then for all p > 0, the uler scheme () converges to the true solution (1) in L p -sense with rate 1/, i.e. sup X(t) X n (t) p Cn p 0 t T where the constant C > 0 does not depend on n. Proof. In order to prove this, one observes that the estimates of the rate of convergence of terms of II in (17) is improved if C-5 is used along with C-4. For this, one could replace the right hand side of (16) by the following, e n (s), b n (s) = e n (s),b(s,x(s)) b(s,x n (s)) + e n (s),b(s,x n (s)) b(s,x n (κ n (s))) + e n (s),b(s,x n (κ n (s))) b n (s,x n (κ n (s))) 19

20 which on using (1) becomes e n (s), b n (s) = e n (s),β(s,y(s),x(s)) β(s,y(s),x n (s)) + e n (s),β(s,y(s),x n (s)) β(s,y(s),x n (κ n (s))) + e n (s),β(s,y(s),x n (κ n (s))) β(s,y n (s),x n (κ n (s))). By using C-4 and Cauchy-Schwartz and Young s inequalities, one obtains e n (s), b n (s) (C +1) e n (s) + 1 β(s,y(s),x n(s)) β(s,y(s),x n (κ n (s))) + 1 β(s,y(s),x n(κ n (s))) β(s,y n (s),x n (κ n (s))) which due to C-4 and C-5 yields e n (s), b n (s) (C +1) e n (s) + C {1+ X n(s) l + X n (κ n (s))) l } X n (s) X n (κ n (s))) + C {1+ Y(s) l 1 + Y n (s) l 1 } Y(s) Y n (s). Therefore, one estimates II in (17) by the following, II := (3T) p ½ {CR f(r)} u νnr e n (s), b n (s) p ds { u (96T) p (C +1) p sup e n (s) ds p rτ rτ t s ( )p C u + {1+ X n (s) l + X n (κ n (s))) l } p X n (s) X n (κ n (s))) p ds rτ ( )p C } u + {1+ Y(s) l 1 + Y n (s) l 1 } p Y(s) Y n (s) p ds rτ which on the application of Hölder s inequality and Lemma 1 yields II (96T) p ( C + )p ˆK { (C +1) p (r+1)τ rτ u rτ sup e n (s) ds p rτ t s X } n (s) X n (κ n (s))) p + Y(s) Y n (s) p ds. (9) 0

21 where ˆK is a constant which does not depend on n. Therefore, in (15), one replaces the estimate (17) of II by the estimate (9) whereas III and IV remain the same. Then the first term of (9) is incorporated with similar terms of III and IV so as to apply Gronwall s inequality. Furthermore, one observes that the second and third terms of (9) are of (improved) order O(n p/ ) due to Lemma and inductive assumptions respectively. Thus, by adopting same arguments as adopted in the proof of Corollary 1 with η = n p, one obtains the desired rate (1/). 6 Numerical xamples In this section, we shall illustrate our findings with the help of numerical examples. We consider the following SDD given by dz(t) = az(t)+b{z(t τ)} l 1 dt+β 1 +β Z(t)+β 3 {Z(t τ)} l dw(t) for t 0,τ (30) with initial data ξ(t) when t τ,0, τ > 0 is a fixed delay, l 1,l > 0, a,b,β 1,β,β 3 R. When t 0,τ, then explicit solution of SD (30) is given by, Z(t) = Φ 0,t {ξ(0)+ t Φ 1 0 0,s {b{ξ(s τ)}l 1 β (β 1 +β 3 {ξ(s τ)} l )}ds + t 0 } Φ 1 0,s {β 1 +β 3 {ξ(s τ)} l }dw(s), (31) where, ({ Φ 0,t = exp a β } ) t+β W(t). And when t τ,τ, then the explicit solution is given by Z(t) = Φ τ,t {Z(τ)+ t τ Φ 1 τ,s b{z(s τ)}l 1 β (β 1 +β 3 {Z(s τ)} l )ds + t τ Φ 1 τ,s (β 1 +β 3 {Z(s τ)} l )dw(s) } (3) 1

22 where, ({ Φ τ,t = exp a β } ) (t τ)+β (W(t) W(τ). For uler scheme, we divide the interval 0,τ into sub-intervals of width τ N for some positive integer N. Then solution at kh-th grid is given by Z n ((k +1)h) = az n (kh)+b{ξ(kh τ)} l 1 h + β 1 +β Z n (kh)+β 3 {ξ(kh τ)} l (W((k +1)h) W(kh)) with ξ(t) for t τ,0, k = 1,..., N. We use this solution to find the solution of the uler scheme in the interval τ,τ. We take values of the parameters of SDD (30) given in Table 1. From Table, we observe that as step size is decreased, the error is decreased too. In Figure 1, the reference line has a slope of 0.5. Table 1: Values of Parameters p τ a b β 1 β β 3 ξ(t) for t τ, t+1 Table : ffect of decreasing step size on rror ( Z(T) Zn (T) ) h l 1 = 1/,l = 1/ l 1 = 1,l = 1 l 1 =,l = slope

23 Figure 1: Strong order of convergence of the uler scheme. (a) l 1 = 1,l = 1 (b) l 1 =,l = 3 References 1 Arriojas, M., Hu, Y., Mohammed, S.-. and Pap, G A Delayed Black and Scholes Formula, Stochastic Analysis and Applications, 5(): urich, C. W. and Milton, J. G Noise-induced transitions in human postural sway. Physical Review, 54(6): Federico, S. and Øksendal B. K Optimal stopping of stochastic differential equations with delay driven by a Lêvy noise. Potential Analysis, 34: Federico, S A stochastic control problem with delay arising in a pension fund model. Finance and Stochastics, 15(3): Goldobin, D., Rosenblum, M. and Pikovsky, A Coherence of noisy oscillators with delayed feedback. Physica A: Statistical Mechanics and its Applications, 37(1): Huber, D. and Tsimring, L. S Dynamics of an ensemble of noisy bistable elements with global time delayed coupling. Physical Review Letters, 91(6): Mao, X Stochastic Differential quations and Applications. Horwood Publishing, Chichester. 3

24 8 Mao, X. and Sabanis, S. 01. Delay geometric Brownian motion in financial option valuation. Stochastics: An Interntional Journal of Probability and Stochastic Processes, ifirst: McWilliams, N. and Sabanis, S Arithmetic Asian options under stochastic delay models. Applied Mathematical Finance, 18(5): Øksendal, B., Sulem, A. and Zhang, T Optimal control of stochastic delay equations and time-advanced backward stochastic differential equations. Advances in Applied Probability, 43: Baker, C. T. H. and Buckwar, Numerical analysis of explicit one-step methods for stochastic delay differential equations. LMS Journal of Computation and Mathematics, 3: Gyöngy, I. and Sabanis, S. 01. A note on uler approximation for stochastic differential equations with delay. ArXiv, 13 Hu,Y., Mohammed, S.-. A. and Yan, F Discrete-time approximations of stochastic delay equations: The Milstein scheme. Annals of Probability, 3(1A): Küchler, U. and Platen, Strong discrete time approximation of stochastic differential equations with time delay. Mathematics and Computers in Simulation, 54: Mao, X Numerical solutions of stochastic functional differential equations. LMS Journal of Computation and Mathematics, 6: Mao, X. and Sabanis, S Numerical solutions of stochastic differential delay equations under local Lipschitz condition. Journal of Computational and Applied Mathematics, 151: Buckwar,., Kuske, R., Mohammed, S.-. and Shardlow, T Weak convergence of the uler scheme for stochastic differential delay equations. LMS journal of Computation and Mathematics, 11: Buckwar,. and Shardlow, T Weak approximation of stochastic differential delay equations. IMA Journal of Numerical Analysis, 5(1):

25 19 Küchler, U. and Platen,. 00. Weak discrete time approximation of stochastic differential equations with time delay. Mathematics and Computers in Simulation, 59(6): Kloeden, P.. and Neuenkirch, A The path-wise convergence of approximation schemes for stochastic differential equations. LMS journal of Computation and Mathematics, 10: Gyöngy, I A note on ulers approximations, Potential Analysis, 8: Renesse, M. K. Von and Scheutzow, M xistence and uniqueness of solutions of stochastic functional differential equations. Random Operators and Stochastic quations, 18: Gyöngy, I. and Krylov, N. V On stochastic equations with respect to semimartingales. Stochastics, 4(1): Krylov, N. V A simple proof of the existence of a solution of Itô s equation with monotone coefficients. Theory of Probability and Its Applications, 35: Higham, D. J., Mao, X. and M. Stuart. 00. Strong convergence of uler-type methods for non-linear stochastic differential equations. SIAM Journal on Numerical Analysis, 40(3):