ETNA Kent State University
|
|
- Annabelle Watts
- 6 years ago
- Views:
Transcription
1 Electronic Transactions on Numerical Analysis. Volume 35, pp. 2-26, 29. Copyright 29,. ISSN ETNA GAUSSIAN DIRECT QUADRATURE METHODS FOR DOUBLE DELAY VOLTERRA INTEGRAL EQUATIONS ANGELAMARIA CARDONE, IDA DEL PRETE, AND CLAUDIA NITSCH Abstract. In this paper we consider Volterra integral equations with two constant delays. We construct Direct Quadrature methods based on Gaussian formulas, combined with a suitable interpolation technique. We study the convergence and the stability properties of the methods and we carry out some numerical experiments that confirm our theoretical results. Key words. Volterra integral equations, Direct Quadrature method, Gaussian quadrature formulas, convergence, stability AMS subject classifications. 65R2. Introduction. In this paper we consider double delays Volterra integral equations (VIEs) of the type (.) y(t) = f(t) + t τ t τ 2 k(t τ)g(y(τ))dτ, t [τ 2, T], with y(t) = φ(t), t [, τ 2 ], τ, τ 2 R +, where φ(t) is a known function such that (.2) φ(τ 2 ) = f(τ 2 ) + τ 2 τ k(τ 2 τ)g(φ(τ))dτ. We assume that the functions f(t), k(t), and φ(t) are at least continuous on [, T], on [τ, τ 2 ] and on [, τ 2 ], respectively, and that g(y) satisfies the Lipschitz condition. These assumptions ensure existence, uniqueness, and continuity of the solution of (.) [7]. By successively differentiating (.) it is easy to verify that y (l), l =, 2,... presents some points, θ,..., θ Z, of primary discontinuities (θ := τ 2 for y ; θ := τ 2, θ 2 := τ 2 + τ, θ 3 := 2τ 2 for y ;...), and it is continuous in ]lτ 2, T]. Double delay VIEs arise in the mathematical modeling of population dynamics, whose present history depends only on a finite and variable part of the past history. For example, equations of the form (.) model the growth of a population structured by age with a finite life span [, 4]. The numerical treatment of (.) has been carried out only recently and in the specialized literature a few papers can be found on this topic. To our knowledge, the only numerical methods for the equation (.) have been constructed in [6, 7], where Direct Quadrature (DQ) methods based on Newton-Cotes formulas have been proposed. The aim of our research is to extend the class of numerical methods for solving equation (.) to DQ methods based on Gaussian quadrature formulas. It is known that Gaussian formulas ensure a higher order of accuracy and have better stability properties than Newton- Cotes formulas [2]. On the other hand, the use of DQ methods based on Gaussian formulas Received July 7, 28. Accepted for publication May 5, 29. Published online November 6, 29. Recommended by S. Ehrich. Dipartimento di Matematica e Informatica, Università degli Studi di Salerno, Via Ponte don Melillo I-8484 Fisciano (Sa), Italy (ancardone@unisa.it). Dipartimento di Matematica e Applicazioni, R. Caccioppoli, Università degli Studi di Napoli Federico II, Via Cintia, Monte S. Angelo, I-826 Napoli, Italy (ida.delprete@unina.it,clanitsch@yahoo.it). 2
2 22 A. CARDONE, I. DEL PRETE, AND C. NITSCH produces several problems. As a matter of fact, such methods require the knowledge of the solution at some points not belonging to the mesh. In order to overcome this difficulty we have chosen to use an interpolation technique. We have studied the convergence properties of the constructed methods and we have proven that their order of convergence is the minimum between the order of convergence of the Gaussian formula and the degree of accuracy of the interpolating polynomial. The study of numerical stability of our methods has been carried out on the following test equation, introduced in [7]: t τ y(t) = + (λ + µ(t s))y(s)ds, t τ 2 t [τ 2, T], λ, µ R. We have found sufficient conditions under which the numerical solution produced by our method shows the same behavior as the analytical one. In particular we have determined the bound and the limiting value of the numerical solution. Section 2 contains the construction of the method and in Section 3 the convergence analysis is carried out. The numerical stability of our method is treated in Section 4. In Section 5, we report some numerical experiments that confirm the theoretical results stated in Sections 3 and 4. In Section 6, some concluding remarks and future developments are reported. 2. The method. Let Π N = {t j : < t < t < < t N = T } be a partition of the time interval [, T] with constant stepsize h = t j+ t j, j =,...,N, and assume that there exist n and n 2, positive integers, such that (2.) h = τ n = τ 2 n 2. In the following we denote by y j an approximation to the exact solution y(t j ) of (.). Let {ξ k } m k= and {ω k} m k= be the nodes and the weights of an m-point Gaussian quadrature formula on [, ]. Then, the m-point Gaussian quadrature formula on [, h] has nodes {ξ k h} m k= and weights {hω k } m k=, (2.2) h Φ(ξ)dξ h k= m ω k Φ(ξ k h), where Φ(ξ) is any continuous integrand function. The integral equation (.) at the mesh points is (2.3) tj τ y(t j ) = f(t j ) + k(t j τ)g(y(τ))dτ, t j τ 2 j = n 2 +,...,N. The integral in (2.3) can be written as (2.4) tj τ n 2 k(t j τ)g(y(τ))dτ = t j τ 2 r= h k(q r h τ)g(y(t j qr + τ))dτ, where n 2 := n 2 n, q r := n 2 r +. We discretize each of the integrals on [, h] by the quadrature rule (2.2), thus obtaining (2.5) n 2 m y(t j ) f(t j ) + h ω k k(q r h ξ k h)g(y(t j qr + ξ k h)), j = n 2 +,...,N.
3 GAUSSIAN DIRECT QUADRATURE METHODS FOR DELAY VOLTERRA INTEGRAL EQUATIONS 23 The points, t j qr + ξ k h, k =,...,m, r =,...,n 2, do not belong to the mesh Π N. In order to overcome this problem, we adopt an interpolation technique, similar to that used in [5] to discretize a linear delay integro-differential equation. We construct the Lagrange interpolating polynomial P(x) of degree s for the data points, (t j qr s, y j qr s ),..., (t j qr, y j qr ),..., (t j qr+, y j qr+ ), with s := s +, and s, N, that is, (2.6) P(x) = L l (x)y j qr+l, where L l is the lth fundamental Lagrange polynomial with nodes t j qr s,..., t j qr+. We replace y(t j qr +ξ k h) by P(t j qr +ξ k h) in (2.5), thus obtaining the numerical method, (2.7) n 2 m y j = f(t j ) + h ω k k(q r h ξ k h)g(p(t j qr + ξ k h)), j = n 2 +,...,N. In the following it will be useful to observe that L l (t j qr + xh) = P l (x), where P l is the lth fundamental Lagrange polynomial determined by the nodes, s,...,, namely, (2.8) P l (x) = i= s i l x i l i. Thus the method (2.7) can be written equivalently as n 2 s m + y j = f(t j ) + h ω k k(q r h ξ k h)g P l (ξ k )y j qr+l, (2.9) j = n 2 +,...,N. The method (2.7) depends on the parameters s, and m, which have to satisfy some suitable conditions for its applicability. First, in order not to require values of the solution outside [, T], we have to require that s. In addiction, a necessary and sufficient condition to avoid the use of future mesh points (where the numerical solution is not yet known) is (2.) n +. Thus, we have s and n +. Finally, we observe that, if (2.) holds and n +, then the method (2.7) is explicit.
4 24 A. CARDONE, I. DEL PRETE, AND C. NITSCH 3. Convergence analysis. In the previous section we have constructed our method by using two kinds of approximations, one arising from the discretization of the integral in (.) by the Gaussian quadrature formula (2.2), and another one arising from the computation of the unknown values of the solution by the interpolating polynomial (2.6). Therefore, in the convergence analysis, we have to take into account the contributions to the errors of both approximations. In the following we assume that h satisfies the condition (2.), which implies that the discontinuity points θ,..., θ Z of order s+ are all included in the mesh Π N. Furthermore, we assume that (3.) j, r z : either t j qr s,..., t j qr+ [θ z, θ z+ ] t j qr s θ Z or t j qr+ θ hold. Condition (3.) may be satisfied by a suitable choice of s and. Now we are able to prove the following theorem. THEOREM 3.. Let y j be the numerical solution of (.) obtained by the method (2.7) with s and n +. Let p = min(2m, ) and q = max(2m, ), where 2m is the order of the quadrature formula (2.2) and s is the degree of the interpolating polynomial (2.6). Assume that f C s+ ([, T]), k C q ([τ, τ 2 ]), φ C s+ ([, τ 2 ]), g C q (R), and g satisfies the Lipschitz condition. Then, for sufficiently small step size h, the error e j = y(t j ) y j satisfies max e j Ch p, j N for some finite C not depending on h. REMARK 3.2. From the smoothness hypotheses on φ, f, and k, the exact solution y(t) of (.) is at least times continuously differentiable on [θ z, θ z + ], z =,..., Z, and on [, θ ] and [θ Z, T]. From the expression for y (ν) (t), ν =,...,, obtained by successively differentiating (.) with respect to t, it is readily seen that both the left and right limits of y (ν) (t) as t θ z exist and are finite. Proof of Theorem 3.. Take t j [θ z, θ z+ ], with z,..., Z. We have [ n 2 ] h m e j = k(q r h τ)g(y(t j qr + τ))dτ h w k k(q r h ξ k h)g(p(t j qr + ξ k h)). r= k= n 2 We can rewrite the error as e j = (B jr + D jr ), where and D jr := B jr := h h r= k(q r h τ)g(y(t j qr + τ))dτ k(q r h τ)g(p(t j qr + τ))dτ h h or k(q r h τ)g(p(t j qr + τ))dτ m w k k(q r h ξ k h)g(p(t j qr + ξ k h)). k= Thus, letting L be the Lipschitz constant for g, we have B jr L +L h h k(q r h τ) I jr (τ) dτ k(q r h τ) P l (τ/h)(y(t j qr+l) y j qr+l) dτ,
5 GAUSSIAN DIRECT QUADRATURE METHODS FOR DELAY VOLTERRA INTEGRAL EQUATIONS 25 where I jr (τ) = y(t j qr + τ) P l (τ/h)y(t j qr+l) is the interpolation error at t j qr +τ. By condition (3.) and Remark 3.2, there exists C s >, such that I jr (τ) C s h s+, τ [, h]. We observe that the constant C s when s is large has an exponential behaviour, since it depends on the Lebesgue constant. Therefore, (3.2) B jr C s h s+2 + C h e j qr+l, where C s = L C s max k(t) and C = L max max P l(x) max k(t). D jr is t [τ,τ 2] l { s,...,} x [,] t [τ,τ 2] the Gauss-Legendre quadrature error in [, h] for the function k(q r h τ)g(p(t j qr + τ). Then, there exists C m > such that [3, (2.7.2), p. 98] (3.3) D jr C m h 2m+. By (3.2) and (3.3) it follows that n 2 s + e j Cs h s+2 + C h e j qr+l + C m h 2m+ r= Therefore, = (τ 2 τ )( Cs h s+ + C m h 2m) + C h C s,m h p + C h j e i. i= e j C s,m C h hp + C j C h h e i. i= e j qr+l Now we apply the Gronwall-type inequality [2, p. 4], and, since there are no starting errors, we get e j C s,m C h hp e C C h T. Therefore, e j = O(h p ) as h and the theorem follows. Theorem 3. allows us to choose the parameters of the method by balancing efficiency and accuracy, as shown in the following example. EXAMPLE 3.3. The best way to achieve order 2 is the -node Gaussian quadrature formula (m = ) combined with the linear interpolation (s = ). The choices s =, = yield (3.4) (3.5) ω =, ξ = /2, P(x) = t j q r+ x y j qr + x t j q r y j qr+. h h
6 26 A. CARDONE, I. DEL PRETE, AND C. NITSCH 4. Numerical stability. Our next task is to study the stability properties of our numerical methods with respect to test equation introduced in [7], (4.) t τ y(t) = + (λ + µ(t τ))y(τ)dτ, t τ 2 t [τ 2, T], λ, µ R, according to the following definition. DEFINITION 4.. A numerical method is stable with respect to (4.) when its application to (4.) gives a numerical solution behaving like the continous one. It is known that the stability analysis for test equations is the starting point for the investigation of the stability properties of the method for more general equations. We make the following definitions: (4.2) ρ := λ(τ 2 τ ) + µ 2 (τ2 2 τ2 ), (4.3) α := 2µ (λ + µτ ) 2, β := 2µ (λ + µτ 2) 2, Φ := max [,τ 2] φ(t). In the theorem below we summarize some theoretical results about the bound and the limiting value of the analytical solution of (4.) [7]. THEOREM 4.2. Assume that one of the following sets of conditions holds: a) (λ + µτ )(λ + µτ 2 ), ρ <, b) (λ + µτ )(λ + µτ 2 ), α + β <. Then y(t) is bounded for all t and y(t) ( α + β ) + Φ. Moreover, lim y(t) = t + ρ. The numerical solution of (4.) obtained by the method (2.9) is (4.4) n 2 m y j = + h ω k (λ + µh(q r ξ k )) P l (ξ k )y j qr+l, j = n 2 +,...,N. Let us define the sets of indices (4.5) R = {(r, k) : λ + µh(q r ξ k ) }, R 2 = {(r, k) : λ + µh(q r ξ k ) > }, and (4.6) L := {l : P l (ξ k ), k =,..., m} = {( ), 2, 4,..., 2l,...} L 2 := {l : P l (ξ k ) >, k =,..., m} = {,, 3,..., 2l +,...}. According to formula (2.8) all the zeros of P l (x) are integers. Therefore, in (,) P l (x) cannot vanish, that is P l (x) has a constant sign in (, ). Hence, P l (x) assumes the same sign
7 GAUSSIAN DIRECT QUADRATURE METHODS FOR DELAY VOLTERRA INTEGRAL EQUATIONS 27 at all the nodes ξ k, k =,..., m (which belong to (, ); see formula (2.2)) and equalities (4.6) follow. Moreover, we will use α(h) = h w k (λ + µh(q r ξ k )), r,k R (4.7) β(h) = h w k (λ + µh(q r ξ k )). r,k R 2 We set (4.8) Since (4.9) σ := max k=,...,m P l (ξ k ) =, k =,..., m, we obtain l L P l (ξ k ). P l (ξ k ) 2σ +, k =,..., m. Our first result establishes some sufficient conditions for the boundedness of the numerical solution (4.4). THEOREM 4.3. Assume that one of the following set of conditions holds: a) (λ + µτ )(λ + µτ 2 ), (2σ + ) ρ <, b) (λ + µτ )(λ + µτ 2 ), (2σ + )( α(h) + β(h) ) <. Then (4.) (4.) y j Proof. We rewrite (4.4) as y j + Φ, j =,..., N. (2σ + )( α(h) + β(h) ) n 2 m = + δh w k (λ + µh(q r ξ k ))P qr (ξ k )y j n 2 m +h w k (λ + µh(q r ξ k )) P l (ξ k )y j qr+l, l q r where δ = if = n + (implicit method; see Section 2), otherwise δ =. The sum over l q r contains δ terms. a) We consider the case λ + µx > for x [τ, τ 2 ], since the case λ + µx < can be treated similarly. From (4.) it follows that n 2 m ( δθ) y j + h w k (λ + µh(q r ξ k )) P l (ξ k ) y j qr+l, l q r n 2 m where Θ := h w k (λ + µh(q r ξ k )) P qr (ξ k ). Now we proceed step by step. For j < n 2 +, since y j = y(t j ) = ϕ(t j ), then y j < Φ. For j = n 2 +, since δθ <, we have y n2+ + ((2σ + )ρ δθ)φ. δθ
8 28 A. CARDONE, I. DEL PRETE, AND C. NITSCH For j = n 2 + 2, we have ( δθ) y n2+2 + ((2σ + )ρ δθ)max Now assume that Φ < (2σ+)ρ. Then we have y n2+2 Since (2σ+) δθ δθ <, ( + δθ ( Φ, ) (2σ + )ρ δθ + δθ y n2+2 Φ + Continuing with the same procedure, we can prove (2σ + )ρ. ) + ((2σ + )ρ δθ)φ. δθ ( ) 2 (2σ + )ρ δθ Φ. δθ (4.2) y j Φ +, j =,...,N, (2σ + )ρ which is equivalent to (4.), since in this case α(h) + β(h) = ρ. We arrive at the same result when Φ (2σ+)ρ. b) By (4.) we get (4.3) n 2 m ( δθ) y j + h w k λ + µh(q r ξ k ) P l (ξ k ) y j qr+l, l q r n 2 m where Θ := h w k λ + µh(q r ξ k ) P qr (ξ k ). As in the case a), we proceed step by step. We find that y j Φ, j =,..., n 2, and δθ <. Therefore, for j 2 +, we have (4.4) y n2+ + Φ((2σ + )( α(h) + β(h) ) δθ). δθ For j = n 2 + 2, from (4.3) and (4.4), we obtain [ ] (2σ + )( α(h) + β(h) ) δθ y n2+2 + δθ δθ ( ) 2 (2σ + )( α(h) + β(h) ) δθ + Φ δθ Φ + (2σ + )( α(h) + β(h) ), Here we have assumed that Φ /( (2σ + )( α(h) + β(h) )). The opposite case can be treated in an analogous way. By proceeding step by step the theorem follows. In addition to the previous results, we can show that the limit of the numerical solution is the same as for the analytical solution. THEOREM 4.4. Let the conditions of Theorem 4.3 be satisfied. Then lim j y j = ρ,
9 GAUSSIAN DIRECT QUADRATURE METHODS FOR DELAY VOLTERRA INTEGRAL EQUATIONS 29 with ρ given by (4.2). Proof. If l = liminf j y j and l = limsup y j, then there exist two subsequences j {y ju } ju, {y jv } jv such that l = lim u y j u lim v y j v = l, with (4.5) (4.6) n 2 m y ju = + h w k (λ + µh(q r ξ k )) n 2 m y jv = + h w k (λ + µh(q r ξ k )) l=s P l (ξ k )y ju q r+l, l=s P l (ξ k )y jv q r+l. a) If we assume that λ + µτ > and λ + µτ 2 >, then we can split the sum over l in (4.5) into two parts according to the sign of the Lagrange fundamental polynomials, (4.7) y ju n 2 m = + h w k (λ + µh(q r ξ k )) P l (ξ k )y ju q r+l l L n 2 m +h w k (λ + µh(q r ξ k )) l L 2 P l (ξ k )y ju q r+l. By computing the limit inferior of both sides of (4.7) and by recalling that l = lim y j u = liminf y j u, we get u u (4.8) where l + ᾱ(h)l + β(h)l, n 2 m ᾱ(h) = h w k (λ + µh(q r ξ k )) P l (ξ k ), l L n 2 m β(h) = h w k (λ + µh(q r ξ k )) P l (ξ k ). l L 2 In the same way, by considering (4.6) and by computing the limit superior, we find that (4.9) l + ᾱ(h)l + β(h)l, By substracting (4.8) to (4.9), we find that l l ( β(h) ᾱ(h))(l l ). Since β(h) ᾱ(h) (2σ + ) ρ <, we have l = l and, thus, y j admits a limiting value l = l = y, where y = ρ. As a matter of fact, when j + in the (4.4), we have y = + ρy,
10 2 A. CARDONE, I. DEL PRETE, AND C. NITSCH s + n 2 m by recalling that P l (ξ k ) = and observing that h w k (λ + µh(q r ξ k )) = ρ. In a similar way, we are able to prove the theorem when λ + µτ < and λ + µτ 2 <. b) The equation (4.4) can be written as y j = + h w k (λ + µh(q r ξ k )) P l (ξ k )y j qr+l r,k R l L +h w k (λ + µh(q r ξ k )) P l (ξ k )y j qr+l r,k R l L 2 +h w k (λ + µh(q r ξ k )) P l (ξ k )y j qr+l r,k R 2 l L +h w k (λ + µh(q r ξ k )) P l (ξ k )y j qr+l. r,k R 2 l L 2 Therefore, by computing the limit inferior, we find that l + ᾱ(h)l + β(h)l, where (4.2) ᾱ(h) = h w k (λ + µh(q r ξ k )) P l (ξ k ) r,k R l L 2 +h w k (λ + µh(q r ξ k )) P l (ξ k ), r,k R 2 l L β(h) = h w k (λ + µh(q r ξ k )) P l (ξ k ) r,k R l L +h w k (λ + µh(q r ξ k )) P l (ξ k ), r,k R 2 l L 2 By following the same procedure for y jv, we obtain and l + ᾱ(h)l + β(h)l, l l ( β(h) ᾱ(h))(l l ) = ( β(h) + ᾱ(h) )(l l ). If β(h) + ᾱ(h) (2σ + )( β(h) + α(h) ) <, then l = l = l and lim y j = l = y. j By means of Theorems 4.3 and 4.4, we have shown the following theorem. THEOREM 4.5. Assume that one of the following set of conditions holds: a) (λ + µτ )(λ + µτ 2 ), (2σ + ) ρ <, b) (λ + µτ )(λ + µτ 2 ), (2σ + )( α(h) + β(h) ) <. Then the method (2.7) is stable with respect to the test equation (4.). Here we have used Definition 4. of numerical stability.
11 GAUSSIAN DIRECT QUADRATURE METHODS FOR DELAY VOLTERRA INTEGRAL EQUATIONS 2 In order to establish a connection between the behavior of the numerical solution produced by our method and the behavior of the analytical solution, provided that the stepsize h is sufficiently small, we need this final result. THEOREM 4.6. Let α(h) and β(h) be defined by (4.7. Then where α and β are defined by (4.3). lim α(h) = α and lim β(h) = β, h h Proof. Let us suppose that µ >. Then λ + µx < for x < λ/µ. The opposite case can be treated in the same way. Let q r+ h be the rightmost point, such that λ + µq r. Hence R = {(r, k) : r > r, k =,..., m} { ( r, k) : k = k,..., m } and R 2 = { ( r, k) : k =,..., k } {(r, k) : r =,..., r e k =,..., m} for a suitable k. It follows that with and with α(h) = h = n 2 r= r+ k= q r+h τ m m w k (λ + µh(q r ξ k )) + h w k (λ + µh(q r ξ k )) (λ + µx)dx + h k= k m w k (λ + µh(q r ξ k )) k= k = [ (λ + µq r+ h) 2 (λ + µτ ) 2] + h 2µ = α + p(h), p(h) = 2µ (λ + µq r+h) 2 + h k k= m w k (λ + µh(q r ξ k )) k= k m w k (λ + µh(q r ξ k )), k= k r β(h) = h w k (λ + µh(q r ξ k )) + h k = h w k (λ + µh(q r ξ k )) + k= k k= τ 2 q rh m w k (λ + µh(q r ξ k )) (λ + µx)dx = h w k (λ + µh(q r ξ k )) + [ (λ + µτ2 ) 2 (λ + µq r h) 2] 2µ = β q(h), q(h) = 2µ (λ + µq rh) k 2 h w k (λ + µh(q r ξ k )). According to the definition of r, it follows that q r+ h λ/µ and ( k ) m as h. As a consequence k= lim h p(h) = 2µ (λ + µ ( λ/µ))2 =
12 22 A. CARDONE, I. DEL PRETE, AND C. NITSCH and lim q(h) = lim h h 2µ (λ + µq rh) 2 h = lim h 2µ (λ + µq rh) 2 m w k (λ + µh(q r ξ k )) k= q rh q r+h (λ + µx)dx = lim 2µ (λ + µq r+h) 2 =. h REMARK 4.7. Theorem 4.6 shows that, when the stepsize h is sufficiently small, both sufficient conditions of Theorem 4.5 differ from the hypotheses of Theorem 4.2 by the presence of the factor (2σ + ), if σ. However, it is always possible to determine the parameters of the method in such a way that the numerical solution mimics the behavior of the analytical solution. In the following section, we will illustrate how to choose the best parameters of the method. 5. Numerical experiments. Theorems 3. and 4.5 establish that the parameter s influences both convergence and the stability of the method (2.7). For a concrete application of the method, one has to consider the following points: s should be chosen as s = 2m in order to preserve the order of the Gaussian quadrature formula without loss of efficiency; For any stable equation of type (.) there exists a stable method of order p 2; The parameter σ, which determines the stability of the method, is equal to zero for s = and s =, and increases with s, thus limiting the order of the method. According to Theorem 4.5, as ρ or β α approach one (that is the analytical solution is almost unstable), the maximum order of convergence we may expect is lower. In Table 5. we list some practical choices of s and of s with respect to the value of ρ (the same holds for β α, if h is sufficiently small) and the corresponding order p, which guarantee stability of the method. These values are not the maximal ones, since they come from the sufficient but not necessary conditions for stability of Theorem 4.5. Now we illustrate the performance of our methods for equations of the type (5.) y(t) = f(t) + t τ t τ 2 (λ + µ(t τ))g(y(τ))dτ, t [τ 2, T], with τ =.5, τ 2 =., T = 5., λ =., µ =.2. We have chosen the following methods of type (2.7):. order 2 method, with s =, =, m = (3.4) (3.5), TABLE 5. Practical choices of the parameters of the method. ρ s s p [.,.235] 9 [.235,.47] 7 8 [.47,.65] 5 6 [.65,.857] 3 4 [.857, ] 2
13 GAUSSIAN DIRECT QUADRATURE METHODS FOR DELAY VOLTERRA INTEGRAL EQUATIONS 23 TABLE 5.2 Correct digits for problem (5.) (5.2). Gauss Newton-Cotes h p = 2 p = 4 p = 2 p = TABLE 5.3 Correct digits for problem (5.) (5.3). h p = 2 p = order 4 method, with s =, = 3, m = 2, 3. order 6 method, with s =, = 4, m = 3, 4. order 8 methods, with s =, = 7, or s =, = 6, and m = 4. All four methods are explicit. Analogous tests carried out with implicit methods have produced similar results. We have tested the convergence properties of our methods on the folloowing two problems [7]: problem (5.) with (5.2) g(y) =, f(t) s.t. y(t) = e t, φ(t) = e t ; problem (5.) with (5.3) g(y) = ( + y) 2, f(t) s.t. y(t) = sint, φ(t) = sin t; In Tables the number of correct digits values (cd) of the solution of (5.) obtained by the methods and 2, are listed for different values of h. These tables clearly show that our methods produce the desired order according to Theorem 3.. Moreover, in Table 5.2, we have compared our results with those obtained by the only other numerical approach we know, that is, DQ methods based on Newton-Cotes formulas (in particular the trapezoidal and Simpson 3/8 formulas) [7]: the errors of the two families of methods are very similar, in spite of the error of approximation due to the interpolation technique used for the Gaussian formulas. In order to test the sharpness of the estimates of the stability parameters of our methods, we have carried out a large number of numerical experiments. Here the most significant ones
14 24 A. CARDONE, I. DEL PRETE, AND C. NITSCH 2 order 2 5 order FIG. 5.. On the left: λ =, µ =.2, ρ =.95, (2σ + ) ρ =.95, p = 2, s =, h =.. On the right: λ =, µ =.2, ρ =.95, (2σ + ) ρ =.54, p = 6, s =, h =.. 2 order 8 4 x 4 order FIG On the left: λ =, µ =, ρ =.25, (2σ + ) ρ =.6, p = 8, s =, h =.. On the right: λ =, µ =.2, ρ =.95, (2σ + ) ρ = 2.84, p = 6, s =, h =.. are reported. We consider the performances of the method (2.7) when applied to the test problem (5.4) y(t) = f(t) + t τ t τ 2 (λ + µ(t τ))y(τ)dτ, t [τ 2, T], with τ =.5, τ 2 =. and f such that y(t) = t sin t. In our examples λ and µ satisfy the hypotheses of Theorem 4.2 for the stability of the analytical solution. Problems of type (5.4) may be assimilated to the test equation (4.) with the simplified assumption f(t) f(), as it is usually done in the formulation of test equations; see [6] and references therein. The numerical solution (*-) and the analytical one (-) are compared in Figures , for different values of the parameters λ and µ and for different choices of the parameters of the method. Figures show tests related to the case (λ+µτ )(λ+µτ 2 ) >. When the hypothesis a) of Theorem 4.5 is satisfied, the numerical solution behaves like the analytical one (left plots of Figures 5. and 5.2), while in the other case instability could arise, both for s = and for s = (right plots of Figures ). Similar results are found when (λ + µτ )(λ + µτ 2 ) <. This is illustrated in Figure 5.3. We emphasize that we dispose of high order and stable methods, as shown by test of Figure 5.2, where an order 8 method is
15 GAUSSIAN DIRECT QUADRATURE METHODS FOR DELAY VOLTERRA INTEGRAL EQUATIONS 25 2 order 4 2 x 5 order FIG λ = 8, µ =, β α =.65, (2σ + )(β(h) α(h)) =.95, p = 4, s =, h =.. On the right: λ = 8, µ =, β α =.65, (2σ + )(β(h) α(h)) =.6, p = 6, s =, h =.. used. Our tests confirm what we have noticed in Remark 4.7 about the influence of the factor 2σ + on the numerical stability: for example, in the right plot of Figure 5., ρ < but σ =.83, so that (2σ + ) ρ > and instability occurs. We observe that the estimate on the stability parameter (2σ + ) ρ or (2σ + )(β(h) α(h)) is quite sharp, as shown for example by Figure 5.3, where (2σ + )(β(h) α(h)) varies between the values.95 and.6. In the example illustrated by the right plot of Figure 5.2, ρ is very close to and so, in principle, we expect an order of accuracy not exceeding 2 for a stable method; see Table 5.. Therefore, in order to increase the order, we have applied a cubic spline interpolation technique (which is supposed to have order 4). This method proved to be stable; on the other hand also our DQ method of order 4 with the Lagrange interpolating polynomial is stable and guarantees the same accuracy. 6. Concluding remarks. In this paper we constructed Direct Quadrature methods based on Gaussian quadrature formulas for equation (.). In order to solve the problem of the evaluation of the solution at points not belonging to the mesh, an interpolation technique has been used. In Theorem 3. we have shown that the order of convergence of the method depends both on the order of convergence of the Gaussian formula and on the degree of accuracy of the interpolating polynomial. In order to complete the study of the method proposed, we analyzed the stability with respect to a class of significant test equations introduced in [7]. We found sufficient conditions for numerical stability. These conditions are such that, if the starting problem satisfies the conditions of Theorem 4.2, then it is easy to determine the parameters of the method that secure the stability of the numerical solution. The numerical experiments clearly confirm the theoretical results and show the sharpness of the estimates of the stability parameters. Finally we note that, in order to increase the order of convergence when ρ is close to one, it may be possible to consider another interpolation technique such as one using cubic splines. This new approach, together with other approximation techniques, requires complete analysis of convergence and numerical stability. We intend to dedicate a future work to this topic.
16 26 A. CARDONE, I. DEL PRETE, AND C. NITSCH REFERENCES [] D. BREDA, C. CUSULIN, M. IANNELLI, S. MASET, AND R. VERMIGLIO, Stability analysis of age-structured population equations by pseudospectral differencing methods, J. Math. Biol., 54 (27), pp [2] H. BRUNNER AND P. J. VAN DER HOUWEN, The Numerical Solution of Volterra Equations, Vol. 3, North- Holland, Amsterdam, 986. [3] P. J. DAVIS AND P. RABINOWITZ, Methods of Numerical Integration, Second ed., Academic Press, New York, 984. [4] M. IANNELLI, Mathematical Theory of Age-Structured Populations Dynamics, Applied Mathematics Monographs C.N.R., Vol. 7, Giardini, Pisa, 995. [5] T. LUZYANINA, K. ENGELBORGHS, AND D. ROOSE, Computing stability of differential equations with bounded distributed delays, Numer. Algorithms, 34 (23), pp [6] E. MESSINA, E. RUSSO, AND A. VECCHIO, A stable numerical method for Volterra integral equations with discontinous kernel, J. Math. Anal. Appl., 337 (28), pp [7], A convolution test equation for double delay integral equations, J. Comput. Appl. Math., 228 (29), pp
A convolution test equation for double delay integral equations
Journal of Computational and Applied Mathematics 228 (2009) 589 599 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationSuperconvergence analysis of multistep collocation method for delay Volterra integral equations
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 4, No. 3, 216, pp. 25-216 Superconvergence analysis of multistep collocation method for delay Volterra integral equations
More informationA Note on the Asymptotic Stability of Linear Volterra Difference Equations of Convolution Type
Trinity University Digital Commons @ Trinity Mathematics Faculty Research Mathematics Department -1-007 A Note on the Asymptotic Stability of Linear Volterra Difference Equations of Convolution Type Saber
More informationON CERTAIN CLASSES OF CURVE SINGULARITIES WITH REDUCED TANGENT CONE
ON CERTAIN CLASSES OF CURVE SINGULARITIES WITH REDUCED TANGENT CONE Alessandro De Paris Università degli studi di Napoli Federico II Dipartimento di Matematica e Applicazioni R. Caccioppoli Complesso Monte
More informationCollocation and iterated collocation methods for a class of weakly singular Volterra integral equations
Journal of Computational and Applied Mathematics 229 (29) 363 372 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationPart IB Numerical Analysis
Part IB Numerical Analysis Definitions Based on lectures by G. Moore Notes taken by Dexter Chua Lent 206 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after
More informationSMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm. Abstract
SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS Kai Diethelm Abstract Dedicated to Prof. Michele Caputo on the occasion of his 8th birthday We consider ordinary fractional
More informationJim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More information1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel
1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel Xiaoyong Zhang 1, Junlin Li 2 1 Shanghai Maritime University, Shanghai,
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationStability constants for kernel-based interpolation processes
Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report 59 Stability constants for kernel-based interpolation processes Stefano De Marchi Robert Schaback Dipartimento
More informationOn quadrature formulas for oscillatory evolutionary problems
On quadrature formulas for oscillatory evolutionary problems Angelamaria Cardone, Dajana Conte, Raffaele D Ambrosio, Beatrice Paternoster Abstract Gaussian-type quadrature rules for oscillatory integrand
More informationThe generalized Euler-Maclaurin formula for the numerical solution of Abel-type integral equations.
The generalized Euler-Maclaurin formula for the numerical solution of Abel-type integral equations. Johannes Tausch Abstract An extension of the Euler-Maclaurin formula to singular integrals was introduced
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationJacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation
Wei et al. SpringerPlus (06) 5:70 DOI 0.86/s40064-06-3358-z RESEARCH Jacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation Open Access
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 29 Almost Done! No
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationmultistep methods Last modified: November 28, 2017 Recall that we are interested in the numerical solution of the initial value problem (IVP):
MATH 351 Fall 217 multistep methods http://www.phys.uconn.edu/ rozman/courses/m351_17f/ Last modified: November 28, 217 Recall that we are interested in the numerical solution of the initial value problem
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 35, pp. 148-163, 2009. Copyright 2009,. ISSN 1068-9613. ETNA SPHERICAL QUADRATURE FORMULAS WITH EQUALLY SPACED NODES ON LATITUDINAL CIRCLES DANIELA
More informationA FOURTH-ORDER METHOD FOR NUMERICAL INTEGRATION OF AGE- AND SIZE-STRUCTURED POPULATION MODELS
A FOURTH-ORDER METHOD FOR NUMERICAL INTEGRATION OF AGE- AND SIZE-STRUCTURED POPULATION MODELS Mimmo Iannelli, Tanya Kostova and Fabio Augusto Milner Abstract In many applications of age- and size-structured
More information8.1 Introduction. Consider the initial value problem (IVP):
8.1 Introduction Consider the initial value problem (IVP): y dy dt = f(t, y), y(t 0)=y 0, t 0 t T. Geometrically: solutions are a one parameter family of curves y = y(t) in(t, y)-plane. Assume solution
More informationGALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS
GALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS GEORGIOS AKRIVIS AND CHARALAMBOS MAKRIDAKIS Abstract. We consider discontinuous as well as continuous Galerkin methods for the time discretization
More informationSummer School on Delay Differential Equations and Applications
Summer School on Delay Differential Equations and Applications Dobbiaco (BZ), Italy, June 26 30, 2006 The numerical solution of delay differential equations Page 1 of 211 M. Zennaro Dipartimento di Matematica
More informationA Nodal Spline Collocation Method for the Solution of Cauchy Singular Integral Equations 1
European Society of Computational Methods in Sciences and Engineering (ESCMSE) Journal of Numerical Analysis, Industrial and Applied Mathematics (JNAIAM) vol. 3, no. 3-4, 2008, pp. 211-220 ISSN 1790 8140
More informationQuarter-Sweep Gauss-Seidel Method for Solving First Order Linear Fredholm Integro-differential Equations
MATEMATIKA, 2011, Volume 27, Number 2, 199 208 c Department of Mathematical Sciences, UTM Quarter-Sweep Gauss-Seidel Method for Solving First Order Linear Fredholm Integro-differential Equations 1 E. Aruchunan
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More informationNumerische Mathematik
Numer. Math. (28) 9:43 65 DOI.7/s2-7-25-7 Numerische Mathematik The superconvergence of Newton Cotes rules for the Hadamard finite-part integral on an interval Jiming Wu Weiwei Sun Received: 2 May 27 /
More informationNUMERICAL METHOD FOR THE MIXED VOLTERRA-FREDHOLM INTEGRAL EQUATIONS USING HYBRID LEGENDRE FUNCTIONS
Conference Applications of Mathematics 215, in honor of the birthday anniversaries of Ivo Babuška (9), Milan Práger (85), and Emil Vitásek (85) J. Brandts, S. Korotov, M. Křížek, K. Segeth, J. Šístek,
More informationNUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS
TWMS J Pure Appl Math V5, N2, 24, pp22-228 NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS S ASADI, AH BORZABADI Abstract In this paper, Haar wavelet benefits are applied to the delay
More information1 Solutions to selected problems
Solutions to selected problems Section., #a,c,d. a. p x = n for i = n : 0 p x = xp x + i end b. z = x, y = x for i = : n y = y + x i z = zy end c. y = (t x ), p t = a for i = : n y = y(t x i ) p t = p
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 33 Almost Done! Last
More informationConvergence and Stability of a New Quadrature Rule for Evaluating Hilbert Transform
Convergence and Stability of a New Quadrature Rule for Evaluating Hilbert Transform Maria Rosaria Capobianco CNR - National Research Council of Italy Institute for Computational Applications Mauro Picone,
More informationPOLYNOMIAL INTERPOLATION IN NONDIVISION ALGEBRAS
Electronic Transactions on Numerical Analysis. Volume 44, pp. 660 670, 2015. Copyright c 2015,. ISSN 1068 9613. ETNA POLYNOMIAL INTERPOLATION IN NONDIVISION ALGEBRAS GERHARD OPFER Abstract. Algorithms
More informationNumerical Analysis. A Comprehensive Introduction. H. R. Schwarz University of Zürich Switzerland. with a contribution by
Numerical Analysis A Comprehensive Introduction H. R. Schwarz University of Zürich Switzerland with a contribution by J. Waldvogel Swiss Federal Institute of Technology, Zürich JOHN WILEY & SONS Chichester
More informationOn the Lebesgue constant of subperiodic trigonometric interpolation
On the Lebesgue constant of subperiodic trigonometric interpolation Gaspare Da Fies and Marco Vianello November 4, 202 Abstract We solve a recent conjecture, proving that the Lebesgue constant of Chebyshev-like
More informationPeriodic solutions of weakly coupled superlinear systems
Periodic solutions of weakly coupled superlinear systems Alessandro Fonda and Andrea Sfecci Abstract By the use of a higher dimensional version of the Poincaré Birkhoff theorem, we are able to generalize
More informationRemarks on Bronštein s root theorem
Remarks on Bronštein s root theorem Guy Métivier January 23, 2017 1 Introduction In [Br1], M.D.Bronštein proved that the roots of hyperbolic polynomials (1.1) p(t, τ) = τ m + m p k (t)τ m k. which depend
More informationKey Words: cardinal B-spline, coefficients, moments, rectangular rule, interpolating quadratic spline, hat function, cubic B-spline.
M a t h e m a t i c a B a l a n i c a New Series Vol. 24, 2, Fasc.3-4 Splines in Numerical Integration Zlato Udovičić We gave a short review of several results which are related to the role of splines
More informatione st f (t) dt = e st tf(t) dt = L {t f(t)} s
Additional operational properties How to find the Laplace transform of a function f (t) that is multiplied by a monomial t n, the transform of a special type of integral, and the transform of a periodic
More informationMathematics for chemical engineers. Numerical solution of ordinary differential equations
Mathematics for chemical engineers Drahoslava Janovská Numerical solution of ordinary differential equations Initial value problem Winter Semester 2015-2016 Outline 1 Introduction 2 One step methods Euler
More informationA Product Integration Approach Based on New Orthogonal Polynomials for Nonlinear Weakly Singular Integral Equations
Acta Appl Math (21) 19: 861 873 DOI 1.17/s144-8-9351-y A Product Integration Approach Based on New Orthogonal Polynomials for Nonlinear Weakly Singular Integral Equations M. Rasty M. Hadizadeh Received:
More informationAPPM/MATH Problem Set 6 Solutions
APPM/MATH 460 Problem Set 6 Solutions This assignment is due by 4pm on Wednesday, November 6th You may either turn it in to me in class or in the box outside my office door (ECOT ) Minimal credit will
More informationIterated Defect Correction with B-Splines for a Class of Strongly Nonlinear Two-Point Boundary Value Problems
American Review of Mathematics and Statistics June 2016, Vol. 4, No. 1, pp. 31-44 ISSN: 2374-2348 (Print), 2374-2356 (Online) Copyright The Author(s). All Rights Reserved. Published by American Research
More informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
More information12.0 Properties of orthogonal polynomials
12.0 Properties of orthogonal polynomials In this section we study orthogonal polynomials to use them for the construction of quadrature formulas investigate projections on polynomial spaces and their
More informationOptimal stopping time formulation of adaptive image filtering
Optimal stopping time formulation of adaptive image filtering I. Capuzzo Dolcetta, R. Ferretti 19.04.2000 Abstract This paper presents an approach to image filtering based on an optimal stopping time problem
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Linear Multistep methods Gustaf Söderlind and Carmen Arévalo Numerical Analysis, Lund University Textbooks: A First Course in the
More informationTHE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS. Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands
THE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands 1. Introduction This paper deals with initial value problems for delay
More information1 Solutions to selected problems
Solutions to selected problems Section., #a,c,d. a. p x = n for i = n : 0 p x = xp x + i end b. z = x, y = x for i = : n y = y + x i z = zy end c. y = (t x ), p t = a for i = : n y = y(t x i ) p t = p
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE612B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationEXISTENCE RESULTS FOR AN ELLIPTIC DIRICHLET PROBLEM
LE MATEMATICHE Vol. LXVI (2) Fasc. I pp. 33 4 doi:.448/2.66..3 EXISTENCE RESULTS FOR AN ELLIPTIC DIRICHLET PROBLEM GIUSEPPINA D AGUÌ - GIOVANNI MOLICA BISCI The main purpose of this paper is to present
More informationInitial value problems for ordinary differential equations
Initial value problems for ordinary differential equations Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 IVP
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationNumerical Solution of Two-Dimensional Volterra Integral Equations by Spectral Galerkin Method
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 159-174 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 Numerical Solution of Two-Dimensional Volterra
More informationOrdinary Differential Equations
Chapter 13 Ordinary Differential Equations We motivated the problem of interpolation in Chapter 11 by transitioning from analzying to finding functions. That is, in problems like interpolation and regression,
More informationExistence And Uniqueness Of Mild Solutions Of Second Order Volterra Integrodifferential Equations With Nonlocal Conditions
Applied Mathematics E-Notes, 9(29), 11-18 c ISSN 167-251 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Existence And Uniqueness Of Mild Solutions Of Second Order Volterra Integrodifferential
More informationNumerical Integration (Quadrature) Another application for our interpolation tools!
Numerical Integration (Quadrature) Another application for our interpolation tools! Integration: Area under a curve Curve = data or function Integrating data Finite number of data points spacing specified
More informationAN EXAMPLE OF FUNCTIONAL WHICH IS WEAKLY LOWER SEMICONTINUOUS ON W 1,p FOR EVERY p > 2 BUT NOT ON H0
AN EXAMPLE OF FUNCTIONAL WHICH IS WEAKLY LOWER SEMICONTINUOUS ON W,p FOR EVERY p > BUT NOT ON H FERNANDO FARRONI, RAFFAELLA GIOVA AND FRANÇOIS MURAT Abstract. In this note we give an example of functional
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... find the approximate solutions of derivative (first- and second-order) and antiderivative (definite integral only). Numerical Differentiation
More informationExplicit polynomial expansions of regular real functions by means of even order Bernoulli polynomials and boundary values
Journal of Computational and Applied Mathematics 176 (5 77 9 www.elsevier.com/locate/cam Explicit polynomial expansions of regular real functions by means of even order Bernoulli polynomials and boundary
More informationCONVERGENCE ANALYSIS OF THE JACOBI SPECTRAL-COLLOCATION METHODS FOR VOLTERRA INTEGRAL EQUATIONS WITH A WEAKLY SINGULAR KERNEL
COVERGECE AALYSIS OF THE JACOBI SPECTRAL-COLLOCATIO METHODS FOR VOLTERRA ITEGRAL EQUATIOS WITH A WEAKLY SIGULAR KEREL YAPIG CHE AD TAO TAG Abstract. In this paper, a Jacobi-collocation spectral method
More informationEnergy-Preserving Runge-Kutta methods
Energy-Preserving Runge-Kutta methods Fasma Diele, Brigida Pace Istituto per le Applicazioni del Calcolo M. Picone, CNR, Via Amendola 122, 70126 Bari, Italy f.diele@ba.iac.cnr.it b.pace@ba.iac.cnr.it SDS2010,
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More informationA REACHABLE THROUGHPUT UPPER BOUND FOR LIVE AND SAFE FREE CHOICE NETS VIA T-INVARIANTS
A REACHABLE THROUGHPUT UPPER BOUND FOR LIVE AND SAFE FREE CHOICE NETS VIA T-INVARIANTS Francesco Basile, Ciro Carbone, Pasquale Chiacchio Dipartimento di Ingegneria Elettrica e dell Informazione, Università
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationIntegration. Topic: Trapezoidal Rule. Major: General Engineering. Author: Autar Kaw, Charlie Barker.
Integration Topic: Trapezoidal Rule Major: General Engineering Author: Autar Kaw, Charlie Barker 1 What is Integration Integration: The process of measuring the area under a function plotted on a graph.
More informationINTERPOLATION. and y i = cos x i, i = 0, 1, 2 This gives us the three points. Now find a quadratic polynomial. p(x) = a 0 + a 1 x + a 2 x 2.
INTERPOLATION Interpolation is a process of finding a formula (often a polynomial) whose graph will pass through a given set of points (x, y). As an example, consider defining and x 0 = 0, x 1 = π/4, x
More informationExamination paper for TMA4215 Numerical Mathematics
Department of Mathematical Sciences Examination paper for TMA425 Numerical Mathematics Academic contact during examination: Trond Kvamsdal Phone: 93058702 Examination date: 6th of December 207 Examination
More informationSection 6.6 Gaussian Quadrature
Section 6.6 Gaussian Quadrature Key Terms: Method of undetermined coefficients Nonlinear systems Gaussian quadrature Error Legendre polynomials Inner product Adapted from http://pathfinder.scar.utoronto.ca/~dyer/csca57/book_p/node44.html
More informationNumerical Solution of Integral Equations in Solidification and Melting with Spherical Symmetry
Numerical Solution of Integral Equations in Solidification and Melting with Spherical Symmetry V. S. Ajaev and J. Tausch 2 Southern Methodist University ajaev@smu.edu 2 Southern Methodist University tausch@smu.edu
More informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More informationGALERKIN AND RUNGE KUTTA METHODS: UNIFIED FORMULATION, A POSTERIORI ERROR ESTIMATES AND NODAL SUPERCONVERGENCE
GALERKIN AND RUNGE KUTTA METHODS: UNIFIED FORMULATION, A POSTERIORI ERROR ESTIMATES AND NODAL SUPERCONVERGENCE GEORGIOS AKRIVIS, CHARALAMBOS MAKRIDAKIS, AND RICARDO H. NOCHETTO Abstract. We unify the formulation
More informationInterpolation. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Interpolation 1 / 34
Interpolation Escuela de Ingeniería Informática de Oviedo (Dpto. de Matemáticas-UniOvi) Numerical Computation Interpolation 1 / 34 Outline 1 Introduction 2 Lagrange interpolation 3 Piecewise polynomial
More informationConvergence rates of derivatives of a family of barycentric rational interpolants
Convergence rates of derivatives of a family of barycentric rational interpolants J.-P. Berrut, M. S. Floater and G. Klein University of Fribourg (Switzerland) CMA / IFI, University of Oslo jean-paul.berrut@unifr.ch
More informationDepartment of Applied Mathematics Preliminary Examination in Numerical Analysis August, 2013
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August, 013 August 8, 013 Solutions: 1 Root Finding (a) Let the root be x = α We subtract α from both sides of x n+1 = x
More informationAccuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations
COMMUICATIOS I COMPUTATIOAL PHYSICS Vol. 5, o. -4, pp. 779-79 Commun. Comput. Phys. February 009 Accuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations Tao
More informationBindel, Spring 2012 Intro to Scientific Computing (CS 3220) Week 12: Monday, Apr 16. f(x) dx,
Panel integration Week 12: Monday, Apr 16 Suppose we want to compute the integral b a f(x) dx In estimating a derivative, it makes sense to use a locally accurate approximation to the function around the
More informationNumerical Methods for Differential Equations Mathematical and Computational Tools
Numerical Methods for Differential Equations Mathematical and Computational Tools Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 Part 1. Vector norms, matrix norms and logarithmic
More informationCS 257: Numerical Methods
CS 57: Numerical Methods Final Exam Study Guide Version 1.00 Created by Charles Feng http://www.fenguin.net CS 57: Numerical Methods Final Exam Study Guide 1 Contents 1 Introductory Matter 3 1.1 Calculus
More informationExamination paper for TMA4125 Matematikk 4N
Department of Mathematical Sciences Examination paper for TMA45 Matematikk 4N Academic contact during examination: Anne Kværnø a, Louis-Philippe Thibault b Phone: a 9 66 38 4, b 9 3 0 95 Examination date:
More informationA Survey of Numerical Methods in Fractional Calculus
A Survey of Numerical Methods in Fractional Calculus Fractional Derivatives in Mechanics: State of the Art CNAM Paris, 17 November 2006 Kai Diethelm diethelm@gns-mbh.com, k.diethelm@tu-bs.de Gesellschaft
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationScientific Computing: Numerical Integration
Scientific Computing: Numerical Integration Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Fall 2015 Nov 5th, 2015 A. Donev (Courant Institute) Lecture
More informationEngg. Math. II (Unit-IV) Numerical Analysis
Dr. Satish Shukla of 33 Engg. Math. II (Unit-IV) Numerical Analysis Syllabus. Interpolation and Curve Fitting: Introduction to Interpolation; Calculus of Finite Differences; Finite Difference and Divided
More informationApproximate solution of linear integro-differential equations by using modified Taylor expansion method
ISSN 1 746-7233, Engl, UK World Journal of Modelling Simulation Vol 9 (213) No 4, pp 289-31 Approximate solution of linear integro-differential equations by using modified Taylor expansion method Jalil
More informationASYMPTOTICS AND EVALUATIONS OF FPT DENSITIES THROUGH VARYING BOUNDARIES FOR GAUSS-MARKOV PROCESSES
Scientiae Mathematicae Japonicae Online, e-28, 157 182 157 ASYMPTOTICS AND EVALUATIONS OF FPT DENSITIES THROUGH VARYING BOUNDARIES FOR GAUSS-MARKOV PROCESSES Amelia G. Nobile, Enrica Pirozzi, Luigi M.
More informationNumerische Mathematik
Numer. Math. (997) 76: 479 488 Numerische Mathematik c Springer-Verlag 997 Electronic Edition Exponential decay of C cubic splines vanishing at two symmetric points in each knot interval Sang Dong Kim,,
More informationChapter 2 Interpolation
Chapter 2 Interpolation Experiments usually produce a discrete set of data points (x i, f i ) which represent the value of a function f (x) for a finite set of arguments {x 0...x n }. If additional data
More informationThe collocation method for ODEs: an introduction
058065 - Collocation Methods for Volterra Integral Related Functional Differential The collocation method for ODEs: an introduction A collocation solution u h to a functional equation for example an ordinary
More informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationDynamic Systems and Applications 13 (2004) PARTIAL DIFFERENTIATION ON TIME SCALES
Dynamic Systems and Applications 13 (2004) 351-379 PARTIAL DIFFERENTIATION ON TIME SCALES MARTIN BOHNER AND GUSEIN SH GUSEINOV University of Missouri Rolla, Department of Mathematics and Statistics, Rolla,
More informationStability analysis of delayed SIR epidemic models with a class of nonlinear incidence rates
Published in Applied Mathematics and Computation 218 (2012 5327-5336 DOI: 10.1016/j.amc.2011.11.016 Stability analysis of delayed SIR epidemic models with a class of nonlinear incidence rates Yoichi Enatsu
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 20, 2014 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
More information