Jacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation
|
|
- Allen Mosley
- 5 years ago
- Views:
Transcription
1 Wei et al. SpringerPlus (06) 5:70 DOI 0.86/s z RESEARCH Jacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation Open Access Yunxia Wei *, Yanping Chen, Xiulian Shi 3 and Yuanyuan Zhang 4 *Correspondence: yunxiawei@6.com College of Mathematic and Information Science, Shandong Institute of Business and Technology, Yantai 64005, China Full list of author information is available at the end of the article Abstract We present in this paper the convergence properties of Jacobi spectral collocation method when used to approximate the solution of multidimensional nonlinear Volterra integral equation. The solution is sufficiently smooth while the source function and the kernel function are smooth. We choose the Jacobi Gauss points associated with the multidimensional Jacobi weight function ω(x) = d i= ( x i) α ( + x i ) β, <α, β< d (d denotes the space dimensions) as the collocation points. The error analysis in L -norm and L ω-norm theoretically justifies the exponential convergence of spectral collocation method in multidimensional space. We give two numerical examples in order to illustrate the validity of the proposed Jacobi spectral collocation method. Keywords: Multidimensional nonlinear Volterra integral equation, Jacobi collocation discretization, Multidimensional Gauss quadrature formula, Error estimates Mathematics Subject Classification: 65R0, 45J05, 65N Background We observe that there are many numerical approaches for solving one-dimensional Volterra integral equation, such as Runge Kutta method (Brunner 984; Yuan and Tang 990), polynomial collocation method (Brunner 986; Brunner et al. 00; Brunner and Tang 989), multistep method (Mckee 979; Houwen and Riele 985), hp-discontinuous Galerkin method (Brunner and Schötzau 006) and Taylor series method (Goldfine 977). The spectral collocation method is the most popular form of the spectral methods among practitioners. It is convenient to implement for one-dimensional problems and generally leads to satisfactory results an long as the problems possess sufficient smoothness. In the literature (Tang et al. 008), the authors proposed a Legendre spectral collocation method for Volterra integral equation with a regular kernel in one-dimensional space. Subsequently, Chen and Tang (009, 00), Chen et al. (03), developed the spectral collocation method for one-dimensional weakly singular Volterra integral equation. The proofs of the convergence properties of spectral collocation method for 06 The Author(s). This article is distributed under the terms of the Creative Commons Attribution 4.0 International License ( which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
2 Wei et al. SpringerPlus (06) 5:70 Page of 6 Volterra integro-differential equation with a single spatial variable are given in Wei and Chen (0a, b, 03, 04). Nevertheless, to the best of our knowledge, there have been no works regarding the theoretical analysis of the spectral approximation for multidimensional Volterra integral equation (Atdev and Ashirov 977; Beesack 985; Pachpatte 0; Suryanarayana 97), even for the case with smooth kernel. We shall extend to several space dimensions the approximation results in Tang et al. (008) for a single spatial variable. The expansion of Jacobi will be considered. We will be concerned with Sobolev-type norms that are most frequently applied to the convergence analysis of spectral methods. We get the discrete scheme by using multidimensional Gauss quadrature formula for the integral term. We will provide a rigorous verification of the exponential decay of the errors for approximate solution. We study the multidimensional nonlinear Volterra integral equation of the form y(t, t,..., t d ) + t t 0 0 td 0 K (t, s, t, s,..., t d, s d, y(s, s,..., s d )) ds d...ds ds = g(t, t,..., t d ), t i [0, T i ], i =,,..., d, () by the Jacobi spectral collocation method. Here, g :[0, T ] [0, T ] [0, T d ] R and K : D R R (where D := {(t, s, t, s,..., t d, s d ) : 0 s i t i T i, i =,,..., d}) are given smooth functions. If the given functions are smooth on their respective domains, the solution y is also the smooth function (see Brunner 004). This fact will be the standing point of this paper. Discretization scheme We consider now the domain = (, ) d and we denote an element of R d by x = (x, x,..., x d ). Let <α, β< d, if ω = ω(x) = d i= ( x i) α ( + x i ) β denotes a d-dimensional Jacobi weight function on, we denote by L ω ( ) the space of the measurable functions u : R such that u(x) ω(x)dx < +. It is a Banach space for the norm The space L ω ( ) is a Hilbert space for the inner product (u, v) ω = u(x)v(x)ω(x)dx. L ( ) is the Banach space of the measurable functions u : R that are bounded outside a set of measure zero, equipped with the norm Given a multi-index α = (α, α,..., α d ) of nonnegative integers, we set and u L ω ( ) = ( u(x) ω(x)dx). u L ( ) = ess sup x u(x). α =α + α + +α d
3 Wei et al. SpringerPlus (06) 5:70 Page 3 of 6 D α α v v = α x α x α. d x d We define Hω m( )= {v L ω ( ): for each nonnegative multi-index α with α m, the distributional derivative D α v belongs to L ω ( )}. This is a Hilbert space for the inner product (u, v) m,ω = α m D α u(x)d α v(x)ω(x)dx, which induces the norm v H m ω ( ) = D α v L ω ( ). α m Let { x j,0 j N} denote the Jacobi Gauss points on the one-dimensional interval (, ) (see Canuto et al. 006; Shen and Tang 006). We now consider multidimensional Jacobi interpolation. Let P N ( ) be the space of all algebraic polynomials of degree up to N in each variable x i for i =,,..., d. Let us introduce the Jacobi Gauss points in : x j = ( x j, x j,..., x jd ) for j = (j, j,..., j d ) N d, j = max i d j i N, and denote by I N the interpolation operator at these points, i.e., for each continuous function u, I N u P N satisfies (I N u)( x j ) = u( x j ) for all j N d, j N. We can represent I N u as follows: I N u(x) = u( x j )F j (x), j N where F j (x) = F j (x )F j (x )...F jd (x d ), {F j } N j=0 is the Lagrange interpolation basis function associated with the Jacobi collocation points { x j } N j=0. The multidimensional Jacobi Gauss quadrature formula is f (x)dx j N f ( x j, x j,..., x jd )ω j ω j...ω jd. () We use the variable transformations t i = T i ( + x i), x i [, ] and s i = T i ( + τ i), τ i [, x i ], i =,,..., d to rewrite () as follows u(x, x,..., x d ) + x x xd ˆK (x, τ, x, τ,..., x d, τ d, u(τ, τ,..., τ d ))dτ d...dτ dτ = f (x, x,..., x d ). (3) Here,
4 Wei et al. SpringerPlus (06) 5:70 Page 4 of 6 ( T f (x, x,..., x d ) = g ( + x ), T ( + x ),..., T d ( T Td K T ˆK (x, τ, x, τ,..., x d, τ d, u) = T T ( + τ ),..., T d ( + x d), T d ( + τ d), u and u(x, x,..., x d ) = y( T ( + x ), T ( + x ),..., T d ) ( + x d), ) ( + x d ) is the smooth solution of problem (3). Firstly, Eq. (3) holds at the collocation points x j = ( x j, x j,..., x jd ) on, i.e., ), ( + x ), T ( + τ ), T ( + x ), xj xj xjd u( x j, x j,..., x jd ) + dτ d dτ dτ = f ( x j, x j,..., x jd ). ˆK ( x j, τ, x j, τ,..., x jd, τ d, u(τ, τ,..., τ d )) (4) In order to obtain high order accuracy for the problem (4), we transfer the integral domain [, x j ] [, x j ] [, x jd ] to a fixed interval u( x j, x j,..., x jd ) + K ( x j, τ ( x j, θ ), x j, τ ( x j, θ ),..., x jd, τ d ( x jd, θ d ), u(τ ( x j, θ ), τ ( x j, θ ),..., τ d ( x jd, θ d )))dθ d dθ dθ = f ( x j, x j,..., x jd ), by using the following transformation (5) τ i = τ i ( x ji, θ i ) = + x j i θ i + x j i, i =,,..., d, (6) where K ( x j, τ, x j, τ,..., x jd, τ d, u) = + x j + x j + x j d ˆK ( x j, τ, x j, τ,..., x jd, τ d, u). Next, let u j j j d be the approximation of the function value u( x j ) and use Legendre Gauss quadrature formula, (5) becomes u j j j d + k N K ( x j, τ ( x j, θ k ), x j, τ ( x j, θ k ),..., x jd, τ d ( x jd, θ kd ), u(τ ( x j, θ k ), τ ( x j, θ k ),..., τ d ( x jd, θ kd )))ω k ω k...ω kd = f ( x j, x j,..., x jd ). (7) Here, {θ k, k N} denotes the Legendre Gauss points on the multidimensional space and {ω k, k N} denotes the corresponding weights. Let u N (x, x,..., x d ) = i N u i i...i d F i (x )F i (x )...F id (x d ). Now, we use u N to approximate the solution u. Then, the Jacobi spectral collocation method is to seek u N such that u i i i d satisfy the following collocation equation:
5 Wei et al. SpringerPlus (06) 5:70 Page 5 of 6 u j j j d + K ( x j, τ ( x j, θ k ), x j, τ ( x j, θ k ),..., x jd, τ d ( x jd, θ kd ), i N k N u i i...i d F i (τ ( x j, θ k ))F i (τ ( x j, θ k ))...F id (τ d ( x jd, θ kd )))ω k ω k...ω kd = f ( x j, x j,..., x jd ). (8) We can get the values of u i i i d by solving (8) and obtain the expressions of u N (x) accordingly. Let the error function of the solution be written as e u (x) := u(x) u N (x). Since the exact solution of the problem () can be written as y(t) = u(x) (t i = T i ( + x i), t i [0, T i ], x i [, ]), we can define its approximate solution y N (t) = u N (x). Then the corresponding error function satisfy ( ) ε y (t) := y(t) y N (t) = e u (x) = e u t, t,..., t d. T T T d Remark In our work, we let the multidimensional Jacobi weight function ω(x) = d i= ( x i) α ( + x i ) β, <α, β< d. So ω(x) = ( x)α ( + x) β, <α, β< for d =. In Tang et al. (008), the authors choose α = β = 0. Some lemmas The following result can be found in Canuto et al. (006). Lemma Assume that Gauss quadrature formula is used to integrate the product uφ, where u H m ( ) for some m > d and φ P N ( ). Then there exists a constant C independent of N such that (u, φ) (u, φ) N CN m u H ( ) φ L ( ), (9) where (, ) represents the continuous inner product in L ( ) space and (u, φ) N = u(θ j, θ j,..., θ jd )φ(θ j, θ j,..., θ jd )ω j ω j...ω jd. j N The seminorm is defined as m u H ( ) = k=min(m,n+) i= d k u x k i L ( ). Note that only pure derivatives in each spatial direction appear in this expression. From Fedotov (004), we have the following result on the Lebesgue constant for the Lagrange interpolation polynomials associated with the Jacobi-Gauss points. Lemma Let I N := max x k N F k (x )F k (x ) F kd (x d ), we have
6 Wei et al. SpringerPlus (06) 5:70 Page 6 of 6 O ( (log N) d), if <α, β (, O (N max(α,β)+ ) d), if <α, β< d I N = (, O (N α+ ) d), if <β, <α< d (, O (N β+ ) d), if <α, <β< d. (0) Lemma 3 Assume that u(x) Hω m( ) for m > d and denote (I N u)(x) its interpolation polynomial associated with the multidimensional Jacobi Gauss points { x j, j N}. Then the following estimates hold u I N u L ω ( ) CN m u H ω ( ), () u I N u L ( ) CN d+ m u H ω ( ). () Proof The inequality () can be found in Canuto et al. (006). We now prove (). From Canuto et al. (006), we have u I N u H l ω ( ) CN l m u H ω ( ), 0 l m. We know that H l ω ( ) is embedded in C( ) for l > d, namely, u I N u L ( ) C u I N u H l ω ( ) CN l m u H ω ( ) { CN d+ m u H, when d is an even number, ω ( ) CN d+ m u H, when d is an odd number. ω ( ) CN d+ m u H ω ( ). The following Gronwall Lemma, whose proof can be found in Headley (974), will be essential for establishing our main results. Lemma 4 Suppose M 0, a nonnegative integrable function E(x) satisfies E(x, x,..., x d ) M x x xd + G(x, x,..., x d ), (x, x,..., x d ), where G(x) is also an integrable function, we have E(τ, τ,..., τ d )dτ d...dτ dτ E L ω ( ) C G L ω ( ), E L ( ) C G L ( ). (3) (4) From Theorem in Nevai (984), we have the following mean convergence result of Lagrange interpolation based at the multidimensional Jacobi-Gauss points.
7 Wei et al. SpringerPlus (06) 5:70 Page 7 of 6 Lemma 5 For every bounded function v(x), there exists a constant C independent of v such that sup v( x j )F j (x) C max v(x). (5) N x j N L ω ( ) For r 0 and κ (0, ), C r,κ ( ) will denote the space of functions whose r-th derivatives are Hölder continuous with exponent κ, endowed with the norm: v C r,κ ( ) = max max α v(x) α r x x α xα xα d d α v(x ) x α + max sup xα xα d α v(x ) d x α xα xα d d α r x =x ( (x x ) + (x x ) + +(x d x d )) κ. We shall make use of a result of Ragozin (970, (97) in the following lemma. Lemma 6 For nonnegative integer r and κ (0, ), there exists a constant C r,κ > 0 such that for any function v C r,κ ( ), there exists a polynomial function T N v P N such that v T N v L ( ) C r,κ N (r+κ) v C r,κ ( ), (6) Actually, T N is a linear operator from C r,κ ( ) into P N. Lemma 7 Assume there are constants L 0, L, L,..., L d such that ˆK (x, τ, x, τ,..., x d, τ d, v ) ˆK (x, τ, x, τ,..., x d, τ d, v ) L 0 v v, ˆK xi (x, τ, x, τ,..., x d, τ d, v ) ˆK xi (x, τ, x, τ,..., x d, τ d, v ) L i v v, i =,,..., d. Let M v,v be defined by M v,v (x) = x x xd [ ˆK (x, τ, x, τ,..., x d, τ d, v (τ, τ,..., τ d )) ˆK (x, τ, x, τ,..., x d, τ d, v (τ, τ,..., τ d ))]dτ d...dτ dτ. (7) Then, for any κ (0, ) and v, v C( ), there exists a positive constant C L 0, L, L,..., L d such that M v,v (x ) M v,v (x ) ( (x x ) + (x x ) + +(x d x d )) κ C max v (x) v (x), x (8) for any x, x and x = x. This implies that
8 Wei et al. SpringerPlus (06) 5:70 Page 8 of 6 M v,v C 0,κ ( ) C max v (x) v (x). x (9) Proof For ease of exposition, and without essential loss of generality, we will proof this lemma for d = and assume x < x, x < x, Here, M v,v (x ) M v,v (x ) x = [ ˆK (x, τ, x, τ, v (τ, τ )) ˆK (x, τ, x, τ, v (τ, τ ))]dτ dτ x [ ˆK (x, τ, x, τ, v (τ, τ )) ˆK (x, τ, x, τ, v (τ, τ ))]dτ dτ E + E. E = x [ ˆK (x, τ, x, τ, v (τ, τ )) ˆK (x, τ, x, τ, v (τ, τ ))]dτ dτ [ ˆK (x, τ, x, τ, v (τ, τ )) ˆK (x, τ, x, τ, v (τ, τ ))]dτ dτ (0) x x ˆK (x, τ, x, τ, v (τ, τ )) ˆK (x, τ, x, τ, v (τ, τ )) dτ dτ CL 0 v v L ( )(x x ) CL 0 v v L ( )(x x )+κ (x x ) κ C v v L ( )[(x x ) + (x x ) ] κ, () E = where P = = [ ˆK (x, τ, x, τ, v (τ, τ )) ˆK (x, τ, x, τ, v (τ, τ ))]dτ dτ [ ˆK (x, τ, x, τ, v (τ, τ )) ˆK (x, τ, x (P + P + P 3 )dτ C max(p + P + P 3 ), x [ ˆK (x, τ, x, τ, v (τ, τ )) ˆK (x, τ, x, τ, v (τ, τ ))]dτ [ ˆK (x, τ, x, τ, v (τ, τ )) ˆK (x, τ, x, τ, v (τ, τ ))]dτ dτ, τ, v (τ, τ ))]dτ [ ˆK x (x, τ, ξ, τ, v (τ, τ )) ˆK x (x, τ, ξ, τ, v (τ, τ ))](x x )dτ CL v v L ( )(x x ) C v v L ( )(x x )+κ (x x ) κ C v v L ( )[(x x ) + (x x ) ] κ, ξ (x, x ). () (3)
9 Wei et al. SpringerPlus (06) 5:70 Page 9 of 6 similarly, P = = [ ˆK (x, τ, x, τ, v (τ, τ )) ˆK (x, τ, x, τ, v (τ, τ ))]dτ [ ˆK (x, τ, x, τ, v (τ, τ )) ˆK (x, τ, x, τ, v (τ, τ ))]dτ [ ˆK x (η, τ, x, τ, v (τ, τ )) ˆK x (η, τ, x, τ, v (τ, τ ))](x x )dτ CL v v L ( )(x x ) C v v L ( )[(x x ) + (x x ) ] κ, η (x, x ). (4) P 3 = [ ˆK (x, τ, x, τ, v (τ, τ )) ˆK (x, τ, x, τ, v (τ, τ ))]dτ [ ˆK (x, τ, x, τ, v (τ, τ )) ˆK (x, τ, x, τ, v (τ, τ ))]dτ x L 0 v v dτ C v v L ( )(x x )+κ (x x ) κ C v v L ( )[(x x ) + (x x ) ] κ. The estimate (8) for d = is obtained by combining (0) (4). (5) Error estimates Theorem Let u(x) be the exact solution of the multidimensional nonlinear Volterra integral equation (3), which is smooth. u N (x) is the approximate solution, i.e., u(x) u N (x). Assume that k θ k i K (x, θ, x, θ,..., x d, θ d, v ) k Then there is a constant C such that the errors satisfy for m > d +, θ k i L ik v v, i =,,..., d; k =,,..., m, L = max L ik. i d, k m K (x, θ, x, θ,..., x d, θ d, v ) u u N L ( ) CN m (log N) d K + N d+ u H ω ( ), if <α, β ( ), dk N max(α,β)+ + N d+ u H ω ( ), if <α, β< d ( ), dk N α+ + N d+ u H ω ( ), if <β, <α< d ( ), dk N β+ + N d+ u H ω ( ), if <α, <β< d. (6)
10 Wei et al. SpringerPlus (06) 5:70 Page 0 of 6 where K = max j N K ( x j, θ, x j, θ,..., x jd, θ d, u(θ, θ,..., θ d )) H ( ), C d, L, L 0, L, L,..., L d. Proof We subtract (8) from (5) to get the error equation u( x j, x j,..., x jd ) u j j...j d + where Using the variable transformation (6), we have [ K ( x j, τ ( x j, θ ), x j, τ ( x j, θ ),..., x jd, τ d ( x jd, θ d ), u(τ ( x j, θ ), τ ( x j, θ ),..., τ d ( x jd, θ d ))) K ( x j, τ ( x j, θ ), x j, τ ( x j, θ ),..., x jd, τ d ( x jd, θ d ), u N (τ ( x j, θ ), τ ( x j, θ ),..., τ d ( x jd, θ d )))]dθ d dθ dθ = I( x j, x j,..., x jd ), I( x j, x j,..., x jd ) = k N K ( x j, τ ( x j, θ k ), x j, τ ( x j, θ k ),..., x jd, τ d ( x jd, θ kd ), u N (τ ( x j, θ k ), τ ( x j, θ k ),..., τ d ( x jd, θ kd )))ω k ω k...ω kd K ( x j, τ ( x j, θ ), x j, τ ( x j, θ ),..., x jd, τ d ( x jd, θ d ), u N (τ ( x j, θ ), τ ( x j, θ ),..., τ d ( x jd, θ d )))dθ d...dθ dθ. u( x j, x j,..., x jd ) u j j...j d + xj xj... xjd [ ˆK ( x j, τ, x j, τ,..., x jd, τ d, u(τ, τ,..., τ d )) ˆK ( x j, τ, x j, τ,..., x jd, τ d, u N (τ, τ,..., τ d ))]dτ d...dτ dτ = I( x j, x j,..., x jd ). (7) Multiplying F j (x )F j (x )...F jd (x d ) on both sides of Eq. (7) and summing up j N yield x x e u (x, x,..., x d ) + Consequently, xd [ ˆK (x, τ, x, τ,..., x d, τ d, u(τ, τ,..., τ d )) ˆK (x, τ, x, τ,..., x d, τ d, u N (τ, τ,..., τ d ))]dτ d...dτ dτ = J (x, x,..., x d ) + J (x, x,..., x d ) + J 3 (x, x,..., x d ). (8) x x xd e u (x, x,..., x d ) L 0... e u (τ, τ,..., τ d ) dτ d...dτ dτ + J (x, x,..., x d ) + J (x, x,..., x d ) + J 3 (x, x,..., x d ), (9) where
11 Wei et al. SpringerPlus (06) 5:70 Page of 6 J (x) = I( x j, x j,..., x jd )F j (x )F j (x )...F jd (x d ), j N J (x) = u(x, x,..., x d ) (I N u)(x, x,..., x d ), J 3 (x) = x x xd [ ˆK (x, τ, x, τ,..., x d, τ d, u(τ, τ,..., τ d )) ˆK (x, τ, x, τ,..., x d, τ d, u N (τ, τ,..., τ d ))]dτ d...dτ dτ I N x x xd [ ˆK (x, τ, x, τ,..., x d, τ d, u(τ, τ,..., τ d )) ˆK (x, τ, x, τ,..., x d, τ d, u N (τ, τ,..., τ d ))]dτ d...dτ dτ. It follows from the Gronwall inequality in Lemma 4 that e u L ( ) C ( J L ( ) + J L ( ) + J 3 L ( )). Using (9) and (0), we have ( ) J L ( ) C I N max I( x j, x j,..., x jd ) j N C I N N m K ( x j, θ, x j, θ,..., x jd, θ d, u N (θ, θ,..., θ d )) H ( ), C I N N m (K + K ( x j, θ, x j, θ,..., x jd, θ d, u N (θ, θ,..., θ d )) K ( x j, θ, x j, θ,..., x jd, θ d, u(θ, θ,..., θ d )) H ( ) ). (30) (3) A straightforward computation shows that K ( x j, θ, x j, θ,..., x jd, θ d, u N (θ, θ,..., θ d )) K ( x j, θ, x j, θ,..., x jd, θ d, u(θ, θ,..., θ d )) H ( ) ( m d k K θi k ( x j, θ, x j, θ,..., x jd, θ d, u N (θ, θ,..., θ d )) k= i= k θ k i K ( x j, θ, x j, θ,..., x jd, θ d, u(θ, θ,..., θ d )) L ( ) ) ( m k= i= d L ik u N u L ( ) ) L m d u N u L ( ) C e u L ( ). (3) Due to Lemma 3, J L ( ) CN d+ m u H ω ( ). (33) By virtue of Lemmas 6 and 7,
12 Wei et al. SpringerPlus (06) 5:70 Page of 6 J 3 L ( ) = (I I N )M u,un L ( ) = (I I N )(M u,un T N M u,un ) L ( ) ( + I N ) M u,un T N M u,un L ( ) C I N N κ M u,un C 0,κ ( ) C I N N κ e u L ( ) C(log N) d N κ e u L ( ), if <α, β, C(N max(α,β)+ ) d N κ e u L ( ), if <α, β< d, C(N α+ ) d N κ e u L ( ), if <β, <α< d, C(N β+ ) d N κ e u L ( ), if <α, <β< d. (34) We now obtain the estimate for e u L ( ) by using (30) (34), e u L ( ) CN m( ) I N K + N d+ u H, ω ( ) where in last step we have used the following assumption, 0 <κ<, if <α, β, (max(α, β) + )d <κ<, if <α, β< d, (α + )d <κ<, if <β, <α< d, (β + )d <κ<, if <α, <β< d. This completes the proof of the theorem. (35) Theorem If the hypotheses given in Theorem hold and κ satisfies (35), then Proof m u u N L ω ( ) CN ( + N κ (log N) d )K + ( + N d+ κ ) u H ω ( ), ( if <α, β, + N κ( ) ) N max(α,β)+ d K + ( + N d+ κ ) u H ω ( ), if <α, β< d, ( + N κ( N α+ ) d) K + ( + N d+ κ ) u H ω ( ), if <β, <α< d, ( + N κ( N β+ ) d ) K + ( + N d+ κ ) u H ω ( ), if <α, <β< d. By using (8) and Gronwall inequality in Lemma 4, we obtain that ) e u L ω ( ) ( J C L ω ( ) + J L ω ( ) + J 3 L ω ( ). (36) (37) Using Lemmas, 5 and (3) we have for J L ω ( ) C max I(x) CN m (K + e u L x ω ( ) ). (38) Due to Lemma 3,
13 Wei et al. SpringerPlus (06) 5:70 Page 3 of 6 J L ω ( ) CN m u H ω ( ). (39) By virtue of Lemmas 6 and 7, J 3 L ω ( ) = (I I N )M u,un L ω ( ) = (I I N )(M u,un T N M u,un ) L ω ( ) M u,un T N M u,un L ω ( ) + I N (M u,un T N M u,un ) L ω ( ) C M u,un T N M u,un L ( ) CN κ e u L ( ) CN m κ( ) I N K + N d+ u H. ω ( ) (40) The desired estimate (36) is obtained by combining (37) (40) and using the same technique as in the proof of Theorem. Numerical results We give two numerical examples to confirm our analysis. To examine the accuracy of the results, L ω and L errors are employed to assess the efficiency of the method. All the calculations are supported by the software Matlab. Example We consider the following two-dimensional Volterra integral equation x y u(x, y) + cos(x + y)e ξη u(ξ, η)dηdξ = e xy + cos(x + y)(x + )(y + ). (4) The corresponding exact solution is given by u(x, y) = e xy. We select α = 3, β =. Table shows the errors u u N L ω ( ) and u u N L ( ) obtained by using the spectral collocation method described above. Furthermore, the numerical results are plotted for N in Fig.. It is observed that the desired exponential rate of convergence is obtained. Example Consider the equation with Table The errors u u N L ω ( ) and u u N L ( ) N 4 6 L -error 9.373e e e 008 L ω-error.85e e e 009 N 8 0 L -error e e e 04 L ω-error e 03.30e 03.34e 03
14 Wei et al. SpringerPlus (06) 5:70 Page 4 of L error L ω error errors in log scale N Fig. The errors u u N versus the number of collocation points in L and L ω norms x y v(x, y) + cos(x + ξ)v(ξ, η)dηdξ = sin(x + y) 4 sin(3x + y) + 4 sin(x + y ) (x + )cos(x y) + (x + )cos(x + ) + 4 sin(3x ) sin(x 3). 4 (4) The corresponding exact solution is given by v(x, y) = sin(x + y). We select α = 3, β = 3 4. Table shows the errors v v N L ω ( ) and v v N L ( ). The numerical results are plotted for N in Fig.. Conclusions In this paper, we proposed a spectral collocation method based on Jacobi orthogonal polynomials to obtain approximate solution for multidimensional nonlinear Volterra integral equation. The most important contribution of this work is that we are able to Table The errors v v N L ω ( ) and v v N L ( ) N 4 6 L -error.0746e e e 006 L ω-error 4.797e e e 006 N 8 0 L -error.503e e 0.99e 03 L ω-error e e e 04
15 Wei et al. SpringerPlus (06) 5:70 Page 5 of 6 0 L error L ω error 0 4 errors in log scale N Fig. The errors v v N versus the number of collocation points in L and L ω norms demonstrate rigorously that the errors of spectral approximations decay exponentially in both L ( ) norm and L ω ( ) norm on d-dimensional space, which is a desired feature for a spectral method. Authors contributions YW and YC carried out the spectral collocation method studies, performed the error analysis and drafted the manuscript. XS participated in the numerical experiments. YZ helped to draft the manuscript. All authors read and approved the final manuscript. Author details College of Mathematic and Information Science, Shandong Institute of Business and Technology, Yantai 64005, China. School of Mathematical Sciences, South China Normal University, Guangzhou 5063, China. 3 School of Mathematics and Statistics, Zhaoqing University, Zhaoqing 5606, China. 4 Department of Mathematics and Information Science, Yantai University, Yantai 64005, China. Acknowledgements This work is supported by National Natural Science Foundation of China (40347, , 745, 64055, 4693). Competing interests The authors declare that they have no competing interests. Received: 6 June 06 Accepted: September 06 References Atdev S, Ashirov S (977) Solutions of multidimensional nonlinear Volterra operator equations. Ukr Math J 9: Beesack PR (985) Systems of multidimensional Volterra integral equations and inequalities. Nonlinear Anal Theory 9: Brunner H (984) Implicit Runge Kutta methods of optimal order for Volterra integro-differential equations. Math Comput 4:95 09 Brunner H (986) Polynomial spline collocation methods for Volterra integro-differential equations with weakly singular kernels. IMA J Numer Anal 6: 39 Brunner H (004) Collocation methods for Volterra integral and related functional equations. Cambridge University Press, Cambridge Brunner H, Tang T (989) Polynomial spline collocation methods for the nonlinear Basset equation. Comput Math Appl 8: Brunner H, Schötzau D (006) hp-discontinuous Galerkin time-stepping for Volterra integrodifferential equations. SIAM J Numer Anal 44:4 45 Brunner H, Pedas A, Vainikko G (00) Piecewise polynomial collocation methods for linear Volterra integro-differential equations with weakly singular kernels. SIAM J Numer Anal 39: Canuto C, Hussaini MY, Quarteroni A et al (006) Spectral methods fundamentals in single domains. Springer, Berlin
16 Wei et al. SpringerPlus (06) 5:70 Page 6 of 6 Chen Y, Tang T (009) Spectral methods for weakly singular Volterra integral equations with smooth solutions. J Comput Appl Math 33: Chen Y, Tang T (00) Convergence analysis of the Jacobi spectral-collocation methods for Volterra integral equations with a weakly singular kernel. Math Comput 79:47 67 Chen Y, Li X, Tang T (03) A note on Jacobi spectral-collocation methods for weakly singular Volterra integral equations with smooth solutions. J Comput Math 3:47 56 Fedotov AI (004) Lebesgue constant estimation in multidimensional Sobolev space. J Math 4:5 3 Goldfine A (977) Taylor series methods for the solution of Volterra integral and integro-differential equations. Math Comput 3: Headley VB (974) A multidimensional nonlinear Gronwall inequality. J Math Anal Appl 47:50 55 Mckee S (979) Cyclic multistep methods for solving Volterra integro-differential equations. SIAM J Numer Anal 6:06 4 Nevai P (984) Mean convergence of Lagrange interpolation. Trans Am Math Soc 8: Pachpatte BG (0) Multidimensional integral equations and inequalities. Springer, Berlin 9 Ragozin DL (970) Polynomial approximation on compact manifolds and homogeneous spaces. Trans Am Math Soc 50:4 53 Ragozin DL (97) Constructive polynomial approximation on spheres and projective spaces. Trans Am Math Soc 6:57 70 Shen J, Tang T (006) Spectral and high-order methods with applications. Science Press, Beijing Suryanarayana MB (97) On multidimensional integral equations of Volterra type. Pac J Math 4: Tang T, Xu X, Chen J (008) On spectral methods for Volterra integral equations and the convergence analysis. J Comput Math 6: van der Houwen PJ, te Riele HJJ (985) Linear multistep methods for Volterra integral and integro-differential equations. Math Comput 45: Wei YX, Chen Y (0a) Convergence analysis of the spectral methods for weakly singular Volterra integro-differential equations with smooth solutions. Adv Appl Math Mech 4: 0 Wei YX, Chen Y (0b) Legendre spectral collocation methods for pantograph Volterra delay-integro-differential equation. J Sci Comput 53: Wei YX, Chen Y (03) A spectral method for neutral Volterra integro-differential equation with weakly singular kernel. Numer Math Theory Methods Appl 6: Wei YX, Chen Y (04) Legendre spectral collocation method for neutral and high-order Volterra integro-differential equation. Appl Numer Math 8:5 9 Yuan W, Tang T (990) The numerical analysis of implicit Runge Kutta methods for a certain nonlinear integro-differential equation. Math Comput 54:55 68
1462. 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 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 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 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 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 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 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 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 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 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 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 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 informationNew hybrid conjugate gradient methods with the generalized Wolfe line search
Xu and Kong SpringerPlus (016)5:881 DOI 10.1186/s40064-016-5-9 METHODOLOGY New hybrid conjugate gradient methods with the generalized Wolfe line search Open Access Xiao Xu * and Fan yu Kong *Correspondence:
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 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 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 informationSobolev Spaces 27 PART II. Review of Sobolev Spaces
Sobolev Spaces 27 PART II Review of Sobolev Spaces Sobolev Spaces 28 SOBOLEV SPACES WEAK DERIVATIVES I Given R d, define a multi index α as an ordered collection of integers α = (α 1,...,α d ), such that
More informationSolving a class of nonlinear two-dimensional Volterra integral equations by using two-dimensional triangular orthogonal functions.
Journal of Mathematical Modeling Vol 1, No 1, 213, pp 28-4 JMM Solving a class of nonlinear two-dimensional Volterra integral equations by using two-dimensional triangular orthogonal functions Farshid
More informationResearch Article A Legendre Wavelet Spectral Collocation Method for Solving Oscillatory Initial Value Problems
Applied Mathematics Volume 013, Article ID 591636, 5 pages http://dx.doi.org/10.1155/013/591636 Research Article A Legendre Wavelet Spectral Collocation Method for Solving Oscillatory Initial Value Problems
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 informationL p Bounds for the parabolic singular integral operator
RESEARCH L p Bounds for the parabolic singular integral operator Yanping Chen *, Feixing Wang 2 and Wei Yu 3 Open Access * Correspondence: yanpingch@26. com Department of Applied Mathematics, School of
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 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 informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationMultiplesolutionsofap-Laplacian model involving a fractional derivative
Liu et al. Advances in Difference Equations 213, 213:126 R E S E A R C H Open Access Multiplesolutionsofap-Laplacian model involving a fractional derivative Xiping Liu 1*,MeiJia 1* and Weigao Ge 2 * Correspondence:
More informationResearch Article Existence of Periodic Positive Solutions for Abstract Difference Equations
Discrete Dynamics in Nature and Society Volume 2011, Article ID 870164, 7 pages doi:10.1155/2011/870164 Research Article Existence of Periodic Positive Solutions for Abstract Difference Equations Shugui
More informationVariational Theory of Solitons for a Higher Order Generalized Camassa-Holm Equation
International Journal of Mathematical Analysis Vol. 11, 2017, no. 21, 1007-1018 HIKAI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2017.710141 Variational Theory of Solitons for a Higher Order Generalized
More informationOn Solutions of Evolution Equations with Proportional Time Delay
On Solutions of Evolution Equations with Proportional Time Delay Weijiu Liu and John C. Clements Department of Mathematics and Statistics Dalhousie University Halifax, Nova Scotia, B3H 3J5, Canada Fax:
More informationThe Generalized Viscosity Implicit Rules of Asymptotically Nonexpansive Mappings in Hilbert Spaces
Applied Mathematical Sciences, Vol. 11, 2017, no. 12, 549-560 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.718 The Generalized Viscosity Implicit Rules of Asymptotically Nonexpansive
More informationApproximation algorithms for nonnegative polynomial optimization problems over unit spheres
Front. Math. China 2017, 12(6): 1409 1426 https://doi.org/10.1007/s11464-017-0644-1 Approximation algorithms for nonnegative polynomial optimization problems over unit spheres Xinzhen ZHANG 1, Guanglu
More informationResearch Article Solvability of a Class of Integral Inclusions
Abstract and Applied Analysis Volume 212, Article ID 21327, 12 pages doi:1.1155/212/21327 Research Article Solvability of a Class of Integral Inclusions Ying Chen and Shihuang Hong Institute of Applied
More informationA LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.
A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE THOMAS CHEN AND NATAŠA PAVLOVIĆ Abstract. We prove a Beale-Kato-Majda criterion
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 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 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 informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationApplication of the perturbation iteration method to boundary layer type problems
DOI 10.1186/s40064-016-1859-4 RESEARCH Open Access Application of the perturbation iteration method to boundary layer type problems Mehmet Pakdemirli * *Correspondence: mpak@cbu.edu.tr Applied Mathematics
More informationPermutation transformations of tensors with an application
DOI 10.1186/s40064-016-3720-1 RESEARCH Open Access Permutation transformations of tensors with an application Yao Tang Li *, Zheng Bo Li, Qi Long Liu and Qiong Liu *Correspondence: liyaotang@ynu.edu.cn
More informationRunge-Kutta Method for Solving Uncertain Differential Equations
Yang and Shen Journal of Uncertainty Analysis and Applications 215) 3:17 DOI 1.1186/s4467-15-38-4 RESEARCH Runge-Kutta Method for Solving Uncertain Differential Equations Xiangfeng Yang * and Yuanyuan
More informationCubic spline collocation for a class of weakly singular Fredholm integral equations and corresponding eigenvalue problem
Cubic spline collocation for a class of weakly singular Fredholm integral equations and corresponding eigenvalue problem ARVET PEDAS Institute of Mathematics University of Tartu J. Liivi 2, 549 Tartu ESTOIA
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 informationResearch Article A New Fractional Integral Inequality with Singularity and Its Application
Abstract and Applied Analysis Volume 212, Article ID 93798, 12 pages doi:1.1155/212/93798 Research Article A New Fractional Integral Inequality with Singularity and Its Application Qiong-Xiang Kong 1 and
More informationResearch Article On a Quasi-Neutral Approximation to the Incompressible Euler Equations
Applied Mathematics Volume 2012, Article ID 957185, 8 pages doi:10.1155/2012/957185 Research Article On a Quasi-Neutral Approximation to the Incompressible Euler Equations Jianwei Yang and Zhitao Zhuang
More informationTHE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 4, Number 1, Winter 1992 THE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS J.-P. KAUTHEN ABSTRACT. We present a method of lines
More informationGlobal Solutions for a Nonlinear Wave Equation with the p-laplacian Operator
Global Solutions for a Nonlinear Wave Equation with the p-laplacian Operator Hongjun Gao Institute of Applied Physics and Computational Mathematics 188 Beijing, China To Fu Ma Departamento de Matemática
More informationG-SPLINE INTERPOLATION FOR APPROXIMATING THE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS USING LINEAR MULTI- STEP METHODS
Journal of Al-Nahrain University Vol.0(), December, 00, pp.8- Science G-SPLINE INTERPOLATION FOR APPROXIMATING THE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS USING LINEAR MULTI- STEP METHODS Osama H.
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationSpectral collocation method for compact integral operators
Wayne State University Wayne State University Dissertations --0 Spectral collocation method for compact integral operators Can Huang Wayne State University, Follow this and additional works at: http://digitalcommons.wayne.edu/oa_dissertations
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 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 informationResearch Article Hermite Wavelet Method for Fractional Delay Differential Equations
Difference Equations, Article ID 359093, 8 pages http://dx.doi.org/0.55/04/359093 Research Article Hermite Wavelet Method for Fractional Delay Differential Equations Umer Saeed and Mujeeb ur Rehman School
More informationA Study on Linear and Nonlinear Stiff Problems. Using Single-Term Haar Wavelet Series Technique
Int. Journal of Math. Analysis, Vol. 7, 3, no. 53, 65-636 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.988/ijma.3.3894 A Study on Linear and Nonlinear Stiff Problems Using Single-Term Haar Wavelet Series
More informationAn Efficient Multiscale Runge-Kutta Galerkin Method for Generalized Burgers-Huxley Equation
Applied Mathematical Sciences, Vol. 11, 2017, no. 30, 1467-1479 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.7141 An Efficient Multiscale Runge-Kutta Galerkin Method for Generalized Burgers-Huxley
More informationOn the existence of solution to a boundary value problem of fractional differential equation on the infinite interval
Shen et al. Boundary Value Problems 5 5:4 DOI.86/s366-5-59-z R E S E A R C H Open Access On the existence of solution to a boundary value problem of fractional differential equation on the infinite interval
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 informationA NEW COLLOCATION METHOD FOR SOLVING CERTAIN HADAMARD FINITE-PART INTEGRAL EQUATION
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 6, Number, Pages 4 54 c 9 Institute for Scientific Computing and Information A NEW COLLOCATION METHOD FOR SOLVING CERTAIN HADAMARD FINITE-PART
More informationAdaptation of Taylor s Formula for Solving System of Differential Equations
Nonlinear Analysis and Differential Equations, Vol. 4, 2016, no. 2, 95-107 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/nade.2016.51144 Adaptation of Taylor s Formula for Solving System of Differential
More informationWELL-POSEDNESS OF THE TWO-DIMENSIONAL GENERALIZED BENJAMIN-BONA-MAHONY EQUATION ON THE UPPER HALF PLANE
WELL-POSEDNESS OF THE TWO-DIMENSIONAL GENERALIZED BENJAMIN-BONA-MAHONY EQUATION ON THE UPPER HALF PLANE YING-CHIEH LIN, C. H. ARTHUR CHENG, JOHN M. HONG, JIAHONG WU, AND JUAN-MING YUAN Abstract. This paper
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationResearch Article Bounds of Solutions of Integrodifferential Equations
Abstract and Applied Analysis Volume 211, Article ID 571795, 7 pages doi:1.1155/211/571795 Research Article Bounds of Solutions of Integrodifferential Equations Zdeněk Šmarda Department of Mathematics,
More informationStrong convergence theorems for total quasi-ϕasymptotically
RESEARCH Open Access Strong convergence theorems for total quasi-ϕasymptotically nonexpansive multi-valued mappings in Banach spaces Jinfang Tang 1 and Shih-sen Chang 2* * Correspondence: changss@yahoo.
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
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 informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
More informationResearch Article The Existence of Countably Many Positive Solutions for Nonlinear nth-order Three-Point Boundary Value Problems
Hindawi Publishing Corporation Boundary Value Problems Volume 9, Article ID 575, 8 pages doi:.55/9/575 Research Article The Existence of Countably Many Positive Solutions for Nonlinear nth-order Three-Point
More informationComputing the Determinant and Inverse of the Complex Fibonacci Hermitian Toeplitz Matrix
British Journal of Mathematics & Computer Science 9(6: -6 206; Article nobjmcs30398 ISSN: 223-085 SCIENCEDOMAIN international wwwsciencedomainorg Computing the Determinant and Inverse of the Complex Fibonacci
More informationOn boundary value problems for fractional integro-differential equations in Banach spaces
Malaya J. Mat. 3425 54 553 On boundary value problems for fractional integro-differential equations in Banach spaces Sabri T. M. Thabet a, and Machindra B. Dhakne b a,b Department of Mathematics, Dr. Babasaheb
More informationExistence and Uniqueness of Anti-Periodic Solutions for Nonlinear Higher-Order Differential Equations with Two Deviating Arguments
Advances in Dynamical Systems and Applications ISSN 973-531, Volume 7, Number 1, pp. 19 143 (1) http://campus.mst.edu/adsa Existence and Uniqueness of Anti-Periodic Solutions for Nonlinear Higher-Order
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 informationASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS. Tian Ma. Shouhong Wang
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 11, Number 1, July 004 pp. 189 04 ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS Tian Ma Department of
More informationApproximate Solution of BVPs for 4th-Order IDEs by Using RKHS Method
Applied Mathematical Sciences, Vol. 6, 01, no. 50, 453-464 Approximate Solution of BVPs for 4th-Order IDEs by Using RKHS Method Mohammed Al-Smadi Mathematics Department, Al-Qassim University, Saudi Arabia
More informationA collocation method for solving some integral equations in distributions
A collocation method for solving some integral equations in distributions Sapto W. Indratno Department of Mathematics Kansas State University, Manhattan, KS 66506-2602, USA sapto@math.ksu.edu A G Ramm
More informationSome new sequences that converge to the Ioachimescu constant
You et al. Journal of Inequalities and Applications (06) 06:48 DOI 0.86/s3660-06-089-x R E S E A R C H Open Access Some new sequences that converge to the Ioachimescu constant Xu You *, Di-Rong Chen,3
More informationL ϕ -Embedding inequalities for some operators on differential forms
Shi et al. Journal of Inequalities and Applications 2015 2015:392 DOI 10.1186/s13660-015-0887-x R S A R C H Open Access L -mbedding inequalities for some operators on differential forms Guannan Shi 1,
More informationA general theorem on the stability of a class of functional equations including quadratic additive functional equations
Lee and Jung SpringerPlus 20165:159 DOI 10.1186/s40064-016-1771-y RESEARCH A general theorem on the stability of a class of functional equations including quadratic additive functional equations Yang Hi
More informationHankel Operators plus Orthogonal Polynomials. Yield Combinatorial Identities
Hanel Operators plus Orthogonal Polynomials Yield Combinatorial Identities E. A. Herman, Grinnell College Abstract: A Hanel operator H [h i+j ] can be factored as H MM, where M maps a space of L functions
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
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 informationBulletin of the. Iranian Mathematical Society
ISSN: 1017-060X (Print) ISSN: 1735-8515 (Online) Bulletin of the Iranian Mathematical Society Vol. 42 (2016), No. 1, pp. 129 141. Title: On nonlocal elliptic system of p-kirchhoff-type in Author(s): L.
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 informationNontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation
Advances in Dynamical Systems and Applications ISSN 973-532, Volume 6, Number 2, pp. 24 254 (2 http://campus.mst.edu/adsa Nontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
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 NONLINEAR DIFFERENTIAL EQUATION INVOLVING REFLECTION OF THE ARGUMENT
ARCHIVUM MATHEMATICUM (BRNO) Tomus 40 (2004), 63 68 A NONLINEAR DIFFERENTIAL EQUATION INVOLVING REFLECTION OF THE ARGUMENT T. F. MA, E. S. MIRANDA AND M. B. DE SOUZA CORTES Abstract. We study the nonlinear
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationAn interface problem with singular perturbation on a subinterval
Xie Boundary Value Problems 2014, 2014:201 R E S E A R C H Open Access An interface problem with singular perturbation on a subinterval Feng Xie * * Correspondence: fxie@dhu.edu.cn Department of Applied
More informationAppendix A Functional Analysis
Appendix A Functional Analysis A.1 Metric Spaces, Banach Spaces, and Hilbert Spaces Definition A.1. Metric space. Let X be a set. A map d : X X R is called metric on X if for all x,y,z X it is i) d(x,y)
More informationADOMIAN-TAU OPERATIONAL METHOD FOR SOLVING NON-LINEAR FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS WITH PADE APPROXIMANT. A. Khani
Acta Universitatis Apulensis ISSN: 1582-5329 No 38/214 pp 11-22 ADOMIAN-TAU OPERATIONAL METHOD FOR SOLVING NON-LINEAR FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS WITH PADE APPROXIMANT A Khani Abstract In this
More informationSolving Nonlinear Two-Dimensional Volterra Integral Equations of the First-kind Using the Bivariate Shifted Legendre Functions
International Journal of Mathematical Modelling & Computations Vol. 5, No. 3, Summer 215, 219-23 Solving Nonlinear Two-Dimensional Volterra Integral Equations of the First-kind Using the Bivariate Shifted
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 informationAdvances in Difference Equations 2012, 2012:7
Advances in Difference Equations This Provisional PDF corresponds to the article as it appeared upon acceptance. Fully formatted PDF and full text (HTML) versions will be made available soon. Extremal
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationGLOBAL EXISTENCE OF SOLUTIONS TO A HYPERBOLIC-PARABOLIC SYSTEM
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 35, Number 4, April 7, Pages 7 7 S -99396)8773-9 Article electronically published on September 8, 6 GLOBAL EXISTENCE OF SOLUTIONS TO A HYPERBOLIC-PARABOLIC
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 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 informationON PERIODIC SOLUTIONS TO SOME LAGRANGIAN SYSTEM WITH TWO DEGREES OF FREEDOM
ON PERIODIC SOLUTIONS TO SOME LAGRANGIAN SYSTEM WITH TWO DEGREES OF FREEDOM OLEG ZUBELEVICH DEPT. OF THEORETICAL MECHANICS, MECHANICS AND MATHEMATICS FACULTY, M. V. LOMONOSOV MOSCOW STATE UNIVERSITY RUSSIA,
More informationPower Series Solutions We use power series to solve second order differential equations
Objectives Power Series Solutions We use power series to solve second order differential equations We use power series expansions to find solutions to second order, linear, variable coefficient equations
More information