1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel
|
|
- Hilary Douglas
- 5 years ago
- Views:
Transcription
1 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, China 2 Taiyuan University of Science and Technology, Taiyuan, China 1 Corresponding author 1 zxyo20@163.com, 2 lijunlin9726@163.com (Received 7 November 2013; received in revised form 1 November 2014; accepted 20 November 2014) Abstract. The Jacobi pseudo-spectral Galerkin method for the Volterra integro-differential equations of the second kind with a weakly singular kernel is proposed in this paper. We provide a rigorous error analysis for the proposed method, which indicates that the numerical errors (in the, -norm and the -norm) will decay exponentially provided that the source function is sufficiently smooth. Numerical examples are given to illustrate the theoretical results. Keywords: Volterra integro-differential equation, Jacobi pseudo-spectral method, weakly singular kernel, convergence. 1. Introduction In practical applications one frequently encounters the Volterra integro-differential equations of the second kind with a weakly singular kernel of the form: =()() +() + ( ) (, )(), 0<, 0<<1. (1) 0 With the given initial condition (0) =. Where the unknown function () is defined in 0< <. (), () are two given source functions and (, ) is a given kernel. Equations of this type arise as model equations for describing turbulent diffusion problems. The numerical treatment of the Volterra integro-differential Eq. (1) is not simple, mainly due to the fact that the solutions of Eq. (1) usually have a weak singularity at = 0, as discussed in [1], the second derivative of the solution () behaves like ()~. We point out that for Eq. (1) without the singular kernel (i.e., = 0) spectral methods and the corresponding error analysis have been provided recently [2]; see also [3] and [4] for spectral methods to Volterra integral equations and pantograph-type delay differential equations. In both cases, the underlying solutions are smooth. In this work, we will consider a special case, namely, the exact solutions of Eq. (1) are smooth (see also [5]). In this case, the collocation method and product integration method can be applied directly. But the main approach used there is the spectral-collocation method which is similar to a finite-difference approach. Consequently, the corresponding error analysis is more tedious as it does not fit in a unified framework. However, with a finite-element type approach, as will be performed in this work, it is natural to put the approximation scheme under the general Jacobi-Galerkin type framework. As demonstrated in the recent book of Shen etc. [16], there is a unified theory with Jacobi polynomials to approximate numerical solutions for differential and integral equations. It is also rather straightforward to derive the pseudo-spectral Jacobi-Galerkin method from the corresponding continuous version. The relevant convergence theories under the unified framework, as will be seen from Section 4, are cleaner and more reasonable than those obtained in [7]. The purpose of this work is to provide numerical methods for the second kind Volterra integro differential equations based on pseudo-spectral Galerkin methods. Spectral methods are a class of JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
2 techniques used in applied mathematics and scientific computing to numerically solve certain partial differential equations (PDEs) (see e.g. [8-10] and the references therein), often involving the use of the Fast Fourier Transform. Where applicable, spectral methods have excellent error properties, with the so called exponential convergence being the fastest possible. The paper is organized as follows. In Section 2, we introduce the Jacobi pseudo-spectral Galerkin approaches for the Volterra integro-differential Eq. (3). Some preliminaries and useful lemmas are provided in Section 3. In Section 4, the convergence analysis is given. We prove the error estimates in the -norm and, -norm. The numerical experiments are carried out in Section 5, which will be used to verify the theoretical results obtained in Section 4. The final section contains conclusions. 2. Jacobi pseudo-spectral Galerkin method In this section, we formulate the Jacobi pseudo-spectral Galerkin schemes for problem Eq. (1) For this purpose, let, =(1 ) (1+) be a weight function in the usual sense, for, > 1., (), =0, 1,, denote the Jacobi polynomials. The set of Jacobi polynomials {, ()} forms a complete, ( 1,1)-orthogonal system. Before using pseudo-spectral methods, we need to restate problem Eq. (1). The usual way (see [1]) to deal with the original problem is: writing =()+ ( ) (, )(), Eq. (1) is equivalent to a linear Volterra integral equations of the second kind with respect to, : () = + {()() + ()}, () = () + ( ) (, )(). For the sake of applying the theory of orthogonal polynomials conveniently, by the linear transformation: = (1 + ), = 2 Letting: (1 + ). 2 (2) (1 +) (1 +) () =, () =, 2 2 () = + {()() + ()}, 2 () = () + 2 ( ) (, )(). (1 +) () =, 2 Λ = [ 1,1], (3) The weak form of Eq. (3) is to find,, (Λ), (Λ), such that: (, ), = + {()() + ()}, 2,, (, ), =() + 2 ( ) (, )(),,,,, (Λ), (Λ), (4) 3808 JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
3 where (.,.), denotes the usual inner product in the, -space. Now, let be any positive integer and (Λ) be the set of all algebraic polynomials of degree at most. Obviously, the Jacobi polynomials, (),, (),...,, () are the basis functions of (Λ). Next, we denote the collocation points by { }, which is the set of (+ 1) Jacobi Gauss point. We also define the Jacobi interpolating polynomial, (Λ), satisfying:, ( ) =( ), 0. It can be written as an expression of the form:, () =( )( ), (5) where ( ) is the Lagrange interpolation basis function associated with the Jacobi collocation points { }. Now we describe the Jacobi pseudo-spectral Galerkin method. For this purpose, set: (, ) = , [ 1,1]. 2 We define that: () = () () = (, )(, ), (6) () = () = (, ), (7) () = 2 ( ) (, )() = (1 ), (, )(, ). Using ( +1)-point Gauss-Jacobi quadrature formula with weight, to approximate Eqs. (6)-(8) yields: () ():= ,,, ( ), (9) () ():= ,,, (10) () ():= (8),,,,, (11) where { } are the ( + 1)-degree Jacobi-Gauss points associated with,, and { } are the corresponding Jacobi weights. On the other hand, instead of the continuous inner product, the discrete inner product will be implemented by the following equality: JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
4 (, ) =, (12) as a result (, ), =(,), if (Λ). By the definition of,, we have: (, ) =(,, ). (13) The Jacobi pseudo-spectral Galerkin method is to find: () =, (), () =, () (Λ), such that: (, ) =( + +, ), (, ) =(()+, ), (, ) (Λ) (Λ), (14) where and are determined by:,,,,,,,,,, =,,,,, +,,, =(),,. (15) Denoting = [,,...,,,,..., ], Eq. (14) yields a equation of the matrix form: =, (16) where:,,,,,,, 0, 0,,,,, , 0, (, ) =,,,, 0, ,,,,, , , () =,,, 0, (),,, JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
5 3. Some useful lemmas We first introduce some Hilbert spaces. For simplicity, denote () = (/ )(), etc. For a nonnegative integer, define:, ( 1,1):=: (), ( 1,1), 0, with the semi-norm and the norm as:, = (),, = (),, respectively. It is convenient sometime to introduce the semi-norms:,, = (), (,) (). For bounding some approximation error of Jacobi polynomials, we need the following nonuniformly-weighted Sobolev spaces:, ( 1,1):=: (), ( 1,1), 0, equipped with the inner product and the norm as: (, ), =(, ),,, = (, ),. Next, we define the orthogonal projection : (Λ) (Λ) as: (, ) = 0, (Λ), where possesses the following approximation properties ((5.4.11), (5.4.12) and (5.4.24) on p in Ref. ([11]): () (), (17) and:,. (18) We have the following optimal error estimate for the interpolation polynomials based on the Jacobi Gauss points (c.f. [7]). Lemma 3.1 For any function satisfying,, ( 1,1), we have:,, (),, (19) for the Jacobi Gauss points and Jacobi Gauss-Radau points. JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
6 Lemma 3.2 If,, ( 1,1), for some 1 and (Λ), then for the Jacobi Gauss and Jacobi Gauss-Radau integration we have (cf. [7]): (, ), (,),,,,,. (20) We have the following result on the Lebesgue constant for the Lagrange interpolation poly-nomials associated with the zeros of the Jacobi polynomials; (cf. [7]). Lemma 3.3 Let () be the -th Lagrange interpolation polynomials associated with the Gauss, or Gauss-Radau, or Gauss-Lobatto points of the Jacobi polynomials. Then:, := max () = log, 1, 1 [,], = max(, ), otherwise. (21) We now introduce some notation. For 0 and [0,1],, ([ 1,1]) will denote the space of functions whose -th derivatives are Holder continuous with exponent, endowed with the usual norm.,. When = 0,, ([ 1,1]) denotes the space of functions with continuous derivatives on [0, ], also denoted by ([ 1,1]), and with norm.. We will make use of a result of Ragozin ([12-13]), which states that, for each nonnegative Integer and [0,1], there exists a constant, > 0 such that for any function, ([ 1,1]) there exists a polynomial function such that:, (),, (22) where. is the norm of the space ([ 1,1]), and when the function ([ 1,1]). Actually, is a linear operator from, ([ 1,1]) to. We will need the fact that, which be defined by Eq. (11), is compact as an operator from ([0, ]) to, ([ 1,1]) for any 0 <<1 [14]. Lemma 3.4 Let 0 <<1, then, for any function ([ 1,1]), there exists a positive constant such that: ( ) ( ) max (), under the assumption that 0 <<1 for any, [ 1,1] and. This implies that:, (), 0<<1. (23) Clearly, and also satisfy Eq. (24). In our analysis, we shall apply the generalization of Gronwalls lemma. We call such a function () locally integrable on the interval [0, ] if for each [0, ], its Lebesgue integral () is finite. The following result can be found in [15]. Lemma 3.5 Suppose that: () () +() (, )(), [0, ], where, and are locally integrable on the interval [0, ]. Here, all the functions are assumed to be nonnegative. Then: 3812 JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
7 () () +() (, )() (, ) (), [0, ]. Lemma 3.6 Assume that is a nonnegative, locally integrable function defined on [0, ] and satisfying: () () + ( ) (), [0, ], where is a positive constant and () is a nonnegative and continuous function defined on [0, ]. Then, there exists a constant such that: () () +( ) (), [0, ]. Proof. We note that when 0 <<1, the integral ( ) is finite. From Lemma 3.5, we obtain the desired result directly. Lemma 3.7 Assume that is a nonnegative, locally integrable function defined on [ 1, 1] and satisfying: () () + ( ) (), where is a positive constant and () is a nonnegative and continuous function defined on [ 1, 1]. Then, there exists a constant such that: () () + ( ) (). Proof. Let =2 ( 1), =2 ( 1), we obtain that: () () + 2 ( ) (), [0, ]. Using Lemma 3.6 leads to: () () + ( ) (), [0, ]. By the linear transformation =2 ( + 1), =2 ( + 1), desired result follows. Obviously, when = 0, the lemma 3.7 also holds. To prove the error estimate in the weighted -norm, we need the generalized Hardys inequality with weights (see, e.g., [16, 17]). Lemma 3.8 For all measurable function 0, the following generalized Hardys inequality: ()() holds if and only if: () () (), JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
8 Sup () () <, = 1, for the case 1 < <. Here, is an operator of the form: ()() = (, ) (), with (, ) a given kernel,, weight functions, and. We will need the following estimate for the Lagrange interpolation associated with the Jacobi Gaussian collocation points. Lemma 3.9 For every bounded function, there exists a constant independent of such that:, () = () log, 1 <, 1, = max(, ), otherwise, where () is the Lagrange interpolation basis function associated with the Jacobi collocation points. Proof. It is obvious that:, () = () max () [,] max (). [,] By the Lemma 3.3, we obtain the desired result. Lemma 3.10 For every bounded function, there exists a constant independent of such that:, (),, where () is the Lagrange interpolation basis function associated with the Jacobi collocation points. Proof. It is obvious that:, (), =,, = ( ) where =,,,,. As a consequence: sup, (),, =, 3814 JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
9 with =. 4. Convergence for Jacobi pseudo-spectral-galerkin method As, is the interpolation operator which is based on the ( +1)-degree Jacobi-Gauss points with weight,, in terms of Eqs. (13) and (14), the pseudo-spectral Galerkin solution, satisfies: (, ),, +,, =,,,, (, ),,,, =, (),,, (, ) (Λ) (Λ), (24) where: = ( )= (), with: () = = (, ) (, ) ,,, ( ) = (, ),, (, ), (, ),, (, ), (25) in which (, ), represents the continuous inner product with respect to, and (, ) is the corresponding discrete inner product defined by the Gauss-Jacobi quadrature formula. Similar to Eq. (25), we have that: = = (), with: () = = (, ) ,, = ,, (, ), ,, (, ), (26) and: JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
10 = = (), with: () = = ,,,, = , (, ),, (, ) +1 2, (, ),, (, ) (1 ), (, ) (, )., (27) The combination of Eqs. (24)-(27) yields: +, (), +, (),,, =,,,,, (, ), +, (),,, =, (),,, which gives rise to: +, (), +, (), =,, +, (), =, (). (28) By the discussion above, Eqs. (14), (24) and (28) are equivalent. We first consider an auxiliary problem. We want to find, (Λ) such that: (, ) (, ), =(, ), (, ), =((), ), (, ) (Λ) (Λ), (29) where, and are the integral operators defined in Section 2, and (.,.) is still the discrete inner product based on the ( + 1)-degree Jacobi-Gauss points. In terms of the definition of,, Eq. (29) can be written as: (, ),,,, =,,, (, ),, =, (),, (30) which is equivalent to:,, =,,, =, (). (31) When ==0, Eq. (31) can be written as: 3816 JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
11 ,, =0,, =0. In terms of the fact that:,, = +, +,,, = +,. Suppose that: max 2 (,, 2 (), 2. It is clear that from Eq. (6)-(8): = 2 () () + () +, +,, = 2 ( ) (, ) () +,, which yields: ( + ) (( ) +1)( () + () ) + + +, where =,, =,, =,. Using Lemma 3.7 leads to: ( + ) (( ) +1)( + + ) ( + + ), namely: + ( + + ). (32) We now estimate, and. By virtue of Eqs. (22), (23) and Lemma 3.9, we obtain that:, =, =, (1 +, ) log log, 1 < 1, 1 2 < <0. Similarly: JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
12 log,, 1 < 1, 1 2 < <0, and:, log, 1 < 1, 1 < < 0. These, together with Eq. (32), give: log +, 1 < 1 + +, 1 2 < <0, which implies, taking, (0,1 ) such that +>1 2, when is large enough, = = 0. Hence, the and are existent and unique as (Λ) is finite-dimensional. Lemma 4.1. Suppose that, (Λ) and: max 2 (, ), then we have: 2 (), 2, log, +,, 1 < 1 +, +,, 1 2 < <0, (33) +,, +, +log, +,, 1 < 1, 2 ( +, ), +, +,, 1 2 < <0. (34) Proof. Subtracting Eq. (31) from Eq. (3) yields: () +, +, () =,, +, =(), (). (35) Set =(), =(). Direct computation shows that: 3818 JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
13 , +, =, +, ( ) +, +, ( ) =, +( ) ( ), ( ) +, +( ) ( ), ( ) =, ( ) +( ) ( ), ( ) +( ) ( ), ( ) =, + [, ] +,, (36) and:, =, +, ( ) =, +( ) ( ), ( ) = (), ()+( ) ( ), ( ) =, () +, () +,. (37) The insertion of Eqs. (36), (37) into Eq. (35) yields: =, + [, ] +,, =, +,, which implies that: ((. ) +1)( () + () ), (38) where: =,, =,, =,, =,, =,. Using Lemma 3.7 gives: (( ) +1)( ). (39) Then, it follows from Eq. (39) that: +, +, +, +, +,. (40) By using Eq. (18), Lemma 3.9, we obtain that:, =, ( ) log,, 1 < 1 1 +,,, 1 2 < <0, (41) JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
14 and: log,, 1 < 1,,, 1 (42) 2 < <0. We now estimate it is clear that [0, ]. Consequently, using Eqs. (22), (23) and Lemma 3.9 it follows that: =, 1+, log + 1+,, 1 < 1, +, 1 2 < 0, (43) where (0,1 ) and (Λ). Eq. (43) also holds for and. Taking, (0,1 ) such that +>1 2, the estimate Eq. (33) follows from Eqs. (40)-(43), provider that is large enough. Next we prove Eq. (34). Using the generalized Gronwall inequality (Lemma 3.8), we have from Eq. (38) that: +, +, +,, +, + +, + + +,. +, +, +, +, (44) We obtain that from Eqs. (22), (23) and Lemma 3.10: =,, =,,, +, +,, +,. These result, together with the estimates Eqs. (33), (44) and (19), yields (34). Now subtracting Eq. (28) from Eq. (31) leads to:, () +,,, () +,, =0,, (), +, =0, which can be simplified as, by setting =, = :, (), (),, =0,, (), =0. (45) Let = and = be the error corresponding to the Jacobi pseudo-spectral Galerkin solution, of Eq. (14). Now we are prepared to get our global convergence result for problem Eq. (3). Theorem 4.1. Suppose that: 3820 JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
15 max 2 (,, 2 (), 2, and the solution of Eq. (3) is sufficiently smooth. For the Jacobi pseudo spectral Galerkin solution defined in Eq. (14), we have the following estimates: 1) norm of + satisfies: log, +, +log ( + ), 1 < 1, + 2, +, + ( + ), 1. 2 < 0 (46) 2), norm of + satisfies: +, log + +log, +, +, +,, 1 < 1, 2 + +, +, +, +,, 1. 2 < 0 (47) Proof. We first prove the existence and uniqueness of the Jacobi pseudo-spectral Galerkin solution,. As the dimension of (Λ) is finite and Eqs. (14) and (28) are equivalent, we only need to prove that the solution of Eq. (28) is = =0 when = =0. For this purpose, we consider the equation: +, (), +, (), =0, +, (), =0. (48) Obviously Eq. (48) can be written as: =, +,, (), () = + + +, and: =,, () = +, namely: JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
16 = 2 () () + () , = 2 ( ) (, ) () + +. (49) With: =,, =,, =, (), =, (), =,, =, (). Using Eq. (49) gives: ( () + () ), + + ( ) ( () + () ). Namely: ( + ) (( ) +1)( () + () ). (50) Using Lemma 3.7 yields: ( + ) (51) On the other hand, according to Lemma 3.9: =, () (log) 1 (), 1 < (52) 1 (), 2 < <0. By the expression of () in Eq. (25), Lemma 3.2, we have: () +1 2 (, ), (),,,, which, together with Eq. (52), gives: log ( + ), 1 < 1 ( + ), 1 (53) 2 < < JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
17 Similarly, Eq. (53) holds for and. The combination of Eqs. (43) and (53) yields: (log +log ) ( + ), 1 < 1 ( + ) ( + ), 1 (54) 2 < <0. + Based on Eq. (54) and Lemma 3.4 with +>1 2, when is large enough, = =0. As a result, the existence and uniqueness of the Jacobi pseudo-spectral Galerkin solutions, are proved. Now we turn to the error estimate. Actually Eq. (45) can be transformed into: = 2 (() () + ()) +,, +, +, () +, () = 2 ( ) (, )() +, +, (), (55) which yields: + (( ) +1) ( + ) , (56) with =,, =,, =,. Similar to Eq. (52), it follows from Eq. (56) and Lemma 3.7 that: + ( ). (57) Similar to the estimate of Eq. (43), we obtain: log, 1 < < 1, 1. (58) 2 < <0. It also holds for and. In terms of Eqs. (53), (57) and (58), when is large enough, we obtain: log ( + ) log ( ( + ) + ( + ) ), 1 < 1, 2 + ( + ) ( ( + ) + ( + ) ), 1. 2 < <0 (59) By the triangular inequality: JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
18 = + + +, (60) as well as Eqs. (59), (60) and Lemma 4.1, we can obtain the estimated Eq. (46) provided is sufficiently large. Next we prove Eq. (47). Using Lemma 3.7 and the generalized Hardy inequality (Lemma 3.8, == 2), one obtains that from Eq. (56): +, +, + +, +, + +, +, + +, (61) The combination of Eqs. (53), (58) and (59) yields: +, log ( + ) + ( + ), 1 < 1 ( ( + ) + + ), 1 2 < <0. (62) By the triangular inequality again: +, ( + ), + +. (63), In terms of Eqs. (46), (62), (63) and Lemma 4.1, we obtain the desired result. 5. Numerical results We give two numerical examples to confirm our analysis. Example 1. Consider the Volterra integro-differential equation: () = 2() + (1 2) 4 3 (1 + ) + ( ) (). The exact solution is () =. Fig. 1 shows the errors of approximate solution in and Fig. 2 shows the errors weighted, norms obtained by using the Pseudo-spectral methods described above. It is observed that the desired exponential rate of convergence is obtained. Example 2. Consider the Volterra equation integro-differential equation: () = () 16 5 (1 + ) +4(1+) sin() + cos()+ ( ) cos() (). The corresponding exact solution is given by () = cos(). Fig. 3 and Fig. 4 plot the errors for 2 14 in and, norms. Once again the desired spectral accuracy is obtained JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
19 Fig. 1. error of Example 1 Fig. 2., error of Example 1 Fig. 3. error of Example 2 Fig. 4., error of Example 2 6. Concluding remarks This work is concerned with the Jacobi pseudospectral-galerkin methods for solving Volterratype integro-differential equation and the error analysis. To facilitate the use of the methods, we first restate the original integro-differential equation as two simple integral equations of the second kind, then the spectral accuracy associated with and, error estimates are demonstrated theoretically. These results are confirmed by some numerical experiments. We only investigated the case when the solution is smooth in the present work, with the availability of this methodology, it will be possible to extend the results of this paper to the weakly singular VIDEs with nonsmooth solutions which will be the subject of our future work. Acknowledgements This work was financially supported by the National Natural Science Foundation of China ( ), Shanghai Natural Science Foundation (12ZR ), Innovation Program of Shanghai Municipal Education Commission (13ZZ124). References [1] Brunner H. Collocation Methods for Volterra Integral and Related Functional Equations. Cambridge University Press, Cambridge, [2] Chen Y., Tang T. Spectral methods for weakly singular Volterra integral equations with smooth solutions. Journal of Computational and Applied Mathematics, Vol. 233, 2009, p [3] Tao Xia, Xie Ziqing, Zhou Xiaojun Spectral Petrov-Galerkin methods for the second kind Volterra type integro-differential equations. Numerical Mathematics: Theory, Methods and Applications, Vol. 4, Issue 2, 2011, p [4] Ali I., Brunner H., Tang T. A spectral method for pantograph-type delay differential equations and its convergence analysis. Journal of Computational Mathematics, Vol. 27, Issues 2-3, 2009, p JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
20 [5] Brunner H., Schotzau D. Hp-discontinuous Galerkin time-stepping for Volterra integro-differential equations. SIAM Journal on Numerical Analysis, Vol. 44, 2006, p [6] Shen J., Tang T., Wang L.-L. Spectral Methods: Algorithms, Analysis and Applications. Springer Series in Computational Mathematics. Springer, New York, [7] Wang W. Mechanical algorithm for solving the second kind of Volterra integral equation. Applied Mathematics and Computation, Vol. 173, 2006, p [8] Bock I., Lovisek J. On a reliable solution of a Volterra integral equation in a Hilbert space. Applications of Mathematics, Vol. 48, Issue 6, 2003, p [9] Brunner H. The numerical solutions of weakly singular Volterra integral equations by collocation on graded meshes. Mathematics of Computation, Vol. 45, 1985, p [10] Brunner H. Polynomial spline collocation methods for Volterra integro-differential equations with weakly singular kernels. IMA Journal on Numerical Analysis, Vol. 6, 1986, p [11] Canuto C., Hussaini M. Y., Quarteroni A., Zang T. A. Spectral Methods Fundamentals in Single Domains. Springer-Verlag, Berlin, [12] Ragozin D. L. Polynomial approximation on compact manifolds and homogeneous spaces. Transactions of the American Mathematical Society, Vol. 150, 1970, p [13] Ragozin D. L. Constructive polynomial approximation on spheres and projective spaces. Transactions of the American Mathematical Society, Vol. 162, 1971, p [14] Chen Y., Tang T. Convergence analysis of the Jacobi spectral-collocation methods for Volterra integral equations with a weakly singular kernel. Mathemacics of Computation, Vol. 79, 2010, p [15] Willett D. A linear generalization of Gronwall s inequality. Proceedings of the American Mathematical Society, Vol. 16, Issue 4, 1965, p [16] Gogatishvili A., Lang J. The generalized Hardy operator with kernel and variable integral limits in Banach function spaces. Journal of Inequalities and Applications, Vol. 4, Issue 1, 1999, p [17] Kufner A., Persson L. E. Weighted Inequalities of Hardy Type. World Scientific, River Edge, NJ, [18] Atkinson K. E. A Survey of Numerical Methods for the Solution of Fredholm Integral Equations of the Second Kind. SIAM, Philadelphia, [19] Lubich Ch. Fractional linear multi-step methods for Abel-Volterra integral equations of the second kind. Mathematics of Computation, Vol. 45, 1985, p [20] Federson M., Bianconi R., Barbanti L. Linear Volterra integral equations as the limit of discrete systems. Cadernos de Matematica, Vol. 4, 2003, p [21] Shen J., Tang T. Spectral and High-Order Methods with Applications. Science Press, Beijing, [22] Tang T. A note on collocation methods for Volterra integro-differential equations with weakly singular kernels. IMA Journal of Numerical Analysis, Vol. 13, 1993, p [23] Tang T., Xu X., Chen J. On spectral methods for Volterra type integral equations and the convergence analysis. Journal of Computational Mathematics, Vol. 26, 2008, p [24] Jiang Y. J. On spectral methods for Volterra-type integro-differential equations. Journal of Computational and Applied Mathematics, Vol. 230, 2009, p [25] Wan Zhengsu, Chen Yanping, Huang Yunqing Legendre spectral Galerkin method for second-kind Volterra integral equations. Frontiers of Mathematics in China, Vol. 4, Issue 1, 2009, p Xiao-Yong Zhang received Ph.D. degree in Shanghai University, Shanghai, China, in Now he works at College of Arts and Sciences of Shanghai Maritime University. His current research interests include the spectral method, the finite element method. Junlin Li received Ph.D. degree in Beijing University of Aeronautics and Astronautics, Beijing, China, in Now he works at Taiyuan University of Science and Technology. His current research interests include differential equations, automatic control JVE INTERNATIONAL LTD. JOURNAL OF VIBROENGINEERING. DEC 2014, VOLUME 16, ISSUE 8. ISSN
Jacobi 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 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 informationJACOBI SPECTRAL GALERKIN METHODS FOR VOLTERRA INTEGRAL EQUATIONS WITH WEAKLY SINGULAR KERNEL
Bull. Korean Math. Soc. 53 (016), No. 1, pp. 47 6 http://dx.doi.org/10.4134/bkms.016.53.1.47 JACOBI SPECTRAL GALERKIN METHODS FOR VOLTERRA INTEGRAL EQUATIONS WITH WEAKLY SINGULAR KERNEL Yin Yang Abstract.
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 informationConvergence Analysis of Jacobi Pseudo-Spectral Method for the Volterra Delay Integro-Differential Equations
Appl. Math. Inf. Sci. 9, o., 35-45 (5) 35 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/.785/amis/98 Convergence Analysis of Jacobi Pseudo-Spectral Method for the
More informationT u
WANG LI-LIAN Assistant Professor Division of Mathematical Sciences, SPMS Nanyang Technological University Singapore, 637616 T 65-6513-7669 u 65-6316-6984 k lilian@ntu.edu.sg http://www.ntu.edu.sg/home/lilian?
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
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 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 informationhttp://www.springer.com/3-540-30725-7 Erratum Spectral Methods Fundamentals in Single Domains C. Canuto M.Y. Hussaini A. Quarteroni T.A. Zang Springer-Verlag Berlin Heidelberg 2006 Due to a technical error
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 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 informationA Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations
Mathematics A Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations Mohamed Meabed KHADER * and Ahmed Saied HENDY Department of Mathematics, Faculty of Science,
More informationInstructions for Matlab Routines
Instructions for Matlab Routines August, 2011 2 A. Introduction This note is devoted to some instructions to the Matlab routines for the fundamental spectral algorithms presented in the book: Jie Shen,
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationSuperconvergence Results for the Iterated Discrete Legendre Galerkin Method for Hammerstein Integral Equations
Journal of Computer Science & Computational athematics, Volume 5, Issue, December 05 DOI: 0.0967/jcscm.05.0.003 Superconvergence Results for the Iterated Discrete Legendre Galerkin ethod for Hammerstein
More informationA Modification in Successive Approximation Method for Solving Nonlinear Volterra Hammerstein Integral Equations of the Second Kind
Journal of Mathematical Extension Vol. 8, No. 1, (214), 69-86 A Modification in Successive Approximation Method for Solving Nonlinear Volterra Hammerstein Integral Equations of the Second Kind Sh. Javadi
More informationc 2003 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 41, No. 6, pp. 333 349 c 3 Society for Industrial and Applied Mathematics CONVERGENCE ANALYSIS OF SPECTRAL COLLOCATION METHODS FOR A SINGULAR DIFFERENTIAL EQUATION WEIZHANG HUANG,
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 informationBernstein-Szegö Inequalities in Reproducing Kernel Hilbert Spaces ABSTRACT 1. INTRODUCTION
Malaysian Journal of Mathematical Sciences 6(2): 25-36 (202) Bernstein-Szegö Inequalities in Reproducing Kernel Hilbert Spaces Noli N. Reyes and Rosalio G. Artes Institute of Mathematics, University of
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 informationDiscrete Projection Methods for Integral Equations
SUB Gttttingen 7 208 427 244 98 A 5141 Discrete Projection Methods for Integral Equations M.A. Golberg & C.S. Chen TM Computational Mechanics Publications Southampton UK and Boston USA Contents Sources
More informationNUMERICAL SOLUTION FOR CLASS OF ONE DIMENSIONAL PARABOLIC PARTIAL INTEGRO DIFFERENTIAL EQUATIONS VIA LEGENDRE SPECTRAL COLLOCATION METHOD
Journal of Fractional Calculus and Applications, Vol. 5(3S) No., pp. 1-11. (6th. Symp. Frac. Calc. Appl. August, 01). ISSN: 090-5858. http://fcag-egypt.com/journals/jfca/ NUMERICAL SOLUTION FOR CLASS OF
More informationApplied Mathematics Letters. Nonlinear stability of discontinuous Galerkin methods for delay differential equations
Applied Mathematics Letters 23 21 457 461 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Nonlinear stability of discontinuous Galerkin
More informationConvergence of a Gauss Pseudospectral Method for Optimal Control
Convergence of a Gauss Pseudospectral Method for Optimal Control Hongyan Hou William W. Hager Anil V. Rao A convergence theory is presented for approximations of continuous-time optimal control problems
More informationA Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions
Applied Mathematical Sciences, Vol. 5, 2011, no. 23, 1145-1152 A Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions Z. Avazzadeh
More information(This is a sample cover image for this issue. The actual cover is not yet available at this time.)
(This is a sample cover image for this issue. The actual cover is not yet available at this time.) This article appeared in a journal published by Elsevier. The attached copy is furnished to the author
More informationOn the simplest expression of the perturbed Moore Penrose metric generalized inverse
Annals of the University of Bucharest (mathematical series) 4 (LXII) (2013), 433 446 On the simplest expression of the perturbed Moore Penrose metric generalized inverse Jianbing Cao and Yifeng Xue Communicated
More informationarxiv: v1 [math.na] 4 Nov 2015
ON WELL-CONDITIONED SPECTRAL COLLOCATION AND SPECTRAL METHODS BY THE INTEGRAL REFORMULATION KUI DU arxiv:529v [mathna 4 Nov 25 Abstract Well-conditioned spectral collocation and spectral methods have recently
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationOn the decay of elements of inverse triangular Toeplitz matrix
On the decay of elements of inverse triangular Toeplitz matrix Neville Ford, D. V. Savostyanov, N. L. Zamarashkin August 03, 203 arxiv:308.0724v [math.na] 3 Aug 203 Abstract We consider half-infinite triangular
More informationON THE HÖLDER CONTINUITY OF MATRIX FUNCTIONS FOR NORMAL MATRICES
Volume 10 (2009), Issue 4, Article 91, 5 pp. ON THE HÖLDER CONTINUITY O MATRIX UNCTIONS OR NORMAL MATRICES THOMAS P. WIHLER MATHEMATICS INSTITUTE UNIVERSITY O BERN SIDLERSTRASSE 5, CH-3012 BERN SWITZERLAND.
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 informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationGaussian interval quadrature rule for exponential weights
Gaussian interval quadrature rule for exponential weights Aleksandar S. Cvetković, a, Gradimir V. Milovanović b a Department of Mathematics, Faculty of Mechanical Engineering, University of Belgrade, Kraljice
More informationOptimal Spectral-Galerkin Methods Using Generalized Jacobi Polynomials
Journal of Scientific Computing ( 2006) DO: 10.1007/s10915-005-9055-7 Optimal Spectral-Galerkin Methods Using Generalized Jacobi Polynomials Ben-Yu Guo, 1 Jie Shen, 2 and Li-Lian Wang 2,3 Received October
More informationarxiv: v1 [math.na] 24 Feb 2016
Newton type method for nonlinear Fredholm integral equations arxiv:16.7446v1 [math.na] 4 Feb 16 Mona Nabiei, 1 Sohrab Ali Yousefi. 1, Department of Mathematics, Shahid Beheshti University, G. C. P.O. Box
More informationFractional Spectral and Spectral Element Methods
Fractional Calculus, Probability and Non-local Operators: Applications and Recent Developments Nov. 6th - 8th 2013, BCAM, Bilbao, Spain Fractional Spectral and Spectral Element Methods (Based on PhD thesis
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 informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
More informationChebyshev polynomials for solving two dimensional linear and nonlinear integral equations of the second kind
Volume 31, N. 1, pp. 127 142, 2012 Copyright 2012 SBMAC ISSN 0101-8205 / ISSN 1807-0302 (Online) www.scielo.br/cam Chebyshev polynomials for solving two dimensional linear and nonlinear integral equations
More informationSPECTRAL METHODS: ORTHOGONAL POLYNOMIALS
SPECTRAL METHODS: ORTHOGONAL POLYNOMIALS 31 October, 2007 1 INTRODUCTION 2 ORTHOGONAL POLYNOMIALS Properties of Orthogonal Polynomials 3 GAUSS INTEGRATION Gauss- Radau Integration Gauss -Lobatto Integration
More informationVolterra integral equations solved in Fredholm form using Walsh functions
ANZIAM J. 45 (E) ppc269 C282, 24 C269 Volterra integral equations solved in Fredholm form using Walsh functions W. F. Blyth R. L. May P. Widyaningsih (Received 8 August 23; revised 6 Jan 24) Abstract Recently
More informationWojciech Czernous PSEUDOSPECTRAL METHOD FOR SEMILINEAR PARTIAL FUNCTIONAL DIFFERENTIAL EQUATIONS
Opuscula Mathematica Vol. 30 No. 2 2010 http://dx.doi.org/10.7494/opmath.2010.30.2.133 Wojciech Czernous PSEUDOSPECTRAL METHOD FOR SEMILINEAR PARTIAL FUNCTIONAL DIFFERENTIAL EQUATIONS Abstract. We present
More informationThe Approximate Solution of Non Linear Fredholm Weakly Singular Integro-Differential equations by Using Chebyshev polynomials of the First kind
AUSTRALIA JOURAL OF BASIC AD APPLIED SCIECES ISS:1991-8178 EISS: 2309-8414 Journal home page: wwwajbaswebcom The Approximate Solution of on Linear Fredholm Weakly Singular Integro-Differential equations
More informationExistence of solution and solving the integro-differential equations system by the multi-wavelet Petrov-Galerkin method
Int. J. Nonlinear Anal. Appl. 7 (6) No., 7-8 ISSN: 8-68 (electronic) http://dx.doi.org/.75/ijnaa.5.37 Existence of solution and solving the integro-differential equations system by the multi-wavelet Petrov-Galerkin
More informationWe consider in this paper Legendre and Chebyshev approximations of the linear hyperbolic
LEGEDRE AD CHEBYSHEV DUAL-PETROV-GALERKI METHODS FOR HYPERBOLIC EQUATIOS JIE SHE & LI-LIA WAG 2 Abstract. A Legendre and Chebyshev dual-petrov-galerkin method for hyperbolic equations is introduced and
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 informationA Spectral Method for Elliptic Equations: The Neumann Problem
A Spectral Method for Elliptic Equations: The Neumann Problem Kendall Atkinson Departments of Mathematics & Computer Science The University of Iowa Olaf Hansen, David Chien Department of Mathematics California
More informationarxiv:gr-qc/ v1 6 Sep 2006
Introduction to spectral methods Philippe Grandclément Laboratoire Univers et ses Théories, Observatoire de Paris, 5 place J. Janssen, 995 Meudon Cede, France This proceeding is intended to be a first
More informationCOMPOSITION OPERATORS ON HARDY-SOBOLEV SPACES
Indian J. Pure Appl. Math., 46(3): 55-67, June 015 c Indian National Science Academy DOI: 10.1007/s136-015-0115-x COMPOSITION OPERATORS ON HARDY-SOBOLEV SPACES Li He, Guang Fu Cao 1 and Zhong Hua He Department
More informationSPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM. Mirjana Stojanović University of Novi Sad, Yugoslavia
GLASNIK MATEMATIČKI Vol. 37(57(2002, 393 403 SPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM Mirjana Stojanović University of Novi Sad, Yugoslavia Abstract. We introduce piecewise interpolating
More informationNumerical Study of Generalized Three-Dimensional MHD Flow over a Porous Stretching Sheet
Journal of Applied Fluid Mechanics, Vol. 7, No. 3, pp. 473-483, 2014. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. Numerical Study of Generalized Three-Dimensional MHD Flow
More informationA SPECTRAL COLLOCATION METHOD FOR EIGENVALUE PROBLEMS OF COMPACT INTEGRAL OPERATORS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 5, Number 1, Spring 013 A SPECTRAL COLLOCATION METHOD FOR EIGENVALUE PROBLEMS OF COMPACT INTEGRAL OPERATORS CAN HUANG, HAILONG GUO AND ZHIMIN ZHANG
More informationEXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS WITH UNBOUNDED POTENTIAL. 1. Introduction In this article, we consider the Kirchhoff type problem
Electronic Journal of Differential Equations, Vol. 207 (207), No. 84, pp. 2. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS
More informationhal , version 1-22 Nov 2009
Author manuscript, published in "Kinet. Relat. Models 1, 3 8) 355-368" PROPAGATION OF GEVREY REGULARITY FOR SOLUTIONS OF LANDAU EQUATIONS HUA CHEN, WEI-XI LI AND CHAO-JIANG XU Abstract. By using the energy-type
More informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
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 informationA Chance-Constrained Programming Model for Inverse Spanning Tree Problem with Uncertain Edge Weights
A Chance-Constrained Programming Model for Inverse Spanning Tree Problem with Uncertain Edge Weights 1 Xiang Zhang, 2 Qina Wang, 3 Jian Zhou* 1, First Author School of Management, Shanghai University,
More informationCalculus in Gauss Space
Calculus in Gauss Space 1. The Gradient Operator The -dimensional Lebesgue space is the measurable space (E (E )) where E =[0 1) or E = R endowed with the Lebesgue measure, and the calculus of functions
More informationEFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL EQUATIONS BASED ON THE GENERALIZED LAGUERRE POLYNOMIALS
Journal of Fractional Calculus and Applications, Vol. 3. July 212, No.13, pp. 1-14. ISSN: 29-5858. http://www.fcaj.webs.com/ EFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL
More informationAn Approximate Solution for Volterra Integral Equations of the Second Kind in Space with Weight Function
International Journal of Mathematical Analysis Vol. 11 17 no. 18 849-861 HIKARI Ltd www.m-hikari.com https://doi.org/1.1988/ijma.17.771 An Approximate Solution for Volterra Integral Equations of the Second
More informationQuadrature approaches to the solution of two point boundary value problems
Quadrature approaches to the solution of two point boundary value problems Seth F. Oppenheimer Mohsen Razzaghi Department of Mathematics and Statistics Mississippi State University Drawer MA MSU, MS 39763
More informationLocal Fractional Integral Transforms
From the SelectedWorks of Xiao-Jun Yang 2011 Local Fractional Integral Transforms Yang X Available at: https://works.bepress.com/yang_xiaojun/3/ Progress in Nonlinear Science Science is the moving boundary
More informationPartial Differential Equations and the Finite Element Method
Partial Differential Equations and the Finite Element Method Pavel Solin The University of Texas at El Paso Academy of Sciences ofthe Czech Republic iwiley- INTERSCIENCE A JOHN WILEY & SONS, INC, PUBLICATION
More informationIntroduction. Chapter One
Chapter One Introduction The aim of this book is to describe and explain the beautiful mathematical relationships between matrices, moments, orthogonal polynomials, quadrature rules and the Lanczos and
More informationSpectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions
Applied Mathematical Sciences, Vol. 1, 2007, no. 5, 211-218 Spectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions M. Javidi a and A. Golbabai b a Department
More informationApplied Linear Algebra
Applied Linear Algebra Peter J. Olver School of Mathematics University of Minnesota Minneapolis, MN 55455 olver@math.umn.edu http://www.math.umn.edu/ olver Chehrzad Shakiban Department of Mathematics University
More informationBASES FROM EXPONENTS IN LEBESGUE SPACES OF FUNCTIONS WITH VARIABLE SUMMABILITY EXPONENT
Transactions of NAS of Azerbaijan 43 Bilal T. BILALOV, Z.G. GUSEYNOV BASES FROM EXPONENTS IN LEBESGUE SPACES OF FUNCTIONS WITH VARIABLE SUMMABILITY EXPONENT Abstract In the paper we consider basis properties
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationA Class of Shock Wave Solutions of the Periodic Degasperis-Procesi Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.9(2010) No.3,pp.367-373 A Class of Shock Wave Solutions of the Periodic Degasperis-Procesi Equation Caixia Shen
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More informationOn Multivariate Newton Interpolation at Discrete Leja Points
On Multivariate Newton Interpolation at Discrete Leja Points L. Bos 1, S. De Marchi 2, A. Sommariva 2, M. Vianello 2 September 25, 2011 Abstract The basic LU factorization with row pivoting, applied to
More informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
More informationCONVERGENCE OF HYBRID FIXED POINT FOR A PAIR OF NONLINEAR MAPPINGS IN BANACH SPACES
International Journal of Analysis and Applications ISSN 2291-8639 Volume 8, Number 1 2015), 69-78 http://www.etamaths.com CONVERGENCE OF HYBRID FIXED POINT FOR A PAIR OF NONLINEAR MAPPINGS IN BANACH SPACES
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 informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationExistence of Positive Periodic Solutions of Mutualism Systems with Several Delays 1
Advances in Dynamical Systems and Applications. ISSN 973-5321 Volume 1 Number 2 (26), pp. 29 217 c Research India Publications http://www.ripublication.com/adsa.htm Existence of Positive Periodic Solutions
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 35, pp. 2-26, 29. Copyright 29,. ISSN 68-963. ETNA GAUSSIAN DIRECT QUADRATURE METHODS FOR DOUBLE DELAY VOLTERRA INTEGRAL EQUATIONS ANGELAMARIA CARDONE,
More informationSolution of Non Linear Singular Perturbation Equation. Using Hermite Collocation Method
Applied Mathematical Sciences, Vol. 7, 03, no. 09, 5397-5408 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/ams.03.37409 Solution of Non Linear Singular Perturbation Equation Using Hermite Collocation
More informationResearch Article Existence and Uniqueness of Homoclinic Solution for a Class of Nonlinear Second-Order Differential Equations
Applied Mathematics Volume 2012, Article ID 615303, 13 pages doi:10.1155/2012/615303 Research Article Existence and Uniqueness of Homoclinic Solution for a Class of Nonlinear Second-Order Differential
More informationMultiscale method and pseudospectral simulations for linear viscoelastic incompressible flows
Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of
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 informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More informationContents. 2 Sequences and Series Approximation by Rational Numbers Sequences Basics on Sequences...
Contents 1 Real Numbers: The Basics... 1 1.1 Notation... 1 1.2 Natural Numbers... 4 1.3 Integers... 5 1.4 Fractions and Rational Numbers... 10 1.4.1 Introduction... 10 1.4.2 Powers and Radicals of Rational
More informationSingular Value Inequalities for Real and Imaginary Parts of Matrices
Filomat 3:1 16, 63 69 DOI 1.98/FIL16163C Published by Faculty of Sciences Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Singular Value Inequalities for Real Imaginary
More informationA Stable Mixed Finite Element Scheme for the Second Order Elliptic Problems
A Stable Mixed Finite Element Scheme for the Second Order Elliptic Problems Miaochan Zhao, Hongbo Guan, Pei Yin Abstract A stable mixed finite element method (MFEM) for the second order elliptic problems,
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationSolving Poisson s Equations Using Buffered Fourier Spectral Method
Solving Poisson s Equations Using Buffered Fourier Spectral Method Yinlin Dong Hassan Abd Salman Al-Dujaly Chaoqun Liu Technical Report 2015-12 http://www.uta.edu/math/preprint/ Solving Poisson s Equations
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationAn introduction to some aspects of functional analysis
An introduction to some aspects of functional analysis Stephen Semmes Rice University Abstract These informal notes deal with some very basic objects in functional analysis, including norms and seminorms
More informationOn the split equality common fixed point problem for quasi-nonexpansive multi-valued mappings in Banach spaces
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 (06), 5536 5543 Research Article On the split equality common fixed point problem for quasi-nonexpansive multi-valued mappings in Banach spaces
More informationOptimal L p (1 p ) rates of decay to linear diffusion waves for nonlinear evolution equations with ellipticity and dissipation
Nonlinear Analysis ( ) www.elsevier.com/locate/na Optimal L p (1 p ) rates of decay to linear diffusion waves for nonlinear evolution equations with ellipticity and dissipation Renjun Duan a,saipanlin
More informationON THE CONVERGENCE OF THE FOURIER APPROXIMATION FOR EIGENVALUES AND EIGENFUNCTIONS OF DISCONTINUOUS PROBLEMS
SIAM J. UMER. AAL. Vol. 40, o. 6, pp. 54 69 c 003 Society for Industrial and Applied Mathematics O THE COVERGECE OF THE FOURIER APPROXIMATIO FOR EIGEVALUES AD EIGEFUCTIOS OF DISCOTIUOUS PROBLEMS M. S.
More informationBearing fault diagnosis based on TEO and SVM
Bearing fault diagnosis based on TEO and SVM Qingzhu Liu, Yujie Cheng 2 School of Reliability and Systems Engineering, Beihang University, Beijing 9, China Science and Technology on Reliability and Environmental
More informationA MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC PROBLEMS IN SEMI-INFINITE MEDIA
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 A MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC
More informationNotation. General. Notation Description See. Sets, Functions, and Spaces. a b & a b The minimum and the maximum of a and b
Notation General Notation Description See a b & a b The minimum and the maximum of a and b a + & a f S u The non-negative part, a 0, and non-positive part, (a 0) of a R The restriction of the function
More information