A SPECTRAL COLLOCATION METHOD FOR EIGENVALUE PROBLEMS OF COMPACT INTEGRAL OPERATORS
|
|
- Maximilian Warren
- 6 years ago
- Views:
Transcription
1 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 Communicated by Da Xu ABSTRACT. We propose and analyze a new spectral collocation method to solve eigenvalue problems of compact integral operators, particularly, piecewise smooth operator kernels and weakly singular operator kernels of the form 1/ t s µ, 0 < μ < 1. We prove that the convergence rate of eigenvalue approximation depends upon the smoothness of the corresponding eigenfunctions for piecewise smooth kernels. On the other hand, we can numerically obtain a higher rate of convergence for the above weakly singular kernel for some μ s even if the eigenfunction is not smooth. Numerical experiments confirm our theoretical results. 1. Introduction. We consider numerical approximation of the eigenvalue problem for a compact integral operator T on a Banach space. Recent years have witnessed a revitalization of this field, and various methods are applied to solve the problem. The Galerkin, Petrov-Galerkin, collocation, Nyström and degenerate kernel methods have been extensively studied for the approximation of eigenvalues and eigenvectors of integral operators. The results are well documented in the literature. Here, we mention a few related to our current work. As early as 1967, Atkinson proved a general theorem showing the convergence of numerical eigenvalues and eigenvectors to those of compact integral operators []. In 1975, he further obtained a convergence rate for the approximation [3], based upon which Osborn established a general spectral approximation theory for compact operators, when a sequence of {T n } approximates T in a collectively compact manner. The analysis of [3, 17] covers many methods and provides a basis for the convergence analysis of our method. In [13], Dellwo and Friedman proposed 010 AMS Mathematics subject classification. Primary 47A10, 47A58, 65J99, 65MR0. Keywords and phrases. Eigenvalue problem, spectral collocation method, weakly singular kernel, integral operator, super-geometric convergence. Received by the editors on November 11, 011. DOI:10.116/JIE Copyright c 013 Rocky Mountain Mathematics Consortium 79
2 80 CAN HUANG, HAILONG GUO AND ZHIMIN ZHANG a new approach by solving a polynomial eigenvalue problem of a higher degree, based upon which Alam et al. [1] obtained an accelerated spectral approximation for eigenelements. Kulkarni [16] introduced another method by involving a new approximation operator T n and obtained a high-order convergence rate. In addition, a multiscale method was discussed in [11]. Comprehensive studies for eigenvalue problem can be found in [5, 9, 1]. For the numerical solution of integral equations or integro-differential equations, interested readers are referred to [4, 7]. In this article, we approximate eigenfunctions by some appropriate orthogonal polynomial expansions. In a different manner from previous methods in the literature we find the exact integration when calculating the convolution of the singular kernel with the orthogonal polynomials. The key ingredients here are some special identities. By doing so, we: 1) avoid large numerical quadrature errors accumulated with the singular kernels and thereby obtain higher accuracy for eigenvalue approximations, and ) avoid product integration methods and therefore reduce the computational cost. Furthermore, if the kernel is positive definite and piecewisely smooth, a refined result can be obtained. To fix the idea, we consider problems of the form (1.1) 0 k(t, s)u(s) ds = λu(t), t [0, 1], where k(t, s) = t s μ for 0 <μ<1, k(t, s) is piecewisely smooth or smooth. We will develop algorithms for all three types of problems separately. This paper is organized as follows. In Section, some preliminary knowledge is given. In Section 3, algorithms for all types of equations are listed. Section 4 is devoted to convergence analysis of algorithms. Finally, we illustrate our theories with numerical examples in Section 5. Throughout the paper, C stands for a generic constant that is independent of collocation points p but which may depend upon the index μ and the number of pieces a piecewise kernel has.. Preliminaries. Let T : X X be a compact linear operator on a Banach space X and σ(t )andρ(t ) the spectrum and resolvent of T, respectively. Let λ be a nonzero eigenvalue of T with multiplicity m, and let Γ be a circle centered at λ which lies in ρ(t ) and which encloses
3 COMPACT INTEGRAL OPERATOR EIGENVALUE PROBLEMS 81 no other points in σ(t ). Then the spectral projection associated with T and λ is defined by E = 1 (T zi) 1 dz πi Γ and max z Γ (T zi) 1 C. Let {T n } be a sequence of operators in B(X) thatconvergestot in a collective way, i.e., the set {T n x : x 1, n=1,,...} is sequentially compact. For n large enough, Γ ρ(t n ) and the associated projection, E n = 1 (T n zi) 1 dz πi Γ exists and max z Γ (T n zi) 1 C. Clearly, dim (E) =dim(e n )= m and T n E n = E n T n. Furthermore, the spectrum of T n inside Γ contains m approximations of λ, i.e., λ n,1,λ n,,...,λ n,m, counted according to their algebraic multiplicities [9, 17]. Let λ n = λ n,1 + λ n, + + λ n,m. m Then we have the following theorem. Theorem.1 [17]. For all n sufficiently large, λ λ n C (T T n ) R(E), where R(E) is the range of the projection E. This is a rather general result. We may refine the result if the kernel is positive definite. Let (.1) a(u, v) = 0 0 k(t, s)u(s)v(t) ds dt, b(u, v) = 0 u(t)v(t) dt, where v is a test function in the L space V. If the bilinear operator a(u, v) is coercive, then we can list eigenvalues of T by λ 1 λ λ 3 0,
4 8 CAN HUANG, HAILONG GUO AND ZHIMIN ZHANG with zero the only possible cluster point. Let us consider a numerical approximation of the first eigenpair (λ, u). Let (λ p,v p ) be their Galerkin approximation, and let u p be the Legendre expansion of u. Wehave a(u, u) (.) λ = b(u, u) =sup a(v, v) v V b(v, v), λ p = a(v p,v p ) b(v p,v p ) =max a(v, v) v P p b(v, v). Here P p is the polynomial space with degree no more than p. Denote λ p = a(u p,u p )/b(u p,u p ); then we have the following lemma. Lemma.. Let λ, λ p and λ p be defined as above and a(u, v) coercive. Then (.3) 0 λ λ p λ λ p = λ u u p b u b u u p a u. b Proof. From[5, page 701, Lemma 9.1], we have (.4) 0 ν p ν ν p ν u u p b u a ν u u p a u, a where ν =1/λ, ν p =1/λ p and ν p =1/ λ p. Hence, (.5) 0 λ λ p λ Using the fact that λ p u u p b u a λ p λ u u p a u. a we derive (.3) from (.5). a(u p,u p )=λ p b(u p,u p ), Next, we introduce some identities, which will be essential in this paper. Towards this end, we define the class of Jacobi polynomials P n (α,β) (x). Under the normalization P (α,β) k (1) = ( ) k+α k, one has the expression, namely, (.6) P (α,β) k (x) = 1 k k l=0 ( k + α k l )( k + β l ) (x 1) l (x +1) k l.
5 COMPACT INTEGRAL OPERATOR EIGENVALUE PROBLEMS 83 Jacobi polynomials satisfy the three-term recursive relations: (.7) P (α,β) 0 (x) =1, P (α,β) 1 (x) = 1 [(α β)+(α + β +)x], a 1,k P (α,β) k+1 (x) =a,kp (α,β) k (x) a 3,k P (α,β) k 1 (x), where (.8) a 1,k =(k +1)(k + α + β + 1)(k + α + β), a,k =(k + α + β +1)(α β )+xγ(k + α + β +3)/Γ(k + α + β), a 3,k =(k + α)(k + β)(k + α + β +). Especially if α =0andβ = 0, Jacobi polynomials become Legendre polynomials. Lemma.3 [19]. Let a, b be positive constants and L n (x) the Legendre polynomials with degree n on [ 1, 1]. Then (.9) (.10) b a b ) ds = n! (b a) α P n (α, α) (α) n+1 ( s (s a) α 1 L n b b <a<b, α>0, ) ds = n! (b + a) β P n ( β,β) (β) n+1 s s) a(b β 1 L n( a a <b<a, β>0, ( a b ( b a ), ), where (k) n+1 = k(k +1) (k + n). Specifically, if we choose a =1,b = x, β =1 μ in (.10), then we obtain (.11) x 1 L n (t) (x t) μ dt = n! (1 + x) 1 μ P n (μ 1,1 μ) (x), (1 μ) n+1 and a = x, b =1,α =1 μ in (.9), we arrive at (.1) x L n (t) (t x) μ dt = n! (1 x) 1 μ P n (1 μ,μ 1) (x). (1 μ) n+1
6 84 CAN HUANG, HAILONG GUO AND ZHIMIN ZHANG Remark 1. We use identities (.11) and (.1) in our algorithm for weakly singular kernels after we expand eigenvectors by Legendre polynomials. Lemma.4 [18]. Let α> 1, β> 1 and 0 <ν<1. Then, for 1 <x<1, (.13) where (.14) (.15) 1 (1 t) α (1 + t) β P m (α,β) (t) dt x t ν = cos(πν/)φ 1(x)+cosπ((ν/) β)φ (x), m =0, 1,,..., Γ(ν)cos(πν/) Φ 1 (x) = Φ (x) = Γ(m + α +1)Γ(m + ν)γ(β ν +1)( 1)m α β+ν 1 Γ(m + α + β ν +)m! F 1 ( m + ν, ν m α β 1; β + ν; 1+x Γ(m + β +1)Γ(ν β 1)( 1)m+1 α (1 + x) ν β 1 m! F 1 ( m + β +1, m α; β ν +; 1+x Here, F 1 (a, b; c; z) is known as Gauss s hypergeometric functions. ). ), For the sake of convergence analysis, we need to introduce the error estimate of Gauss quadrature. Lemma.5 [1]. Let f C n, x i and w i be the Gauss points and their corresponding quadrature weights on the interval [a, b]. Then (.16) b n f(x) dx w i f(x i )= (b a)n+1 (n!) 4 (n + 1)[(n)!] 3 f (n) (ξ), ξ (a, b). a i=0
7 COMPACT INTEGRAL OPERATOR EIGENVALUE PROBLEMS Algorithms. In this section, we develop algorithms for eigenproblem with all three kinds of kernels mentioned before. Models that we consider in this article are: (1) Weakly singular kernels (3.1) λy(t) = () Piecewise smooth kernels 0 y(s) ds, 0 <μ<1, t [0, 1]; t s μ (3.) λy(t) = k(t, s)y(s) ds, t [0, 1], 0 { t s/ if 0 t s 1, where k(t, s) = s/ if 0 s<t 1; (3) Smooth kernels (3.3) λy(t) = 0 e st y(s) ds, t [0, 1] The first algorithm for (3.1). It is clear that (3.1) is equivalent to t y(s) 1 (3.4) λy(t) = 0 (t s) μ ds + y(s) t (s t) μ ds. We make a change of variable t =(1+x)/ andobtain (3.5) (1+x)/ ( μ +x 1 ( s) y(s) ds+ s 1+x ) μ y(s) ds=λu(x), 0 (1+x)/ where x [ 1, 1] and u(x) =y((1 + x)/). Next, we make another change of variable, s =(1+τ)/ andreach (3.6) ( ) 1 μ x (x τ) μ u(τ) dτ + 1 ( 1 ) 1 μ x (τ x) μ u(τ) dτ =λu(x), x [ 1, 1].
8 86 CAN HUANG, HAILONG GUO AND ZHIMIN ZHANG Let u p (x) = p c jl j (x) be the approximation of u(x). Obviously, the c j s satisfy the equation (3.7) ( ) 1 μ xi c j 1 L j (τ) (x i τ) μ dτ + ( ) 1 μ 1 c j x i Substituting (.11) and (.1) into (3.7), we obtain L j (τ) (τ x i ) μ dτ = λ p c j L j (x i ). (3.8) [( ) 1 μ 1 j! c j (1 + x i ) 1 μ P (μ 1,1 μ) j (x i ) (1 μ) j+1 ( ) 1 μ ] 1 j! + (1 x i ) 1 μ P (1 μ,μ 1) j (x i ) (1 μ) j+1 = λ p c j L j (x i ), i =0,...,p. If we write ( ) 1 μ 1 j! a ij = (1 + x i ) 1 μ P (μ 1,1 μ) j (x i ) (1 μ) j+1 ( ) 1 μ 1 j! + (1 x i ) 1 μ P (1 μ,μ 1) j (x i ) (1 μ) j+1 b ij = L j (x i ), then we have AC p = λ p BC p, where A = (a ij ), B = b ij,c p = (c 0,c 1,...,c p ) T. 3.. The second algorithm for (3.1). From [0, Theorem 1], we assume that the first true eigenvector is of the form (3.9) y(t) = d 1 t 1 μ + d (1 t) 1 μ + a smoother function φ(t). Hence, we approximate the eigenvector by u p (t) =d 1 t 1 μ +d (1 t) 1 μ + p c jp j (t), where P j (t) is the shifted Legendre polynomial on [0, 1],
9 COMPACT INTEGRAL OPERATOR EIGENVALUE PROBLEMS 87 j =0, 1,...,p. Substituting it into (3.1) and taking the same change of variables as the previous algorithm, we obtain ( ) μ 1 (1 + τ) (3.10) (d 1 μ 1 1 x τ μ dτ + d (1 τ) 1 μ ) 1 x τ μ dτ ( ) 1 μ ( x ) + c j (x τ) μ L j (τ) dτ + (τ x) μ L j (τ) dτ 1 x ( ) 1 μ ( ) 1 μ 1+x 1 x = d 1 + d + c j L j (x), x [ 1, 1]. From Lemmas.3 and.4 and (3.10), we obtain (3.11) ( 1 ) μ ( Γ( μ) μ Γ(3 μ) F 1 μ, μ ; μ 1; 1+x ) i d 1 ( ) μ ( ( 1 1 Γ( μ)γ(μ)γ(1 μ) + Γ(μ) μ F 1 μ, μ ; μ; 1+x i Γ(3 μ) ( Γ(μ 1) μ 1 (1 + x i ) μ 1 F 1 1,μ 1; μ; 1+x i [( ) 1 μ 1 j! + c j (1 + x i ) 1 μ P (μ 1,1 μ) j (x i ) (1 μ) j+1 = λ p ( + ( ) 1 μ 1 j! c j L j (x i )+d 1 ( 1+xi (1 x i ) 1 μ P (1 μ,μ 1) j (x i ) (1 μ) j+1 ) 1 μ ( 1 xi + d i =0,...,p+. ) 1 μ ), ] ) )) d Note that the first hypergeometric function is not well defined when μ =1/. However, the integration of the two singular terms with the
10 88 CAN HUANG, HAILONG GUO AND ZHIMIN ZHANG kernel are simpler, in which case, the linear system is (3.1) ( π(1 + xi ) 1 xi x i 4 ( ) 1 xi log 1+ 1+x ( )) i 1+xi log d 1 4 ( 1+xi + x ( ) i 1 tanh 1 1+xi π(x ) i 1) d 4 [( ) 1 μ 1 j! + c j (1 + x i ) 1 μ P (μ 1,1 μ) j (x i ) (1 μ) j+1 + ( ) 1 μ 1 j! ( 1+xi 1 xi = λ p d 1 + d + (1 x i ) 1 μ P (1 μ,μ 1) j (x i ) (1 μ) j+1 ) c j L j (x i ), i =0,...,p+. ] 3.3. Algorithm for (3.). We make the change of variable as before and let u(x) =y((1 + x)/). We obtain (3.13) λu(x) = 1 1+τ u(τ) dτ x (x τ)u(τ) dτ. Let u p (x) = p c jl j (x) be the approximation of u(x). Then the c j s satisfy (3.14) 1+τ λ p c j L j (x i )= c j L j (τ) dτ + c j (x i τ)l j (τ) dτ, 1 8 x i i =0,...,p,
11 COMPACT INTEGRAL OPERATOR EIGENVALUE PROBLEMS 89 i.e., (3.15) λ p c j L j (x i )= + 1 x i 8 ( c0 4 + c ) 1 1 c j k=0 ( 1+xi w k (x i x k )L j i =0,...,p. + 1 x i x k ), Here, the numerical integration is exact. The scheme is of the form where AC p = λ p BC p, b ij = L j (x i ) (1 x i )/8 p k=0 w k(x i x k )L j ((1+x i )/+(1 x i )/x k ) if j 0, 1, (1 x i )/8 p k=0 a ij = w k(x i x k )L j ((1+x i )/+(1 x i )/x k )+1/4 if j =0; (1 x i )/8 p k=0 w k(x i x k )L j ((1+x i )/+(1 x i )/x k )+1/1 if j = Algorithm for (3.3). Substituting the Legendre expansion y(t) = p y jl j (t) into (3.3) and collocating at n Gaussian points, we have (3.16) y j e sti L j (s) ds = λ y j L j (t i ), 0 i =1,,...,n. The matrix form of (3.16) is (3.17) Ky = λly,
12 90 CAN HUANG, HAILONG GUO AND ZHIMIN ZHANG where (3.18) K ij = 0 e sti L j (s) ds, L ij = L j (t i ), y =(y 1,y,...,y n ) T. K ij can be calculated by the n-point Gaussian quadrature (3.19) e sti L j (S) ds 0 l=0 e s lt i L j (s l )w l, s k = t k. 4. Convergence analysis. Let L k be the standard Legendre polynomial of degree k, and let π p f P p [ 1, 1] interpolate a smooth function f at (p +1)-Gauss points: 1 <x 0 < <x p < 1. Let T k be the first kind Chebyshev polynomial of degree k. Then the remainder of the interpolation is (4.1) f(x) π p f(x) =f[x 0,x 1,...,x p,x]ν(x), where ν(x) =(x x 0 )(x x 1 ) (x x p ). Note that (4.) d n L n (x) = 1 n n! dx n (x 1) n = 1 d n n ( ) n n n! dx n x (n j) ( 1) j j = 1 n ( ) n n (n j)(n j 1) (n j n+1)x n j ( 1) j. n! j From the term with j = 0, we get the leading coefficient ( ) 1 n (4.3) n (n)(n 1) (n n +1)( 1) 0 = (n)! n! 0 n (n!) By the Stirling formula, (4.4) (n)! n (n!) (n/e)n 4πn n [(n/e) n πn] =n.
13 COMPACT INTEGRAL OPERATOR EIGENVALUE PROBLEMS 91 Hence, (4.5) f(x) π p f(x) f[x 0,x 1,...,x p,x] p+1 L p+1 (x). If f C p+1 [ 1, 1], the divided difference (4.6) f[x 0,x 1,...,x p,x]= f (p+1) (ξ x ), ξ x ( 1, 1). (p +1)! The result can be concluded as the following theory. Theorem 4.1. Ck!R k, then (1) If y(t) satisfies condition (K): y (k) L [0,1] (4.7) y π p y L [0,1] C (4R) p+1 ; () If y(t) satisfies condition (M): y (k) L [0,1] CM k, we have (4.8) y π p y L [0,1] C ( ) p+1 em. p +1 4(p +1) Proof. We make the change of variables t = 1+x, s = 1+τ, x,τ [ 1, 1], and let u(x) =y((1 + x)/). Then the result for y under condition (K) follows directly from (4.6) and the fact that dt =(1/)dx. If y satisfies condition (M), by applying the Stirling s formula, (4.9) y π p y L [0,1] = u π p u L [ 1,1] CM p+1 (4R) p+1 (p +1)! cm p+1 π(p + 1)(4(p +1)/e) p+1 = C ( ) p+1 em. p +1 4(p +1)
14 9 CAN HUANG, HAILONG GUO AND ZHIMIN ZHANG For non-smooth functions, we need some other estimates. Theorem 4. [8]. (1) For any f H k ( 1, 1), (4.10) f π p f L ( 1,1) Cp k f H k,p ( 1,1). () For any f Hw k ( 1, 1), (4.11) f π c pf L ( 1,1) Cp k f H k,p w ( 1,1), where two seminorms are defined by ( k f H k,p ( 1,1) = s=min(k,p+1) f H k,p w ( 1,1) = ( k s=min(k,p+1) f (s) L ( 1,1)) 1/, f (s) L w ( 1,1) ) 1/, and the weight w(x) =(1 x) 1/ (1+x) 1/ and πp c is the interpolatory operator on Chebyshev points. Let R(E) andr(e p ) be the range of E and E p, respectively. Define π p : R(E) R(E p ) as an interpolatory projection by π p (t) = p ξ jl j (t) andξ j is determined by ξ j L j (t i )=x(t i ), i =0,...,p. Then our algorithms can be written as (4.1) T p u p = λ p u p, where T p = π p T. Theorem 4.3. Let y be the exact first eigenvector and T acompact operator in (3.1), (3.) or (3.3) and T p defined as above. Then
15 COMPACT INTEGRAL OPERATOR EIGENVALUE PROBLEMS 93 (1) If y H k (0, 1), (4.13) λ λ p C (p) k ; () Furthermore, if y satisfies condition (K), (4.14) λ λ p C (4R) p+1 ; (3) Furthermore, if y satisfies condition (M), (4.15) λ λ p C ( ) p+1 em. p +1 4(p +1) Proof. The result follows directly from Theorems.1, 4.1, 4. and [16, Theorem.]. To make the paper self contained, we put the proof here. Let Êp = E p R(E) : R(E) R(E p ). Then, for large p, Êp is bijective and Ê 1 p [17]. Define T = T R(E),and T p := Ê 1T pêp. Then p (4.16) λ λ p = 1 m trace ( T T p ) T T p = Ê 1 n (ÊpT ÊT p) C (T T p ) R(E) = C (I π p )T R(E). Since Tu is smoother than u, see[10, 15]; then the result follows. Theorem 4.4. Let λ and λ p be the exact eigenvalue and its numerical approximation of a positive definite operator T whose kernel is a piecewise smooth function, respectively. Then (1) if u satisfies condition (K), ( (4.17) λ λ p C 1 (4R) p+ + e p ) ; p p 3/ 6p
16 94 CAN HUANG, HAILONG GUO AND ZHIMIN ZHANG () if u satisfies condition (M), ( ( ) p+ 1 em e p ) (4.18) λ λ p C + ; p +1 4(p +1) p p 3/ 6p (3) if u H k [0, 1], ( 1 (4.19) λ λ p C (p) k + e p ), p p 3/ 6p Proof. By our algorithms, we have (4.0) 0 k(t i,s)u p (s) ds = λ p u p (t i ), where t i are (p + 1)-Gauss points on [0,1]. Multiplying both sides by L j (t i )w i and summing up from 0 to p, we obtain (4.1) k(t i,s)u p (s)l j (t i )w i ds = λ p u p (t i )L j (t i )w i. 0 Here, w i are weights of the Gauss quadrature. If we write Ã=( 0 0 k(t, s)l j(s)l i (t) dsdt) ij, B =( 0 L j(t)l i (t) dt) ij and recall that u p (x) = p i=0 ũil i (x), we obtain Ãũ = λ p Bũ, where ũ =[ũ 0, ũ 1,...,ũ p ] T. However, for most cases, we can only apply numerical quadrature to find elements of à and B. If the kernel is piecewise smooth, we apply the Gauss quadrature piece by piece. Therefore, the system that we actually solve is (4.) Au = λ p Bu. Now we are ready to analyze errors of eigenvalue approximations.
17 COMPACT INTEGRAL OPERATOR EIGENVALUE PROBLEMS 95 First, we analyze the case when the kernel is a linear piecewise polynomial. Noting the fact that (p + 1)-Gauss quadrature is exact for all polynomials of degree less than or equal to p +1,wehave (4.3) A = Ã, and B = B. Here, the integration is piecewise, so is the numerical integration. Denote the arithmetic mean of the approximation of λ by λ p again, if it is a multiple eigenvalue. We derive from Lemma. that (4.4) λ λ p = λ λ p 1/4 p+ if u satisfies condition (K); C 1/(p +1)(eM/(4(p + 1))) p+ if u satisfies condition (M); 1/(p) k if u H k [0, 1]. If the kernel is piecewise smooth or smooth, from the analysis of the previous case, we only need to estimate A Ã since B = B. Ifwewrite the remainder of Gaussian quadrature as ε, then (4.5) A Ã Cε Here, we define E<F if and only if (E) ij < (F ) ij. By the error estimate of the Gauss quadrature and Stirling s formula, we have (4.6) ( [p!] 4 ) ε C (p + 1)[(p)!] 3 ( [ πp(p/e) p ] 4 ) C (p +1)[ π(p)(p/e) p ] 3 ( e p ) C. p p+1/ 6p Hence, (4.7) ( A Ã pe p ) n C, p p+1/ 6p n =1,.
18 96 CAN HUANG, HAILONG GUO AND ZHIMIN ZHANG Clearly, (4.) is equivalent to (4.8) B 1 Au = λ p u. Thus, (4.9) B 1 à B 1 A n B 1 n à A n ( p e p ) C, n =1,, p p+1/ 6p by noting that B = B = diag (1, 1/3,...,1/(p + 1)). Therefore, by a perturbation theory, see [9, page 30], we have ( (4.30) λ p λ p C e p p p 3/ 6p Denote the arithmetic mean of the approximation of λ by λ p again, if it is a multiple eigenvalue. We derive from Lemma. and (4.30) that λ λ p λ λ p + λ p λ p ((1/4 p+ )+(e p /(p p 1/ 6p ))) if u satisfies condition (K); (1/(p +1)(eM/(4(p + 1))) p+ +(e p /(p p 1/ 6p ))) C if u satisfies condition (M); ((1/(p) k )+(e p /(p p 1/ 6p ))) if u H k [0, 1]. ). Remark. Theorem 4.4 shows that, although numerical integration contributes to the error of eigenvalue approximation, it is trivial compared with truncation error for our method. Hence, in our numerical experiments, we ignore it for reference curves. 5. Numerical examples. In this section, we will find numerical approximations to solutions of some examples to demonstrate our theory.
19 COMPACT INTEGRAL OPERATOR EIGENVALUE PROBLEMS Our Method Refence Curve 10 0 Our Method Refence Curve λ λ p λ λ p p p FIGURE 1. (Left) Linear piecewise kernel. (Right) Kernel: e st. TABLE 1. Example 5.1: λ λ p. p error 5.98e e e e-11 p 7 8 error e e-15 TABLE. Example 5.: λ p (The first algorithm). p p Example 5.1 [1]. We consider a problem with form (3.). Then each λ j =1/((j 1) π ), j =1,,..., is an eigenvalue of T of algebraic multiplicity m =. Let λ denote the arithmetic mean of the two eigenvalues of T p to the largest two eigenvalues λ =1/π. Numerical errors are presented in Table 1 and the left part of Figure 1, from which we see that the error decays super-geometrically. Here Reference Curve is the graph of ( ) p+ 1 eπ f(p) =. 100(p +1) 4(p +1)
20 98 CAN HUANG, HAILONG GUO AND ZHIMIN ZHANG 10 3 μ=1/ μ=1/3 Our Method 10 1 Discontinuous collocation method Reference curve Slope is 14/3 Slope is 7/6 λ p λ 10 5 Slope is 10/3 λ p λ 10 8 λ p λ p p Intervals FIGURE. Kernel: t s 1/3. The first algorithm (Left), the second algorithm (Middle) and the discontinuous linear elements method (Right). TABLE 3. Example 5.: λ p (The second algorithm). p p Example 5.. Now let us consider an eigenproblem of form (3.1) with μ = 1/3. From [0], eigenfunctions belong to H (7/6) ε (0, 1), where ε is a sufficiently small positive number and we expect to obtain a convergence rate of O(p 7/3 ) based on Theorem 4.3 for the first algorithm. Here, we apply both our spectral collocation methods and the three-point Gaussian collocation on equally spaced intervals with discontinuous piecewise quadratic elements method. Unfortunately, we do not know the exact eigenvalues for such types of kernels. However, we list some of our numerical approximations in Table, and we use the numerical approximation of the second algorithm for p = 70asour exact value to obtain Figure. It is easy to see that we can only obtain a 7 digit accuracy for the first algorithm; we obtain an 11 digit of accuracy and a convergence rate of O(p 14/3 ) for the second algorithm, see Table 3. However, the convergence rate for the three-point Gaussian collocation with discontinuous piecewise quadratic element is only O(h 7/6 ). This fact also confirms results in [6], which says that the convergence rate for p-version methods doubles the convergence rate for the h-version method if the true solution is singular. Example 5.3. We consider an eigenproblem of form (3.1) with μ =1/. In this case, eigenfunctions belong to H 1 (0, 1). Again, we use both algorithms to solve it and consider the numerical approximation
21 COMPACT INTEGRAL OPERATOR EIGENVALUE PROBLEMS 99 of the second algorithm for p = 70 as the exact first eigenvalue. Numerical results are shown in Table 4, Table 5 and Figure 3. Example 5.4. Consider the eigenvalue problem of the form (3.3). We apply the algorithm in Section 3. Since the kernel is smooth, the first eigenvalue converges very fast, see Table 6 and the right part of Figure 1. In this case, Reference Curve is the graph of f(p) =1/(10(p + 1))(e/((p + 1))) p μ=1/ 10 5 μ=1/ Our Method λ p λ Slope is 3 λ p λ 10 7 Slope is p FIGURE 3. (Left) Kernel: t s 1/ (The first algorithm). (Right) Kernel: t s 1/ (The second algorithm). p TABLE 4. Example 5.3: λ p (The first algorithm). p p TABLE 5. Example 5.3: λ p (The second algorithm). p p TABLE 6. Example 5.4: λ λ p. p error e e e e-11 p 5 6 error e e-16
22 100 CAN HUANG, HAILONG GUO AND ZHIMIN ZHANG Acknowledgments. The work of the third author is supported in part by the U.S. National Science Foundation through grant DMS REFERENCES 1. R.Alam,R.P.KulkarniandB.V.Limaye,Accelerated spectral approximation, Math. Comp. 67 (1998), K.E. Atkinson, The numerical solution of the eigenvalue problem for compact integral operators, Trans. Amer. Math. Soc. 19 (1967), , Convergence rates for approximate eigenvalues of compact integral operators, SIAMJ. Numer. Anal. 1 (1975), , The numerical solution of integral equation of the second kind, Combridge University Press, Cambridge, I. Babu ska and J. Osborn, Eigenvalue problems, in Handbook of numerical analysis, Vol.:Finite element methods, P.G. Ciarlet and J.L. Lions, eds., Elsevier Science Publishers, North-Holland, I. Babu ska and M. Suri, The h-p version of the finite element method with quasiuniform meshes, Rairo Model. Math. Anal. Num. 1 (1987), H. Brunner, Collocation methods for Volterra integral and related functional equations, Cambridge University Press, Cambridge, C. Canuto, M.Y. Hussaini, A. Quarteroni and T.A. Zang, Spectral methodsfundamentals in single domains, Springer, Berlin, F. Chatelin, Spectral approximation of linear operators, Academic Press, New York, Y. Chen and T. Tang, Spectral methods for weakly singular Volterra integral equations with smooth solutions, J. Comp. Appl. Math. 33 (009), PLEASE double-check page numbers.. this does not appear correct Z. Chen, G. Nelakanti, Y. Xu and Y. Zhang, A fast collocation method for eigen-problems of weakly singular integral operators, J.Sci.Comp.41 (009), P.J. Davis and P. Rabinowitz, Methods of numerical integration, nd edition, Comp. Sci. Appl. Math., Academic Press, Orlando, D. Dellwo and M.B. Friedman, Accelerated spectral analysis of compact operators, SIAMJ. Numer. Anal. 1 (1984), R.A. Devore and L.R. Scott, Error bounds for Gaussian quadrature and weighted-l 1 polynomial approximation, SIAMJ.Numer.Anal.1 (1984), C. Huang, T. Tang and Z. Zhang, Geometric convegence of spectral collocation method for Volterra or Fredholm integral equations with weakly singular kernels, J. Comp. Math. 9 (011), R.P. Kulkarni, A new superconvergent collocation method for eigenvalue problems, Math. Comp. 75 (006), J.E. Osborn, Spectral approximation for compact operators, Math. Comp. 9 (1975),
23 COMPACT INTEGRAL OPERATOR EIGENVALUE PROBLEMS I. Podlubny, Riesz potential and Riemann-Liouville fractional integral and derivatives of Jacobi polynomials, Appl. Math. Lett. 10 (1997), A.P. Prudnikov, Yu.A. Brychkov and O.I Marichev, Integrals and series, Volume, Gordon and Breach Science, New York, G. Vainikko and A. Pedas, The properties of solutions of weakly singular integral equations, J. Austral. Math. Soc. (1981), Y.D. Yang, Finite element analysis for eigenvalue problems, People s Publication Press of Guizhou, Guiyang, China, 004 (in Chinese). Department of Mathematics, Wayne State University, Detroit, MI 480 address: huangcan@msu.edu Department of Mathematics, Wayne State University, Detroit, MI 480 address: hlguo@wayne.edu Department of Mathematics, Wayne State University, Detroit, MI 480 and Beijing Computational Science Research Center, Beijing , China address: zzhang@math.wayne.edu
Spectral 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 information1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel
1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel Xiaoyong Zhang 1, Junlin Li 2 1 Shanghai Maritime University, Shanghai,
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More 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 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 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 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 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 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 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 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 informationFinite Element Superconvergence Approximation for One-Dimensional Singularly Perturbed Problems
Finite Element Superconvergence Approximation for One-Dimensional Singularly Perturbed Problems Zhimin Zhang Department of Mathematics Wayne State University Detroit, MI 48202 Received 19 July 2000; accepted
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 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 informationFractional Calculus for Solving Abel s Integral Equations Using Chebyshev Polynomials
Applied Mathematical Sciences, Vol. 5, 211, no. 45, 227-2216 Fractional Calculus for Solving Abel s Integral Equations Using Chebyshev Polynomials Z. Avazzadeh, B. Shafiee and G. B. Loghmani Department
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 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 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 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 informationon! 0, 1 and 2 In the Zienkiewicz-Zhu SPR p 1 and p 2 are obtained by solving the locally discrete least-squares p
Analysis of a Class of Superconvergence Patch Recovery Techniques for Linear and Bilinear Finite Elements Bo Li Zhimin Zhang y Abstract Mathematical proofs are presented for the derivative superconvergence
More informationMulti-Term Linear Fractional Nabla Difference Equations with Constant Coefficients
International Journal of Difference Equations ISSN 0973-6069, Volume 0, Number, pp. 9 06 205 http://campus.mst.edu/ijde Multi-Term Linear Fractional Nabla Difference Equations with Constant Coefficients
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 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 informationSolving Integral Equations by Petrov-Galerkin Method and Using Hermite-type 3 1 Elements
5 1 July 24, Antalya, Turkey Dynamical Systems and Applications, Proceedings, pp. 436 442 Solving Integral Equations by Petrov-Galerkin Method and Using Hermite-type 3 1 Elements M. Karami Department of
More informationSUPERCONVERGENCE POINTS OF FRACTIONAL SPECTRAL INTERPOLATION
SIAM J. SCI. COMPUT. Vol. 38, No. 1, pp. A598 A613 c 16 Society for Industrial and Applied Mathematics SUPERCONVERGENCE POINTS OF FRACTIONAL SPECTRAL INTERPOLATION XUAN ZHAO AND ZHIMIN ZHANG Abstract.
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 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 informationThe generalized Euler-Maclaurin formula for the numerical solution of Abel-type integral equations.
The generalized Euler-Maclaurin formula for the numerical solution of Abel-type integral equations. Johannes Tausch Abstract An extension of the Euler-Maclaurin formula to singular integrals was introduced
More 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 informationSOLUTION OF NONLINEAR VOLTERRA-HAMMERSTEIN INTEGRAL EQUATIONS VIA SINGLE-TERM WALSH SERIES METHOD
SOLUTION OF NONLINEAR VOLTERRA-HAMMERSTEIN INTEGRAL EQUATIONS VIA SINGLE-TERM WALSH SERIES METHOD B. SEPEHRIAN AND M. RAZZAGHI Received 5 November 24 Single-term Walsh series are developed to approximate
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 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 informationSPECTRAL PROPERTIES OF THE SIMPLE LAYER POTENTIAL TYPE OPERATORS
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 43, Number 3, 2013 SPECTRAL PROPERTIES OF THE SIMPLE LAYER POTENTIAL TYPE OPERATORS MILUTIN R. OSTANIĆ ABSTRACT. We establish the exact asymptotical behavior
More informationA NEARLY-OPTIMAL ALGORITHM FOR THE FREDHOLM PROBLEM OF THE SECOND KIND OVER A NON-TENSOR PRODUCT SOBOLEV SPACE
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 27, Number 1, Spring 2015 A NEARLY-OPTIMAL ALGORITHM FOR THE FREDHOLM PROBLEM OF THE SECOND KIND OVER A NON-TENSOR PRODUCT SOBOLEV SPACE A.G. WERSCHULZ
More informationTo link to this article:
This article was downloaded by: [Wuhan University] On: 11 March 2015, At: 01:08 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer
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 informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationNEW ESTIMATES FOR RITZ VECTORS
MATHEMATICS OF COMPUTATION Volume 66, Number 219, July 1997, Pages 985 995 S 0025-5718(97)00855-7 NEW ESTIMATES FOR RITZ VECTORS ANDREW V. KNYAZEV Abstract. The following estimate for the Rayleigh Ritz
More informationEXISTENCE AND UNIQUENESS OF SOLUTIONS FOR FOURTH-ORDER BOUNDARY-VALUE PROBLEMS IN BANACH SPACES
Electronic Journal of Differential Equations, Vol. 9(9), No. 33, pp. 1 8. ISSN: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu EXISTENCE AND UNIQUENESS
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 informationPositive solutions for a class of fractional boundary value problems
Nonlinear Analysis: Modelling and Control, Vol. 21, No. 1, 1 17 ISSN 1392-5113 http://dx.doi.org/1.15388/na.216.1.1 Positive solutions for a class of fractional boundary value problems Jiafa Xu a, Zhongli
More informationarxiv: v1 [math.na] 29 Sep 2017
SUPERCONVERGENCE POINTS FOR THE SPECTRAL INTERPOLATION OF RIESZ FRACTIONAL DERIVATIVES BEICHUAN DENG, ZHIMIN ZHANG, AND XUAN ZHAO arxiv:79.3v [math.na] 9 Sep 7 Abstract. In this paper, superconvergence
More informationSolving Integral Equations of the Second Kind by Using Wavelet Basis in the PG Method
5 1 July 24, Antalya, Turkey Dynamical Systems and Applications, Proceedings, pp. 515 52 Solving Integral Equations of the Second Kind by Using Wavelet Basis in the PG Method K. Maleknejad Department of
More informationSUPER-GEOMETRIC CONVERGENCE OF A SPECTRAL ELEMENT METHOD FOR EIGENVALUE PROBLEMS WITH JUMP COEFFICIENTS *
Journal of Comutational Mathematics Vol.8, No.,, 48 48. htt://www.global-sci.org/jcm doi:.48/jcm.9.-m6 SUPER-GEOMETRIC CONVERGENCE OF A SPECTRAL ELEMENT METHOD FOR EIGENVALUE PROBLEMS WITH JUMP COEFFICIENTS
More informationRoots and Coefficients Polynomials Preliminary Maths Extension 1
Preliminary Maths Extension Question If, and are the roots of x 5x x 0, find the following. (d) (e) Question If p, q and r are the roots of x x x 4 0, evaluate the following. pq r pq qr rp p q q r r p
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 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 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 informationRNDr. Petr Tomiczek CSc.
HABILITAČNÍ PRÁCE RNDr. Petr Tomiczek CSc. Plzeň 26 Nonlinear differential equation of second order RNDr. Petr Tomiczek CSc. Department of Mathematics, University of West Bohemia 5 Contents 1 Introduction
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 informationNORMS ON SPACE OF MATRICES
NORMS ON SPACE OF MATRICES. Operator Norms on Space of linear maps Let A be an n n real matrix and x 0 be a vector in R n. We would like to use the Picard iteration method to solve for the following system
More informationSUPERCONVERGENCE OF DG METHOD FOR ONE-DIMENSIONAL SINGULARLY PERTURBED PROBLEMS *1)
Journal of Computational Mathematics, Vol.5, No., 007, 185 00. SUPERCONVERGENCE OF DG METHOD FOR ONE-DIMENSIONAL SINGULARLY PERTURBED PROBLEMS *1) Ziqing Xie (College of Mathematics and Computer Science,
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
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 informationON PROJECTIVE METHODS OF APPROXIMATE SOLUTION OF SINGULAR INTEGRAL EQUATIONS. Introduction Let us consider an operator equation of second kind [1]
GEORGIAN MATHEMATICAL JOURNAL: Vol. 3, No. 5, 1996, 457-474 ON PROJECTIVE METHODS OF APPROXIMATE SOLUTION OF SINGULAR INTEGRAL EQUATIONS A. JISHKARIANI AND G. KHVEDELIDZE Abstract. The estimate for the
More informationResearch Article Soliton Solutions for the Wick-Type Stochastic KP Equation
Abstract and Applied Analysis Volume 212, Article ID 327682, 9 pages doi:1.1155/212/327682 Research Article Soliton Solutions for the Wick-Type Stochastic KP Equation Y. F. Guo, 1, 2 L. M. Ling, 2 and
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationWavelets Application to the Petrov-Galerkin Method for Hammerstein Equations
Wavelets Application to the Petrov-Galerkin Method for Hammerstein Equations Hideaki Kaneko, Richard D. Noren and Boriboon Novaprateep Department of Mathematics and Statistics Old Dominion University Norfolk,
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationA GALERKIN S PERTURBATION TYPE METHOD TO APPROXIMATE A FIXED POINT OF A COMPACT OPERATOR
International Journal of Pure and Applied Mathematics Volume 69 No. 1 2011, 1-14 A GALERKIN S PERTURBATION TYPE METHOD TO APPROXIMATE A FIXED POINT OF A COMPACT OPERATOR L. Grammont Laboratoire de Mathématiques
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 informationUltraconvergence of ZZ Patch Recovery at Mesh Symmetry Points
Ultraconvergence of ZZ Patch Recovery at Mesh Symmetry Points Zhimin Zhang and Runchang Lin Department of Mathematics, Wayne State University Abstract. The ultraconvergence property of the Zienkiewicz-Zhu
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More 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 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 informationSung-Wook Park*, Hyuk Han**, and Se Won Park***
JOURNAL OF THE CHUNGCHEONG MATHEMATICAL SOCIETY Volume 16, No. 1, June 2003 CONTINUITY OF LINEAR OPERATOR INTERTWINING WITH DECOMPOSABLE OPERATORS AND PURE HYPONORMAL OPERATORS Sung-Wook Park*, Hyuk Han**,
More informationPolynomial Preserving Recovery for Quadratic Elements on Anisotropic Meshes
Polynomial Preserving Recovery for Quadratic Elements on Anisotropic Meshes Can Huang, 1 Zhimin Zhang 1, 1 Department of Mathematics, Wayne State University, Detroit, Michigan 480 College of Mathematics
More informationMean Convergence of Interpolation at Zeros of Airy Functions
Mean Convergence of Interpolation at Zeros of Airy Functions D. S. Lubinsky Abstract The classical Erdős-Turán theorem established mean convergence of Lagrange interpolants at zeros of orthogonal polynomials.
More informationA fast method for solving the Heat equation by Layer Potentials
A fast method for solving the Heat equation by Layer Potentials Johannes Tausch Abstract Boundary integral formulations of the heat equation involve time convolutions in addition to surface potentials.
More informationWEYL S THEOREM FOR PAIRS OF COMMUTING HYPONORMAL OPERATORS
WEYL S THEOREM FOR PAIRS OF COMMUTING HYPONORMAL OPERATORS SAMEER CHAVAN AND RAÚL CURTO Abstract. Let T be a pair of commuting hyponormal operators satisfying the so-called quasitriangular property dim
More informationA C 0 linear finite element method for two fourth-order eigenvalue problems
IMA Journal of Numerical Analysis (2017) 37, 2120 2138 doi: 10.1093/imanum/drw051 Advance Access publication on November 2, 2016 A C 0 linear finite element method for two fourth-order eigenvalue problems
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 informationOn compact operators
On compact operators Alen Aleanderian * Abstract In this basic note, we consider some basic properties of compact operators. We also consider the spectrum of compact operators on Hilbert spaces. A basic
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationExistence and multiple solutions for a second-order difference boundary value problem via critical point theory
J. Math. Anal. Appl. 36 (7) 511 5 www.elsevier.com/locate/jmaa Existence and multiple solutions for a second-order difference boundary value problem via critical point theory Haihua Liang a,b,, Peixuan
More informationCONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS
J. Korean Math. Soc. 34 (1997), No. 3, pp. 515 531 CONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS S. K. CHUNG, A.K.PANI AND M. G. PARK ABSTRACT. In this paper,
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 informationarxiv: v1 [math.na] 1 Sep 2018
On the perturbation of an L -orthogonal projection Xuefeng Xu arxiv:18090000v1 [mathna] 1 Sep 018 September 5 018 Abstract The L -orthogonal projection is an important mathematical tool in scientific computing
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 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 informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
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 informationInternational Journal of Pure and Applied Mathematics Volume 55 No ,
International Journal of Pure and Applied Mathematics Volume 55 No. 9, 47-56 THE HANKEL CONVOLUTION OF ARBITRARY ORDER C.K. Li Department of Mathematics and Computer Science Brandon University Brandon,
More informationON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR TERM
Internat. J. Math. & Math. Sci. Vol. 22, No. 3 (999 587 595 S 6-72 9922587-2 Electronic Publishing House ON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR
More informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
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 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 informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationarxiv: v2 [math.pr] 27 Oct 2015
A brief note on the Karhunen-Loève expansion Alen Alexanderian arxiv:1509.07526v2 [math.pr] 27 Oct 2015 October 28, 2015 Abstract We provide a detailed derivation of the Karhunen Loève expansion of a stochastic
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 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 informationFast Structured Spectral Methods
Spectral methods HSS structures Fast algorithms Conclusion Fast Structured Spectral Methods Yingwei Wang Department of Mathematics, Purdue University Joint work with Prof Jie Shen and Prof Jianlin Xia
More informationEnhancing eigenvalue approximation by gradient recovery on adaptive meshes
IMA Journal of Numerical Analysis Advance Access published October 29, 2008 IMA Journal of Numerical Analysis Page 1 of 15 doi:10.1093/imanum/drn050 Enhancing eigenvalue approximation by gradient recovery
More informationA Spectral Method for Nonlinear Elliptic Equations
A Spectral Method for Nonlinear Elliptic Equations Kendall Atkinson, David Chien, Olaf Hansen July 8, 6 Abstract Let Ω be an open, simply connected, and bounded region in R d, d, and assume its boundary
More informationNonlocal problems for the generalized Bagley-Torvik fractional differential equation
Nonlocal problems for the generalized Bagley-Torvik fractional differential equation Svatoslav Staněk Workshop on differential equations Malá Morávka, 28. 5. 212 () s 1 / 32 Overview 1) Introduction 2)
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
More information