JACOBI SPECTRAL GALERKIN METHODS FOR VOLTERRA INTEGRAL EQUATIONS WITH WEAKLY SINGULAR KERNEL
|
|
- Laura Hubbard
- 5 years ago
- Views:
Transcription
1 Bull. Korean Math. Soc. 53 (016), No. 1, pp JACOBI SPECTRAL GALERKIN METHODS FOR VOLTERRA INTEGRAL EQUATIONS WITH WEAKLY SINGULAR KERNEL Yin Yang Abstract. We propose and analyze spectral and pseudo-spectral Jacobi- Galerkin approaches for weakly singular Volterra integral equations (VIEs). We provide a rigorous error analysis for spectral and pseudospectral Jacobi-Galerkin methods, which show that the errors of the approximate solution decay exponentially in L norm and weighted L - norm. The numerical examples are given to illustrate the theoretical results. 1. Introduction Volterra integral equations (VIEs) arise widely in mathematical models of certain biological and physical phenomena. Due to the wide application of these equations, they must be solved successfully with efficient numerical methods. In this article, we are concerned with the numerical study of the following Volterra integral equations with weakly singular kernel: (1) y(t) = t 0 (t τ) γ K(t,τ)y(τ)dτ +f(t), 0 < γ < 1, t [0,T], where the source function f and the kernel function K are given and the function y(t) is the unknown function. Here, f and K are assumed to be sufficiently smooth on their respective domains [0,T] and 0 τ t T. The numerical treatment of the VIEs (1) is a challenge work, mainly due to the fact that the solutions of (1) usually have a weak singularity at t = 0, even when the inhomogeneous term f(t) is regular. Collocation spectral methods and the corresponding error analysis have been provided recently [14, 16] for (1) without the singular kernel (i.e., γ = 0) in case of the underlying Received February 1, 015; Revised June 7, Mathematics Subject Classification. Primary 65R0, 45J05, 65N1. Key words and phrases. spectral Galerkin methods, Jacobi polynomial, Volterra integral equations with weakly singular kernels. This work was financially supported by NSFC Project ( , , ), and Project of Scientific Research Fund of Hunan Provincial Science and Technology Department (013RS4057). 47 c 016 Korean Mathematical Society
2 48 Y. YANG solutions are smooth. Chen and Tang [4, 17] developed a novel spectral Jacobicollocation method to solve (1) and provided a rigorous error analysis which theoretically justifies the spectral rate of convergence, see also [19] for Jacobi spectral-collocation method for fractional integro-differential equations. In[0], the authors extended the Legendre-collocation methods to nonlinear Volterra integral equations. Recently, in [16], the authors provided a Legendre spectral Galerkin method for second-kind Volterra integral equations, [15, 18] provide general spectral and pseudo-spectral Jacobi-Petrov-Galerkin approaches for the second kind Volterra integro-differential equations. Inspired by those work, we extend spectral Galerkin approach to Volterra integral equations with weakly singular kernels and provide a rigorous convergence analysis for the spectral and pseudo-spectral Jacobi-Galerkin methods, which indicates that the proposed methods converges exponentially provided that the data in the given VIEs are smooth. This paper is organized as follows. In Section, we demonstrate the implementation of the spectral and pseudo-spectral Galerkin approaches for Volterra integral equations with weakly singular kernels. Some lemmas useful for establishing the convergence result will be provided in Section 3. The convergence analysis for both spectral and pseudo-spectral Jacobi-Galerkin methods will be given in Section 4. Numerical results will be carried out in Section 5, which will be used to verify the theoretical result obtained in Section 4. Finally, in Section 6, we end with conclusion and future work.. Spectral and pseudo-spectral Galerkin methods For the sake of applying the theory of orthogonal polynomials, we use the change of variable t = 1 t T(1+x), x = T, τ = 1 τ T(1+s), s = T, and let u(x) = y( 1 T(1+x)), g(x) = f(1 T(1+x)), k(x,s) = ( ) 1 γ T K( 1 T(1+x)1 T(1+s)). Then we get the following one dimension VIEs with weakly singular kernel defined on [,1] () u(x) = x (x s) γ k(x,s)u(s)ds+g(x), 0 < γ < 1, x Λ := [,1]. We formulate the Jacobi-spectral Galerkin schemes and investigate the global convergence properties for the problem (). For this purpose, we first define a
3 JACOBI SPECTRAL GALERKIN METHODS FOR VIES 49 linear integral operator G : C(Λ) C(Λ) by (Gφ)(x) := x (x s) γ k(x,s)φ(s)ds. Then, the problem () reads: find u = u(x) such that (3) u(x) = (Gu)(x)+g(x), x Λ and its weak form is to find u L (Λ) such that (4) (u,v) = (Gu,v)+(g,v), v L (Λ), where (, ) denotes the usual inner product in the L -space. Let us demonstrate the numerical implementation of the spectral Jacobi- Galerkin approach first. Denote by N the set of all nonnegative integers. For any N N, P N denotes the set ofall algebraicpolynomials ofdegree at most N in Λ, φ j (x) is the j-th Jacobi polynomial corresponding to the weight function (x) = (1 x) α (1+x) β. As a result, P N = spanφ 0 (x),φ 1 (x),...,φ N (x)}. Our spectral Jacobi-Galerkin approximation of (3) is now defined as: Find u N P N such that (5) (u N,v N ) = (Gu N,v N ) +(g,v N ), v P N, where (u,v) = 1 u(x)v(x) (x)dx is the continuous inner product. Set u N (x) = N j=0 ξ jφ j (x). Substituting it into (5) and taking v N = φ i (x), we obtain (6) ξ j (φ i,φ j ) = j=0 ξ j (φ i,gφ j ) + (φ i,g), j=0 which leads to an equation of the matrix form (7) where Aξ = Bξ +C, j=0 ξ = [ξ 0,ξ 1,...,ξ N ] T, A i,j = (φ i,φ j ), B i,j = (φ i,gφ j ), C i = (φ i,g). Next we propose the pseudo-spectral Jacobi-Galerkin method. Set s(x,θ) = 1+x θ + x, θ 1,
4 50 Y. YANG it is clear that (8) with Gu(x) = = x 1 (x s) γ k(x,s)u(s)ds (1 θ) γ k(x,s(x,θ))u(s(x,θ))dθ ( ) 1 γ 1+x k(x,s(x,θ)) = k(x,s(x,θ)). Using (N + 1)-point Gauss quadrature formula to approximate (8) yields (9) Gu(x) G N u(x) := k=0 k(x,s(x,θ k ))u(s(x,θ k )) γ,0 k, where θ k } N k=0 are the (N + 1)-degree Jacobi-Gauss points corresponding to the weights γ,0 k } N k=0. On the other hand, instead of the continuous inner product, the discrete inner product will be implemented in (5) and (6), i.e., (10) (u,v) (u,v),n = As a result, k=0 u(x k )v(x k ) k (x k ). (u,v) = (u,v),n, if u,v P N. Substitute (9) and (10) into (5), the pseudo-spectral Jacobi-Galerkin method is to find ū N (x) = ξ j φ j (x) such that (11) j=0 (ū N,v N ),N = (G N ū N,v N ),N +(g,v N ),N, v N P N, ξ j } N j=0 are determined by (1) ξ j (φ i,φ j ),N = ξ j (φ i,g N φ j ),N +(φ i,g),n, j=0 the matrix form (13) where j=0 Ā ξ = B ξ + C, ξ = [ ξ 0, ξ 1,..., ξ N ] T, Ā i,j = (φ i,φ j ),N, B i,j = (φ i,gφ j ),N, Ci = (φ i,g),n.
5 JACOBI SPECTRAL GALERKIN METHODS FOR VIES Some useful lemmas In this section, we will provide some elementary lemmas, which are important for the derivation of the main results in the subsequent section. First we define the Jacobi orthogonal projection operator Π N : L P N which satisfies (14) (Π N u,v N ) = (u,v N ), u L,v N P N, Further, define equipped with the norm L (Λ) = u : u is measurable and u α,β < }, ( 1 u = u (x) (x)dx ) 1 H m (Λ) = u : Dk u L α,β(λ),0 k m}, u H m = ( m d k u dx k k=0 )1. Lemma 3.1 (see [15]). Suppose that u H m (Λ) and m 1. Then (15) u Π N u CN m u H m;n, (16) u Π N u CN 3 4 m u H m;n, where u H m;n denotes the seminorm defined by α,β(λ) m u H m;n = d k u dx k k=min(+1) Lemma 3. (see [6]). Suppose that u L (Λ). Then Π N u C u, Π N u C u. Lemma 3.3 (see []). Assume that an (N + 1)-point Gauss quadrature formula relative to the Jacobi weight is used to integrate the product uϕ, where u H m (Λ) with I for some m 1 and ϕ P N. Then there exists a constant C independent of N such that 1 (17) u(x)ϕ(x) (x)dx (u,ϕ),n CN m u H ϕ.. 1/.
6 5 Y. YANG Lemma 3.4 (see []). Assume that u H m (Λ), m 1 I α,β α,β N u denotes the interpolation operator of u based on (N + 1)-degree Jacobi-Gauss points corresponding to the weight function (x) with < α,β < 1. Then (18a) (18b) u I α,β N u CN m u H u I α,β N u CN 1 m logn u H c, α,β < 1, CN ν+1 m u H c, otherwise, ν = max(α,β),, where c =, denotes the Chebyshev weight function. Lemma 3.5 (see [9]). For every bounded function u, there exists a constant C, independent of u such that I α,β N u(x j) C u, where I α,β N u(x) = N j=0 u(x j)f j (x) is the Lagrange interpolation basis function associated with (N + 1)-degree Jabobi-Gauss points corresponding to the weight function (x). Lemma 3.6 (see [9]). Assume that F j (x)} N j=0 are the N-th degree Lagrange basis polynomials associated with the Gauss points of the Jacobi polynomials. Then I α,β N max F j (x) (19) = x [,1] j=0 O(logN), < α,β 1, O(N ν+1 ), ν = max(α,β), otherwise. Lemma 3.7 (see [11, 1]). For a nonnegative integer r and κ (0,1), there exists a constant C r,κ > 0 such that for any function v C r,κ ([,1]), there exists a polynomial function T N v P N such that (0) v T N v C r,κ N (r+κ) v r,κ, where r,κ is the standard norm in C r,κ ([,1]), T N is a linear operator from C r,κ ([,1]) into P N, as stated in [11, 1]. Lemma 3.8 (see [5]). Let κ (0,1) and let G be defined by (Gv)(x) = x (x τ) µ K(x,τ)v(τ)dτ. Then, for any function v C([,1]), there exists a positive constant C such that Gv(x ) Gv(x ) x x C max x [,1] v(x),
7 JACOBI SPECTRAL GALERKIN METHODS FOR VIES 53 under the assumption that 0 < κ < 1 µ, for any x,x [,1] and x x. This implies that Gv 0,κ C max v(x), 0 < κ < 1 µ. x [,1] Lemma 3.9 ([7], Gronwall inequality). Suppose L 0 and G(x) is a nonnegative, locally integrable function defined on [, 1] satisfying E(x) G(x)+L Then there exists a constant C such that x E(τ)dτ. E L p (I) L G L p (I), p 1. Here and below, C denotes a positive constant which is independent of N, and whose particular meaning will become clear by the context in which it arises. 4. Convergence analysis for spectral and pseudo-spectral Jacobi-Galerkin method According to (5) and the definition of the projection operator Π N, the spectral Jacobi-Galerkin solution u N satisfies (1) u N = Π N Gu N +Π N g. Theorem 4.1. Suppose that u N is the spectral Jacobi-Galerkin solution determined by (5), if the solution u of () satisfies u H (Λ), then we have the following error estimate () u u N CN 3 4 m u H, u u N CN m( ) CN 3 4 κ +1 Proof. Subtracting (1) from (3), yields u H, κ (0,γ). (3) u u N = Gu Π N Gu N +g Π N g. Set e = u u N, direct computation shows that (4) Gu Π N Gu N = Gu Π N Gu+Π N G(u u N ) = Gu Π N Gu+G(u u N ) [G(u u N ) Π N G(u u N )] = (u g) Π N (u g)+g(u u N ) [G(u u N ) Π N G(u u N )] = (u g) Π N (u g)+ge [Ge Π N Ge].
8 54 Y. YANG The insertion of (4) into (3) yields (5) e(x) = = x x (x s) γ k(x,s)e(s)ds+u Π N u+[π N Ge Ge] (x s) γ k(x,s)e(s)ds+i 1 +I, where I 1 = u Π N u, I = Π N Ge Ge. It follows from the Gronwall inequality that (6) e(x) C I i. By Lemma 3.1, (7) i=1 I 1 CN 3 4 m u H m;n, In the virtue of Lemma 3., Lemma 3.7 and Lemma 3.8, (8) I = (Π N I)Ge = (Π N I)(Ge T N Ge) Π N (Ge T N Ge) + Ge T N Ge C Ge T N Ge CN κ Ge 0,κ, κ (0,1 µ) = (0,γ) CN κ e. Combing (6), (7) and (8), we obtain, when N is large enough, u u N CN 3 4 m u H m;n. Now we investigate the -error estimate. It follows from (5) and Gronwall inequality Lemma 3.9 that (9) e(x) C I i. Due to Lemma 3.1, (30) i=1 I 1 CN m u H. It follows from Lemma 3., Lemma 3.7 and Lemma 3.8 that I = (Π N I)Ge (31) = (Π N I)(Ge T N Ge) Π N (Ge T N Ge) + Ge T N Ge C Ge T N Ge CN κ Ge 0,κ
9 JACOBI SPECTRAL GALERKIN METHODS FOR VIES 55 CN κ e, κ (0,1 µ) = (0,γ). The combination of (9), (30) and (31) yields, u u N CN m( 1+N 3 4 κ) u H, provided N is large enough. Hence, the theorem is proved. As I α,β N is the interpolation operator which is based on the (N +1)-degree Jacobi-Guass points, in terms of (11), the pseudo-spectral Galerkin solution ū N satisfies (3) (ū N,v N ) = (I α,β N G Nū N,v N ) +(I α,β N g,v N). Let (33) I(x) = Gū N G N ū N = 1 (1 θ) γ k(x,s(x,θ))un (s(x,θ))dθ k=0 k(x,s(x,θ k ))u N (s(x,θ k )) γ,0 k. Combing (3) and (33), yields (34) (ū N,v N ) =(I α,β N Gū N I α,β N I(x),v N) +(I α,β N g,v N), v N P N, which gives rise to (35) ū N = I α,β N Gū N I α,β N I(x)+Iα,β N g. We first consider an auxiliary problem, i.e., we want to find û N P N, such that (36) (û N,v N ),N = (Gû N,v N ),N +(g,v N ),N, v N P N. In terms of the definition of I α,β N, (36) can be written as (37) (û N,v N ) = (I α,β N Gû N,v N ) +(I α,β N g,v N), v N P N, which is equivalent to (38) û N = I α,β N Gû N +I α,β N g. Lemma 4.. Suppose û N is determined by (38), ν = max(α,β) min(0, 1 γ) and 0 < κ < 1 γ, if the solution u of () satisfies u H(Λ),
10 56 Y. YANG we have (39) 1 CN m logn u H, ν < 1 u û N c, CN ν+1 m u H, 1 c ν < min(0, 1 γ). ) CN ( u m H +lognn 1 κ u c H m;n, ν < 1 u û N, ) CN ( u m H +N ν+1 κ u c H m;n, 1 ν < min(0,γ 1 ). Proof. Subtracting (38) from (3), yields (40) u û N = Gu I α,β N Gu N +g I α,β N g. Set ε = u û N, direct computation shows that (41) Gu I α,β N Gû N = Gu I α,β N Gu+(Gu Gû N) = (u g) I α,β N (u g)+gε [Gε Iα,β N Gε] = I α,β N u u+iα,β N g g +Gε [Gε Iα,β N Gε]. The insertion of (41) into (40) yields (4) ε = x [ ] (Gu Gû N ) I α,β N (Gu Gû N) (x s) γ k(x,s)ε(s)ds+j 1 +J, where J 1 = u I α,β N u, J = I α,β N Gε Gε. It follows from Gronwall inequality that (43) Due to Lemma 3.4, (44) J 1 = u I α,β N u ε(x) C( J 1 + J ). 1 CN m logn u H, ν < 1 c, CN ν+1 m u H, otherwise. c By virtue of Lemma 3.7, Lamme 3.8 and Lemma 3.6, we have (45) J = (I α,β N I)Gε = (I α,β N I)(Gε T NGε) I α,β N (Gε T NGε) + Gε T N Gε O(logN) Gε T N Gε, < ν 1, O(N ν+1 ) Gε T N Gε, otherwise,
11 JACOBI SPECTRAL GALERKIN METHODS FOR VIES 57 ClogNN κ Gε 0,κ, < ν 1, CN ν+1 κ Gε 0,κ, 1 ν < min(0, 1 0 < κ < 1 γ, γ), ClogNN κ ε, < ν 1, CN ν+1 κ ε, 1 ν < min(0, 1 0 < κ < 1 γ. γ), Combing (43), (44) and (45) we obtain, when N is large enough, 1 CN m logn u H, ν < 1 u û N c, CN ν+1 m u H, 1 c ν < min(0, 1 γ). Now we investigate the -error estimate. It follows from (4) and the Gronwall inequality that (46) By Lemma 3.4, (47) ε(x) C( J 1 + J ). J 1 = u I α,β N u CN m u H m;n It follows from Lemma 3.5, Lemma 3.7 and Lemma 3.8 that (48) J = (I α,β N I)Gε = (I α,β N I)(Gε T NGε). I α,β N (Gε T NGε) + Gε T N Gε C Gε T N Gε CN κ ε, 1 < κ < 1 γ. Combing (46), (47) and (48) we obtain, when N is large enough, ) CN ( u m H +lognn 1 κ u c H m;n, ν < 1 u û N, ) CN ( u m H +N ν+1 κ u c H m;n, 1 ν < min(0,γ 1 ). This completes the proof of the lemma. Theorem 4.3. Suppose that the solution u of () satisfies u H (Λ), ν = max(α,β) min(0, 1 γ) and 0 < κ < 1 γ, for the pseudo spectral Jacobi-Galerkin solution ū N, such that (11) holds, we have (49) u ū N [ ] CN m logn N 1 u H +K u, ν < 1 c ], CN ν+1 [N m u 1 H +K u, 1 c ν min(0, 1 γ), u ū N
12 58 Y. YANG CN [(1+logNN m κ )K u + u H c ν [( < 1, CN m 1+N ν+1 )K κ u + u H c 1 ν min(0, 1 γ), where K = max x (,1) k(x,s(x, )) H m;n γ,0. Proof. Now subtracting (35) from (38) leads to ] +lognn 1 κ u H m;n, ] +N ν+1 κ u H m;n, (50) ū N û N = I α,β N (Gū N Gû N ) I α,β N I(x), which can be simplified as, by setting E = ū N û N (51) E = I α,β N GE Iα,β N I(x) = GE GE +I α,β N GE Iα,β N I(x) = GE +Q 1 I α,β N I(x) with Q 1 = I α,β N GE GE. It follows from the Gronwall inequality that (5) Similarly to (45), we have E = Q 1 + I α,β N I(x). Q 1 = I α,β N GE GE CN κ logn E, µ < 1, CN ν+1 κ E, 1 ν < min(0, 1 γ), with 0 < κ < 1 γ. Using Lemma 3.3 and Lemma 3.6, we have ClogN max I(x), ν < I α,β N I(x) x (,1), CN ν+1 max I(x), x (,1) ν < 0, CN m logn max x (,1) k(x,s(x, )) H m;n ū N γ,0, ν < 1 γ,0, CN m N ν+1 max k(x,s(x, )) H x (,1) m;n ū N γ,0, 1 γ,0 ν < 0, CN m logn max k(x,s(x, )) x (,1) H m;n ū N, ν < 1 γ,0, CN m N ν+1 max k(x,s(x, )) x (,1) H m;n ū N, 1 γ,0 ν < 0, CN m logn max k(x,s(x, )) x (,1) H m;n ( u + E ), ν < 1 γ,0, CN m N ν+1 max k(x,s(x, )) x (,1) H m;n ( u + E ), 1 γ,0 ν < 0.
13 JACOBI SPECTRAL GALERKIN METHODS FOR VIES 59 Set K = max x (,1) k(x,s(x, )) H m;n we now obtain the estimate E by γ,0(λ), using (5) CN m lognk u, ν < 1 (53) E, CN ν+1 m K u, 1 ν < min(0, 1 γ). Next, we will give the error estimates in. It follows from (51) and the Gronwall inequality that (54) E = Q 1 + I α,β N I(x). Q 1 can be established in a similar way as (48), Q 1 = I α,β N GE GE CN κ E, 0 < κ < 1 γ. Using Lemma 3.3 and Lemma 3.5, we have I α,β N I(x) C I(x) CN m K û N γ,0 By the convergence result in (53) (m = 1), we have for sufficiently large N. So that E u CN m K ( u + E ). I α,β N I(x) CN m K u. We obtain, when N is large enough, (55) E CN m (1+logNN κ )K u, ν < 1, CN m ( 1+N ν+1 κ )K u, 1 ν < min(0, 1 γ). Finally, it follows from triangular inequality, u ū N u û N + û N ū N, u ū N u û N + û N ū N, as well as Lemma 4., (53) and (55), we can obtain the desired estimated (49). 5. Numerical experiments We give numerical examples to confirm our analysis. To examine the accuracyofthe results, and errorsareemployedtoassessthe efficiency of the method. All the calculations are supported by the software Matlab. Example 5.1. We consider the following the linear weakly singular Volterra integral equation (1) with y(t) = f(t) t 0 (t τ) γ y(τ)dτ, 0 t 6, 0 < γ < 1,
14 60 Y. YANG L Error L α, β Error N 16 Figure 1. Example 5.1: Comparison between approximate solution of pseudo-spectral Jacobi-Galerkin method and exact solution of y(t) = t γ sin(t) with γ = 0.8 (left). The errors versus the number of collocation points in L and weighted L norms (right). with f(t) = sin(t) t γ +π 1 1 t γ Γ(1 γ)b( 1 γ, t )sin(t ), where B(, ) is the Bessel function defined by ( ) k ( ) θ ϑ ϑ 4 B(θ,ϑ) = k!γ(θ +k +1), k=0 The corresponding exact solution is given by y(t) = t γ sin(t). We have reported the obtained numerical results of pseudo-spectral Jacobi- Galerkin for N = 16 and γ = 0.8 in Figure 1 (left). We can see that the numerical result of pseudo-spectral Jacobi-Galerkin approximation solutions are in good agreement with exact solution y(t) = t γ sin(t). Figure 1 (right) illustrates L and weighted L errors of pseudo-spectral Jacobi-Galerkin method versus the number N of the steps. Clearly, these figures show the exponential rate of convergence predicted by the proposed method. In practice, many Volterra equations are usually nonlinear. However, the nonlinearity adds rather little to the difficulty of obtaining a numerical solution. The methods described above remain applicable. Below we will provide a numerical example using the spectral technique proposed in this work. Example 5.. Consider the following the nonlinear weakly singular Volterra integral equation (1) with y(t) = f(t) t 0 (t τ) γ tan(y(t))dτ, 0 t, 0 < γ < 1,
15 JACOBI SPECTRAL GALERKIN METHODS FOR VIES L Error L Error N 16 Figure. Example 5.: Comparison between approximate solution of pseudo-spectral Jacobi-Galerkin method and exact solution of y(t) = atan(t 1 γ ) with γ = 0.4 (left). The errors versus the number of collocation points in L and weighted L norms (right). with f(t) = t γ arctan(t 1 γ )+ πt γ ( t ) γ Γ(1 γ) Γ( 3 γ). This example has a smooth solution y(t) = arctan(t 1 γ ). This is a nonlinear problem. The numerical scheme leads to a nonlinear system, and a proper solver for the nonlinear system (e.g., Newton method) should be used. Figure (left) presents the approximate and exact solution with γ = 0.4, which are found in excellent agreement. Next, Figure (right) illustrates the L and weighted L errors of the pseudo-spectral-galerkin method. These results indicate that the spectral accuracy is obtained for this problem. 6. Conclusions and future work This work has been concerned with the spectral and pseudo-spectral Jacobi- Galerkin analysis of the Volterra integral equations with weakly singular kernel. The most important contribution of this work is that we are able to demonstrate rigorously that the errors of spectral approximations decay exponentially in both infinity and weighted norms, which is a desired feature for a spectral method. Although in this work our convergence theory does not cover the nonlinear case, the methods described above remain applicable, it will be possible to extend the results of this paper to nonlinear case which will be the subject of our future work. Acknowledgment. The author is grateful to the referees for the useful discussion and valuable comments that enabled them to improve some details of the paper.
16 6 Y. YANG References [1] P. Baratella and A. Orsi, A new approach to the numerical solution of weakly singular Volterra integral equations, J. Comput. Appl. Math. 163 (004), no., [] C. Canuto, M. Y. Hussaini, and A. Quarteroni, Spectral Methods, Fundamentals in single domains, Sci. Comput., Springer-Verlag, Berlin, 006 [3] Y. Chen and T. Tang, Spectral methods for weakly singular Volterra integral equations with smooth solutions, J. Comput. Appl. Math. 33 (009), no. 4, [4], Convergence analysis of the Jacobi spectral-collocation methods for Volterra integral equation with a weakly singular kernel, Math. Comp. 79 (010), no. 69, [5] D. Colton and R. Kress, Inverse Acoustic and Electromagnetic Scattering Theory, Appl. Math. Sci., Springer-Verlag, Heidelberg, nd Edition, [6] J. Douglas, T. Dupont, and L. Wahlbin, The stability in L q of the L -projection into finite element function spaces, Numer. Math. 3 (1974), no. 3, [7] D. Henry, Geometric Theory of Semilinear Parabolic Equations, Springer-Verlag, [8] A. Kufner and L. E. Persson, Weighted Inequalities of Hardy Type, World Scientific, New York, 003. [9] G. Mastroianni and D. Occorsto, Optimal systems of nodes for Lagrange interpolation on bounded intervals: a survey, J. Comput. Appl. Math. 134 (001), no. 1-, [10] P. Nevai, Mean convergence of Lagrange interpolation: III, Trans. Amer. Math. Soc. 8 (1984), no., [11] D. L. Ragozin, Polynomial approximation on compact manifolds and homogeneous spaces, Trans. Amer. Math. Soc. 150 (1970), [1], Constructive polynomial approximation on spheres and projective spaces, Trans. Amer. Math. Soc. 16 (1971), [13] S. Jie, T. Tao, and L. Wang, Spectral methods. Algorithms, analysis and applications, Springer Ser. Comput. Math., 41. Springer, Heidelberg, 011. [14] T. Tang, X. Xu, and J. Cheng. On Spectral methods for Volterra integral equation and the convergence analysis, J. Comput. Math. 6 (008), no. 6, [15] X. Tao, Z. Xie, and X. Zhou, Spectral Petrov-Galerkin methods for the second kind Volterra type integro-differential equations, Numer. Math. Theory Methods Appl. 4 (011), no., [16] Z. Wan, Y. Chen, and Y. Huang, Legendre spectral Galerkin method for second-kind Volterra integral equations, Front. Math. China 4 (009), no. 1, [17] Y. Wei and Y. Chen, Convergence analysis of the spectral methods for weakly singular Volterra integro-differential equations with smooth solutions, Adv. Appl. Math. Mech. 4 (01), no. 1, 1 0. [18] Z. Xie and X. Li, and T. Tang, Convergence analysis of spectral Galerkin methods for Volterra type integral equations, J. Sci. Comput. 53 (01), no., [19] Y. Yang, Y. Chen, and Y. Huang, Convergence analysis of the Jacobi spectral-collocation method for fractional integro-differential equations, Acta Math. Sci. Ser. B Engl. Ed. 34 (014), no. 3, [0], Convergence analysis of Legendre-collocation methods for nonlinear Volterra type integro equations, Adv. Appl. Math. Mech. 7 (015), no. 1, Hunan Key Laboratory for Computation and Simulation in Science and Engineering School of Mathematics and Computational Science Xiangtan University Xiangtan , P. R. China address: yangyinxtu@xtu.edu.cn
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 informationJacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation
Wei et al. SpringerPlus (06) 5:70 DOI 0.86/s40064-06-3358-z RESEARCH Jacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation Open Access
More 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 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 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 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 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 informationNON-EXTINCTION OF SOLUTIONS TO A FAST DIFFUSION SYSTEM WITH NONLOCAL SOURCES
Electronic Journal of Differential Equations, Vol. 2016 (2016, No. 45, pp. 1 5. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu NON-EXTINCTION OF
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 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 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 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 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 information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationSEMILINEAR ELLIPTIC EQUATIONS WITH DEPENDENCE ON THE GRADIENT
Electronic Journal of Differential Equations, Vol. 2012 (2012), No. 139, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SEMILINEAR ELLIPTIC
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 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 informationChapter 5 Fast Multipole Methods
Computational Electromagnetics; Chapter 1 1 Chapter 5 Fast Multipole Methods 5.1 Near-field and far-field expansions Like the panel clustering, the Fast Multipole Method (FMM) is a technique for the fast
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 informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
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 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 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 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 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 informationMechanical quadrature methods and extrapolation for solving nonlinear boundary Helmholtz integral equations
Appl. Math. Mech. -Engl. Ed., 32(12), 1505 1514 (2011) DOI 10.1007/s10483-011-1519-7 c Shanghai University and Springer-Verlag Berlin Heidelberg 2011 Applied Mathematics and Mechanics (English Edition)
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationSIMULTANEOUS AND NON-SIMULTANEOUS BLOW-UP AND UNIFORM BLOW-UP PROFILES FOR REACTION-DIFFUSION SYSTEM
Electronic Journal of Differential Euations, Vol. 22 (22), No. 26, pp. 9. ISSN: 72-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SIMULTANEOUS AND NON-SIMULTANEOUS
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 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 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 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 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 informationANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS INVOLVING A NEW FRACTIONAL DERIVATIVE WITHOUT SINGULAR KERNEL
Electronic Journal of Differential Equations, Vol. 217 (217), No. 289, pp. 1 6. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS
More informationUNIFORM BOUNDS FOR BESSEL FUNCTIONS
Journal of Applied Analysis Vol. 1, No. 1 (006), pp. 83 91 UNIFORM BOUNDS FOR BESSEL FUNCTIONS I. KRASIKOV Received October 8, 001 and, in revised form, July 6, 004 Abstract. For ν > 1/ and x real we shall
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 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 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 informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationMemoirs on Differential Equations and Mathematical Physics
Memoirs on Differential Equations and Mathematical Physics Volume 51, 010, 93 108 Said Kouachi and Belgacem Rebiai INVARIANT REGIONS AND THE GLOBAL EXISTENCE FOR REACTION-DIFFUSION SYSTEMS WITH A TRIDIAGONAL
More informationComplete Moment Convergence for Sung s Type Weighted Sums of ρ -Mixing Random Variables
Filomat 32:4 (208), 447 453 https://doi.org/0.2298/fil804447l Published by Faculty of Sciences and Mathematics, Uversity of Niš, Serbia Available at: http://www.pmf..ac.rs/filomat Complete Moment Convergence
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 informationXiyou Cheng Zhitao Zhang. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 34, 2009, 267 277 EXISTENCE OF POSITIVE SOLUTIONS TO SYSTEMS OF NONLINEAR INTEGRAL OR DIFFERENTIAL EQUATIONS Xiyou
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 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 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 informationError estimates for moving least square approximations
Applied Numerical Mathematics 37 (2001) 397 416 Error estimates for moving least square approximations María G. Armentano, Ricardo G. Durán 1 Departamento de Matemática, Facultad de Ciencias Exactas y
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
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 informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationNew Results for Second Order Discrete Hamiltonian Systems. Huiwen Chen*, Zhimin He, Jianli Li and Zigen Ouyang
TAIWANESE JOURNAL OF MATHEMATICS Vol. xx, No. x, pp. 1 26, xx 20xx DOI: 10.11650/tjm/7762 This paper is available online at http://journal.tms.org.tw New Results for Second Order Discrete Hamiltonian Systems
More informationSPECTRAL APPROXIMATION TO A TRANSMISSION EIGENVALUE PROBLEM AND ITS APPLICATIONS TO AN INVERSE PROBLEM
SPECTRAL APPROXIMATIO TO A TRASMISSIO EIGEVALUE PROBLEM AD ITS APPLICATIOS TO A IVERSE PROBLEM JIG A JIE SHE 2,3 Abstract. We first develop an efficient spectral-galerkin method and an rigorous error analysis
More informationColloq. Math. 145(2016), no. 1, ON SOME UNIVERSAL SUMS OF GENERALIZED POLYGONAL NUMBERS. 1. Introduction. x(x 1) (1.1) p m (x) = (m 2) + x.
Colloq. Math. 145(016), no. 1, 149-155. ON SOME UNIVERSAL SUMS OF GENERALIZED POLYGONAL NUMBERS FAN GE AND ZHI-WEI SUN Abstract. For m = 3, 4,... those p m (x) = (m )x(x 1)/ + x with x Z are called generalized
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 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 informationSpectral Collocation Methods for Differential-Algebraic Equations with Arbitrary Index
J Sci Comput (014) 58:499 516 DOI 101007/s10915-013-9755-3 Spectral Collocation Methods for Differential-Algebraic Equations with Arbitrary Index Can Huang Zhimin Zhang Received: 3 December 01 / Revised:
More informationNote on the Chen-Lin Result with the Li-Zhang Method
J. Math. Sci. Univ. Tokyo 18 (2011), 429 439. Note on the Chen-Lin Result with the Li-Zhang Method By Samy Skander Bahoura Abstract. We give a new proof of the Chen-Lin result with the method of moving
More informationNumerische Mathematik
Numer. Math. (997) 76: 479 488 Numerische Mathematik c Springer-Verlag 997 Electronic Edition Exponential decay of C cubic splines vanishing at two symmetric points in each knot interval Sang Dong Kim,,
More 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 informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
More 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 informationApplied Numerical Mathematics
Applied Numerical Mathematics 131 (2018) 1 15 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum The spectral-galerkin approximation of nonlinear eigenvalue
More informationNon-degeneracy of perturbed solutions of semilinear partial differential equations
Non-degeneracy of perturbed solutions of semilinear partial differential equations Robert Magnus, Olivier Moschetta Abstract The equation u + F(V (εx, u = 0 is considered in R n. For small ε > 0 it is
More informationON THE GLOBAL EXISTENCE OF A CROSS-DIFFUSION SYSTEM. Yuan Lou. Wei-Ming Ni. Yaping Wu
DISCRETE AND CONTINUOUS DYNAMICAL SYSTEMS Volume 4, Number 2, April 998 pp. 93 203 ON THE GLOBAL EXISTENCE OF A CROSS-DIFFUSION SYSTEM Yuan Lou Department of Mathematics, University of Chicago Chicago,
More informationFast Numerical Methods for Stochastic Computations
Fast AreviewbyDongbinXiu May 16 th,2013 Outline Motivation 1 Motivation 2 3 4 5 Example: Burgers Equation Let us consider the Burger s equation: u t + uu x = νu xx, x [ 1, 1] u( 1) =1 u(1) = 1 Example:
More informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationarxiv: v1 [math.ap] 18 Jan 2019
manuscripta mathematica manuscript No. (will be inserted by the editor) Yongpan Huang Dongsheng Li Kai Zhang Pointwise Boundary Differentiability of Solutions of Elliptic Equations Received: date / Revised
More informationUNIFORM DECAY OF SOLUTIONS FOR COUPLED VISCOELASTIC WAVE EQUATIONS
Electronic Journal of Differential Equations, Vol. 16 16, No. 7, pp. 1 11. ISSN: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu UNIFORM DECAY OF SOLUTIONS
More informationSeries Solution of Weakly-Singular Kernel Volterra Integro-Differential Equations by the Combined Laplace-Adomian Method
Series Solution of Weakly-Singular Kernel Volterra Integro-Differential Equations by the Combined Laplace-Adomian Method By: Mohsen Soori University: Amirkabir University of Technology (Tehran Polytechnic),
More informationEXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM. Saeid Shokooh and Ghasem A. Afrouzi. 1. Introduction
MATEMATIČKI VESNIK MATEMATIQKI VESNIK 69 4 (217 271 28 December 217 research paper originalni nauqni rad EXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM Saeid Shokooh and Ghasem A.
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
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 informationOn the solution of integral equations of the first kind with singular kernels of Cauchy-type
International Journal of Mathematics and Computer Science, 7(202), no. 2, 29 40 M CS On the solution of integral equations of the first kind with singular kernels of Cauchy-type G. E. Okecha, C. E. Onwukwe
More informationAn Efficient Numerical Method for Solving. the Fractional Diffusion Equation
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 1-12 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 An Efficient Numerical Method for Solving the Fractional
More informationResearch Statement. James Bremer Department of Mathematics, University of California, Davis
Research Statement James Bremer Department of Mathematics, University of California, Davis Email: bremer@math.ucdavis.edu Webpage: https.math.ucdavis.edu/ bremer I work in the field of numerical analysis,
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 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 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 information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
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 informationLOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S PARABOLIC-HYPERBOLIC COMPRESSIBLE NON-ISOTHERMAL MODEL FOR LIQUID CRYSTALS
Electronic Journal of Differential Equations, Vol. 017 (017), No. 3, pp. 1 8. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu LOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S
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 informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationMinimization problems on the Hardy-Sobolev inequality
manuscript No. (will be inserted by the editor) Minimization problems on the Hardy-Sobolev inequality Masato Hashizume Received: date / Accepted: date Abstract We study minimization problems on Hardy-Sobolev
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationAnalysis of Two-Grid Methods for Nonlinear Parabolic Equations by Expanded Mixed Finite Element Methods
Advances in Applied athematics and echanics Adv. Appl. ath. ech., Vol. 1, No. 6, pp. 830-844 DOI: 10.408/aamm.09-m09S09 December 009 Analysis of Two-Grid ethods for Nonlinear Parabolic Equations by Expanded
More informationSome lecture notes for Math 6050E: PDEs, Fall 2016
Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.
More informationExistence and Multiplicity of Solutions for a Class of Semilinear Elliptic Equations 1
Journal of Mathematical Analysis and Applications 257, 321 331 (2001) doi:10.1006/jmaa.2000.7347, available online at http://www.idealibrary.com on Existence and Multiplicity of Solutions for a Class of
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 informationEXISTENCE OF SOLUTIONS FOR CROSS CRITICAL EXPONENTIAL N-LAPLACIAN SYSTEMS
Electronic Journal of Differential Equations, Vol. 2014 (2014), o. 28, pp. 1 10. ISS: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu EXISTECE OF SOLUTIOS
More informationDirect method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions
Journal of Computational and Applied Mathematics 22 (28) 51 57 wwwelseviercom/locate/cam Direct method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions
More informationNon-radial solutions to a bi-harmonic equation with negative exponent
Non-radial solutions to a bi-harmonic equation with negative exponent Ali Hyder Department of Mathematics, University of British Columbia, Vancouver BC V6TZ2, Canada ali.hyder@math.ubc.ca Juncheng Wei
More informationSpectral Representation of Random Processes
Spectral Representation of Random Processes Example: Represent u(t,x,q) by! u K (t, x, Q) = u k (t, x) k(q) where k(q) are orthogonal polynomials. Single Random Variable:! Let k (Q) be orthogonal with
More informationAn Improved F-Expansion Method and Its Application to Coupled Drinfel d Sokolov Wilson Equation
Commun. Theor. Phys. (Beijing, China) 50 (008) pp. 309 314 c Chinese Physical Society Vol. 50, No., August 15, 008 An Improved F-Expansion Method and Its Application to Coupled Drinfel d Sokolov Wilson
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 informationA collocation method for solving the fractional calculus of variation problems
Bol. Soc. Paran. Mat. (3s.) v. 35 1 (2017): 163 172. c SPM ISSN-2175-1188 on line ISSN-00378712 in press SPM: www.spm.uem.br/bspm doi:10.5269/bspm.v35i1.26333 A collocation method for solving the fractional
More informationCLASSROOM NOTES PART II: SPECIAL TOPICS. APM526, Spring 2018 Last update: Apr 11
CLASSROOM NOTES PART II: SPECIAL TOPICS APM526, Spring 2018 Last update: Apr 11 1 Function Space Methods General Setting: Projection into finite dimensional subspaces t u = F (u), u(t = 0) = u I, F : B
More informationA new method for the solution of scattering problems
A new method for the solution of scattering problems Thorsten Hohage, Frank Schmidt and Lin Zschiedrich Konrad-Zuse-Zentrum Berlin, hohage@zibde * after February 22: University of Göttingen Abstract We
More informationYuqing Chen, Yeol Je Cho, and Li Yang
Bull. Korean Math. Soc. 39 (2002), No. 4, pp. 535 541 NOTE ON THE RESULTS WITH LOWER SEMI-CONTINUITY Yuqing Chen, Yeol Je Cho, and Li Yang Abstract. In this paper, we introduce the concept of lower semicontinuous
More information