Superconvergence Results for the Iterated Discrete Legendre Galerkin Method for Hammerstein Integral Equations
|
|
- Ralph Owen
- 5 years ago
- Views:
Transcription
1 Journal of Computer Science & Computational athematics, Volume 5, Issue, December 05 DOI: /jcscm Superconvergence Results for the Iterated Discrete Legendre Galerkin ethod for Hammerstein Integral Equations Payel Das*, Gnaneshwar Nelakanti Department of athematics, Indian Institute of Technology, Kharagpur -7 30, India * Corresponding author daspayel@gmail.com Abstract: In this paper we analyse the iterated discrete Legendre Galerkin method for Fredholm-Hammerstein integral equation with a smooth kernel. Using a sufficiently accurate numerical quadrature rule, we obtain super-convergence rates for the iterated discrete Legendre Galerkin solutions in both infinity and L -norm. Numerical examples are given to illustrate the theoretical results. Keywords: Hammerstein integral equations, Spectral method, Iterated discrete Galerkin, Numerical quadrature, Superconvergence.. Introduction Let X be a Banach space and consider the following Hammerstein integral equation: x(t) Ω k(t, s)ψ(s, x(s))ds = f(t), t Ω, x X, where Ω is a closed bounded region in R. Hammerstein integral equations arise as a reformulation of various physical phenomena in different branches of study such as mathematical physics, vehicular traffic, biology, economics etc. There has been a notable interest in the numerical analysis of solutions of intergal equations (see [-]). The Galerkin, collocation, Petrov-Galerkin, degenerate kernel and Nystro m methods are the most frequently used projection methods for solving the equations of type (). We are mainly interested in iterated discrete Galerkin method in this paper. The projection methods for solving equation () lead to algebraic nonlinear system, in which the coefficients are integrals, it appeared due to inner products and integral operator. Replacement of these integrals by numerical quadrature rule gives rise to the discrete projection methods. The effect of quadrature error on the convergence rates of the approximate solution is considered in these discrete projection methods (see [], [-]). Discrete projection methods for Fredholm nonlinear integral equations with spline bases and their super-convergence results have been studied by many authors such as Atkinson and Potra [3], Atkinson and Flores [5] and many others. Atkinson and Bogomolny [6] have shown that sufficiently accurate numerical quadrature rules can preserve the rates of convergence of the spline based Galerkin method. However, to get better accuracy in spline based discrete projection () methods, the number of partition points should be increased. Hence in such cases, one has to solve a large system of nonlinear equations, which is computationally very much expensive. To overcome the computational complexities encountered in the existing piecewise polynomial based projection methods, we apply polynomially-based projection methods to nonlinear Fredholm integral equations ([, 5, 7, 8]). We choose the approximating subspaces X n to be global polynomial subspaces of degree n which has dimension n +. The advantage of using global polynomials is that the projection method will imply smaller nonlinear systems, something which is highly desirable in practical computations. In particular here, we choose to use Legendre polynomials, which can be generated recursively with ease and possess nice property of orthogonality. In a recent paper [8], we obtain that the discrete Legendre Galerkin solution of the equation () converges with the order O(n r+ ) in both infinity and L -norm, n being the highest degree of the Legendre polynomial employed in the approximation and r is the smoothness of the kernel k the nonlinear function ψ, the right hand side function f and the solution, with n r. In this paper, we investigate the superconvergence property of the iterated discrete Legendre Galerkin method for Fredholm- Hammerstein integral equation (). Our purpose in this paper is to obtain similar superconvergence results in polynomially-based discrete Galerkin method for Fredholm- Hammerstein integral equation () with a smooth kernel as in the case of piecewise polynomial based discrete Galerkin method. By choosing a numerical quadrature rule of certain degree of precesion, we show that the iterated discrete Legendre Galerkin solution converges with the order O(n r ) in both infinity and L -norm. The organization of this paper is as follows. In Section and 3, we discuss the discrete and iterated discrete Legendre spectral Galerkin methods and obtain superconvergence results. In Section, numerical examples are given to illustrate the theoretical results. Throughout this paper, we assume that c is a generic constant.. Discrete Legendre Galerkin ethod In this section, we describe the discrete Galerkin method for solving Fredholm-Hammerstein integral equations using
2 76 Superconvergence Results for the Iterated Discrete Legendre Galerkin ethod for Hammerstein Integral Equations global polynomial basis functions. Let L [,] be the space of real square-integrable functions on [,] and X = C[,] L [,]. Consider the following Hammerstein integral equation x(t) k(t, s)ψ(s, x(s))ds = f(t), t, where k, f and ψ are known real functions and x is the unknown function to be determined. For a fixed t [,], we denote k t (s) = k(t, s). Throughout the paper, the following assumptions are made on f, k(.,. ) and ψ: I. f C[,], II. k(t, s) C([,] [,]) and = sup t,s [,] k(t, s) <, III. IV. Let lim k t t t(. ) k t (. ) = 0, t, t [,], the nonlinear function ψ(s, x) is continuous in s [,] and is Lipschitz continuous in x, i.e., for any x, x R, c > 0 such that ψ(s, x ) ψ(s, x ) c x x, V. the partial derivative ψ (0,) (s, x) of ψ w.r. to the second variable exists and is Lipschitz continuous in x, i.e., for any x, x R, c > 0 such that ψ (0,) (s, x ) ψ (0,) (s, x ) c x x. Ky(t) = k(t, s)y(s)ds, t [,], y X. Note that, using Holder's inequality, we have for any y X and This implies Ky = sup Ky(t) t [,] = sup k(t, s)y(s)ds t [,] sup t,s [,] k(t, s) y(s) ds () y L, (3) Ky L Ky y L. () K L. (5) We will use Kumar and Sloan s [8] technique for finding the approximate solution of the equation (). To do this, we let Then the solution x of Hammerstein integral equation () is obtained by x = f + Kz. (7) For our convenience, we define a nonlinear operator Ψ: X X by Then (6) takes the form Ψ(x)(t) ψ(t, x(t)). (8) z = Ψ(Kz + f). (9) Let T(u) Ψ(Ku + f), u X, then the equation (9) can be written as z = Tz. Theorem. Let X = C[,], f X and k(.,. ) C([,] [,]) with = sup t,s [,] k(t, s) <. Let ψ(s, y(s)) C([,] R) satisfy the Lipschitz condition in the second variable, i.e. ψ(s, y ) ψ(s, y ) c y y, y, y R, with c <. Then the operator equation z = Tz has a unique solution z 0 X, i.e., we have z 0 = Tz 0. Proof of the above theorem can be easily done using similar technique given in Theorem. of [9]. We denote x 0 be the solution of equation (7) corresponding to the solution z 0 of (9), i.e., x 0 = Kz 0 + f. Now, we describe the discrete Legendre Galerkin method for the solution of Hammerstein integral equation (). To do this, we let X n = span{φ 0, φ, φ,, φ n } be the sequence of Legendre polynomial subspaces of X of degree n, where {φ 0, φ, φ,, φ n } forms an orthonormal basis for X n. Here φ i s are given by φ i (s) = i+ L i (s), i = 0,,, n, (0) where L i s are the Legendre polynomials of degree i. These Legendre polynomials can be generated by the following three-term recurrence relation L 0 (s) =, L (s) = s, and for i =,,, n and s [,] s [,], (i + )L i+ (s) = (i + )sl i (s) il i (s). In practice, the integrals in Galerkin method for solving (7) and (9), appearing due to inner products and the integral operator K cannot be evaluated exactly. The replacement of these integrals by numerical quadrature gives rise to the discrete Galerkin method, which we describe below. First we choose a numerical integration scheme z(t) = ψ(t, x(t)), t [,]. (6)
3 Payel Das, Gnaneshwar Nelakanti 77 where (i) (ii) (n) f(t)dt = w p f(t p ), () p= the weights w p are such that w p > 0, p =,,, (n), () the degree of precision d of the quadrature rule is at least n, that is (n) f(t)dt = w p f(t p ), (3) p= for all polynomials of degree n d. Proof: Let us consider the error function E d (f) = f(s)ds w j f(t j ). Using the quadrature rule defined by ()-(3), we see that E d (p) = 0, p P d, where P d is the space of all polynomials of degree d. For all p P d, consider j= (K n K)z(t) = k(t, s)z(s)ds w j k(t, t j )z(t j ) j= = E d (k(t,. )z(. )) E d (p) For the notational convenience, from now on we set (n) =. Using the above quadrature rule ()-(3) (see Golberg [], Sloan [0]), we define the discrete inner product < f, g > = w p f(t p )g(t p ), f, g C[,]. () p= For the approximation of the integral operator K, using the quadrature rule ()-(3), we consider the Nystro m operator K n defined by (K n z)(t) = w p k(t, t p )z(t p ). (5) p= For the rest of the paper, we assume that z 0 C d [,], k(.,.) C d ([,] [,]), ψ C d ([,] R) and f C d [,], where C d [,] denotes the space of d-times continuously differentiable functions on [-,] and d is the degree of precision of the numerical quadrature rule. For z C d [,], let z d, = max{ z (j) : 0 j d}. We denote where k d, = max{ D (l,p) k : 0 l, p d}, l+p (D (l,p) k)(t, s) = t l sp k(t, s), t, s [,]. In the following theorem, we give the error bounds for the integral operator K and the Nystro m operator K n defined by (5). Theorem. Let k(.,. ) C d ([,] [,]), then for any z C d ([,]), there hold (K n K)z cn d k d, z d,, where c is a constant independent of n and d is the degree of precision of the quadrature rule. = E d [k t (. )z(. ) p]. (6) Since weight functions w j > 0, choosing p(x) =, we j= have ds = w j =. Using this we get E d [k t (. )z(. ) p] = [k(t, s)z(s) p(s)]ds w j [k(t, t j )z(t j ) p(t j )] k t (. )z(. ) p [ k t (. )z(. ) p. j= ds + w j ] As p P d is arbitrary, from the above estimate and using Jackson's Theorem [], we obtain j= E d [k t (. )z(. ) p] inf p P d k t (. )z(. ) p c (k t (. )) d z d, d d. (7) Since from (3), we have d n, it follows that E d [k t (. )z(. ) p] c k d, z d, (n) d, (8) where c is a constant independent of n. Now combining estimates (6) and (8), we get (K n K)z cn d k d, z d,, where c is a constant independent of n. This completes the proof. Note that from Theorem, we see that K n converges pointwise in infinity norm. Hence, K n is pointwise bounded and since X = C[,] is a Banach space with. norm, by Uniform Boundedness Principle we have K n is uniformly bounded, i.e.
4 78 Superconvergence Results for the Iterated Discrete Legendre Galerkin ethod for Hammerstein Integral Equations K n p <, (9) where p is a constant independent of n. Discrete orthogonal projection operator: Discrete orthogonal projection namely hyperinterpolation operator Q n : X X n (see Sloan [0]) is defined by and Q n satisfy n Q n u = < u, φ j > φ j, (0) j=0 < Q n u, φ > = < u, φ >, φ X n. () We quote some crucial properties of Q n from Sloan [0] which are needed for the convergence analysis of the approximate solutions. Lemma. Let Q n : X X n be the hyperinterpolation operator defined by (0). Then for any u X, the following result holds < u Q n u, u Q n u > = min χ X n < u χ, u χ >. Lemma. Let Q n : X X n be the hyperinterpolation operator defined by (0). Then the following results hold i. Q n u L p u, where p is a constant ii. iii. independent of n. There exists a constant c > 0 such that for any n N and u X, Q n u u L c inf u φ L 0, as n. φ X n () In particular, if u C r [,], then there holds, Q n u u L cn r u r,, n r, (3) where c is a constant independent u and n. Now the discrete Legendre Galerkin method for (9) is to find z n X n, such that z n = Q n ѱ(k n z n + f). () Letting T n (u) Q n ѱ(k n u + f), u X, the equation () can be written as z n = T n z n. The corresponding discrete approximate solution x n of x is defined by x n = K n z n + f. (5) We quote the following results from [8], which gives the convergence rates of the discrete Legendre Galerkin solutions in both infinity and L -norm. Theorem 3. Let z 0 C r [,], n r, be an isolated solution of the equation (9). Let Q n be the discrete orthogonal projection operator defined by (0). Assume that is not an eigenvalue of the linear operator T (z 0 ), then for sufficiently n, the approximate solution z n defined by the equation () is the unique solution in the sphere B(z 0, δ) = {z: z z 0 δ} for some δ > 0. oreover, there exists a constant 0 < q <, independent of n such that α n + q z n z 0 α n q, where α n = (I T n (z 0 )) (T n (z 0 ) T(z 0 )), and z n z 0 = O(n r+ ). Theorem. Let z 0 C r [,] be an isolated solution of the equation (9) and x 0 be a isolated solution of the equation (7) such that x 0 = Kz 0 + f. x n be the discrete Legendre Galerkin approximation of x 0. Then the following holds x 0 x n, x 0 x n L = O(n r+ ). 3. Iterated Discrete Legendre Galerkin ethod In order to obtain more accurate approximation solution, we further consider the iterated discrete Legendre Galerkin approximation of z 0 in this section. To this end, we define the iterated discrete solution as z n = ѱ(k n z n + f). (6) Applying Q n on both sides of the equation (6), we obtain Q n z n = Q n ѱ(k n z n + f). (7) From equations () and (7), it follows that Q n z n = z n. Using this, we see that the iterated solution z n satisfies the following equation z n = ѱ(k n Q n z n + f). (8) Letting T n(u) ѱ(k n Q n u + f), u X, the equation (8) can be written as z n = T nz n. Corresponding approximate solution x n of x is given by x n = K n z n + f. (9) Next we discuss the existence of the iterated approximate solutions and their error bounds. To do this, we first recall the following definition of υ-convergence and a lemma from []. Definition. Let X be Banach space and BL(X) be space of bounded linear operators from X into X. Let F n, F BL(X). We say F n is υ-convergent to F if F n c <, (F n F)F 0, (F n F)F n 0, as n. Lemma 3. (Ahues et al. []) Let X be a Banach space and F n, F be bounded linear operators on X. If F n is υ-convergent to F and (I F) exists, then (I F n ) exists and uniformly bounded on X, for sufficiently large n. Lemma. Let k(.,. ) C d ([,] [,]), d n > n r, then for any u X, the following hold
5 Payel Das, Gnaneshwar Nelakanti 79 (K n Q n K)u 0, as n. In particular if u C r [,], then (K n Q n K)u = O(n r ). Proof: For any u n X n, it follows that p= ) ( w p [u n (t p )] We denote = ( [u n (s)] e(k t (. )) = (K K n )u n (t) ds) = u n L. (30) = k(t, s)u n (s)ds w p k(t, t p )u n (t p ). p= Using Cauchy-Schwartz inequality and estimate (30), we get e(k t (. )) = k(t, s)u n (s)ds w p k(t, t p )u n (t p ) = ( k(t, s) ds p= p= ) ( un (s) ds) + ( w p k(t, t p ) ) ( w p u n (t p ) ) ( + ( w p ) p= k u n L. ) p= Since for any y X n, e(y) = 0, we have k u n L e(k t (. )) = e(k t (. ) y) inf y X n k y u n L. Using this and Jackson's Theorem [], we obtain (K n K)u n cn r k r, u n L. (3) Since Q n u X n, u X, from (3) we obtain (K K n )Q n u cn r k r, Q n u L cn r k r, p u. (3) Since d r, using estimates (3), (3) and (3), we have for any u C r [,], (K n Q n K)u (K K n )Q n u + K(Q n I)u = O(n r ). (3) Note that form estimates () and (33), it follows that for any u X, (K n Q n K)u 0, as n. (35) This completes the proof. Theorem 5. Let z 0 C d [,], d n > n r, be an isolated solution of the equation (9). Assume that is not an eigenvalue of T (z 0 ). Then for sufficiently large n, the operator I T n (z 0 ) is invertible on X and there exist a constant L > 0 independent of n such that (I T n (z 0 )) L. Proof: Consider We have T n (z 0 ) = ѱ (K n Q n z 0 + f)k n Q n. (36) ѱ (K n Q n z 0 + f) ѱ (K n Q n z 0 + f) ѱ (K n z 0 + f) + ѱ (K n z 0 + f). (37) Form Lemma and Jackson's Theorem [], we have for any z C r [,], < (Q n I)z, (Q n I)z > = min < z χ, z χ > χ X n = min { w p (z χ) (t p )} χ X n p= ( w p ) p= inf χ X n z χ cn r z r,. (38) Using Lipschitz continuity of ψ (0,) (.,. ) and (38), we obtain ѱ (K n Q n z 0 + f) ѱ (K n z 0 + f) c K n (Q n I)z 0 = c sup K n (Q n I)z 0 (t) t [,] = c sup w p k(t, t p )(Q n I)z 0 (t p ) t [,] c sup t [,] p= p= ( w p k(t, t p ) ) ( w p [(Q n I)z 0 (t p )] ) p= cn r k r, p u + (Q n I)u L (33)
6 80 Superconvergence Results for the Iterated Discrete Legendre Galerkin ethod for Hammerstein Integral Equations c ( w p ) p= k < (Q n I)z 0, (Q n I)z 0 > c cn r k z 0 r, 0, as n. (39) Using Theorem, Lipschitz continuity of ψ (0,) (.,. ) and boundedness of ψ (Kz 0 + f), we have ѱ (K n z 0 + f) ѱ (K n z 0 + f) ѱ (Kz 0 + f) + ѱ (Kz 0 + f) c (K n K)z 0 + ѱ (Kz 0 + f) c cn d k d, z 0 d, + ѱ (Kz 0 + f) <. (0) From estimates (37), (39) and (0), it follows that ѱ (K n Q n z 0 + f) B <, () where B is a constant independent of n. Since Q n z X n, from estimate (30), we have p= ) ( w p [Q n z(t p )] = Q n z L. () Hence using Cauchy-Schwartz inequality and the estimate (), we have This implies K n Q n z = sup K n Q n z (t) t [,] = sup w p k(t, t p )Q n z(t p ) t [,] p= p= ( w p ) ( w p [Q n z(t p )] ) p= = Q n z L p z <. (3) K n Q n p. () Combining estimates (36), () and (), we have T n (z 0 ) B <, (5) where B is independent of n. Hence it follows that T n (z 0 ) is uniformly bounded. Using estimates (3), () and (5), we have [T n (z 0 ) T (z 0 )]T n (z 0 ) = [ѱ (K n Q n z 0 + f)k n Q n ѱ (Kz 0 + f)k]t n (z 0 ) ѱ (K n Q n z 0 + f)(k n Q n K)T n (z 0 ) ѱ (K n Q n z 0 + f) (K n Q n K)T n (z 0 ) + ѱ (K n Q n z 0 + f) ѱ (Kz 0 + f) KT n (z 0 ) B (K n Q n K)T n (z 0 ) + B ѱ (K n Q n z 0 + f) ѱ (Kz 0 + f). (6) Using Lipschitz continuity of ψ (0,) (.,. ) and Lemma, we obtain ѱ (K n Q n z 0 + f) ѱ (Kz 0 + f) c (K n Q n K)z 0 ) 0, as n. (7) Next, Let B {x X: x } be the closed unit ball in X. We have T n (z 0 ) = ѱ (K n Q n z 0 + f)k n Q n. Since {K n Q n } is a sequence of compact operators and ѱ (K n Q n z 0 + f) is uniformly bounded, T n (z 0 ) are compact operators. Thus S = {T n (z 0 )x x B, n N} is relatively compact set. Using estimate (35), we can conclude (K n Q n K)T n (z 0 ) = sup{ (K n Q n K)T n (z 0 )x : x B } = sup{ (K n Q n K)y : y S}, 0, as n. (8) Combining the estimates (6), (7) and (8), we have [T n (z 0 ) T (z 0 )]T n (z 0 ) 0, as n. Following the similar steps and using the fact that T (z 0 ) is compact, it can be proved that [T n (z 0 ) T (z 0 )]T (z 0 ) 0, as n. Hence T n (z 0 ) is υ-convergent to T (z 0 ) in. norm. By direct application of Lemma 3, it follows that for sufficiently large n, I T n (z 0 ) is invertible, i.e., there exist a constant L > 0, independent of n such that (I T n (z 0 )) L. This completes the proof. Theorem 6. Let z 0 C d [,], d n > n r, be an isolated solution of the equation (9). Let Q n : X X n be the discrete orthogonal projection operator defined by (0). Assume that is not an eigenvalue of T (z 0 ), then for sufficiently large n, the iterated solution z n defined by (8) is the unique solution in the sphere B(z 0, δ) = {z: z z 0 δ}. oreover, there exists a constant 0 < q <, independent of n such that β n + q z n z 0 β n q, where β n = (I T n (z 0 )) (T n(z 0 ) T(z 0 )). Proof: From Theorem 5, we have (I T n (z 0 )) exists and it is uniformly bounded in infinity norm, i.e., L > 0 such that (I T n (z 0 )) L. + [ѱ (K n Q n z 0 + f) ѱ (Kz 0 + f)]kt n (z 0 )
7 Payel Das, Gnaneshwar Nelakanti 8 Using Lipschitz continuity of ψ (0,) (.,. ) and estimate (3), we have for any z B(z 0, δ), T n(z 0 ) T(z 0 ) = ѱ(k n Q n z 0 + f) ѱ(kz 0 + f) c (K n Q n K)z 0 [T n (z) T n (z 0 )]v = [ѱ (K n Q n z + f) ѱ (K n Q n z 0 + f)]k n Q n v ѱ (K n Q n z + f) ѱ (K n Q n z 0 + f) K n Q n v We have c [ (K n Q n K n )z 0 + (K n K)z 0 ]. (50) This implies c K n Q n (z 0 z) K n Q n v c p z z 0 v c p δ v. sup (I T n (z 0 )) (T n (z) T n (z 0 )) Lc p δ z z 0 δ q, (say), where we choose δ in such a way that q (0,). This proves the estimate (.) of Theorem in [3]. Now using the Lipschtiz continuity of ψ(., x(. )) and Lemma, we have T n(z 0 ) T(z 0 ) = ѱ(k n Q n z 0 + f) ѱ(kz 0 + f) Hence c (K n Q n K)z 0 0, as n. β n = (I T n (z 0 )) (T n(z 0 ) T(z 0 )) Lc (K n Q n K)z 0 0, as n. Choose n large enough such that β n δ( q). We choose δ in such a way that q (0,). This proves the estimate (.5) of Theorem in [3]. Hence applying Theorem of [3], we obtain β n + q z n z 0 β n q, where β n = (I T n (z 0 )) (T n(z 0 ) T(z 0 )). This completes the proof. Theorem 7. Let z 0 C d [,], d n > n r, be an isolated solution of the equation (9) and x 0 be a isolated solution of the equation (7) such that x 0 = Kz 0 + f. x n be the iterated discrete Legendre Galerkin approximation of x 0. Then the following hold x 0 x n L, x 0 x n = O(n r ). Proof: From Theorem 6, we have β n z n z 0 q c (I T n (z 0 )) (T n(z 0 ) T(z 0 )). (9) (K n Q n K n )z 0 = sup (K n Q n K n )z 0 (t) t [,] = sup w p k(t, t p )(Q n I) z 0 (t p ) t [,] p= = sup < k t (. ), (Q n I)z 0 >. (5) t [,] Using the estimate (38), orthogonality property of Q n and Cauchy-Schwartz inequality, we obtain < k t (. ), (Q n I)z 0 > = < (Q n I)k t (. ), (Q n I)z 0 > = w p (Q n I)k(t, t p )(Q n I)z 0 (t p ) p= ( w p [(Q n I)k(t, t p )] p= ) ( w p [(Q n I)z 0 (t p )] ) p= = < (Q n I)k t (. ), (Q n I)k t (. ) > < (Q n I)z 0, (Q n I)z 0 > cn r k r, z 0 r,. (5) Since d r, using Theorem and estimates (9), (50), (5) and (5), we have z n z 0 = O(n min{r,d} )= O(n r ). (53) From estimates (7), (9) and (9), we get x n x 0 = K n z n Kz 0 K n (z n z 0 ) + (K n K)z 0 p (z n z 0 ) + (K n K)z 0. Hence using Theorem and estimate (53) and the fact that d r, we obtain x 0 x n = O(n min{r,d} ) = O(n r ). Also as x 0 x n L x 0 x n, we get x 0 x n L = O(n r ). Now using Lipschtiz continuity of ψ(.,. ), we have This completes the proof.
8 8 Superconvergence Results for the Iterated Discrete Legendre Galerkin ethod for Hammerstein Integral Equations. Conclusions From Theorem, we see that the discrete Legendre Galerkin solution of Hammerstein integral equation converges with the optimal order, O(n r+ ) in both L and infinity norm. Result of Theorem 7 imply that, if the numerical quadrature is of certain degree of precision, the iterated discrete Legendre Galerkin solution converges with the order O(n r ) in both L and infinity norm. This shows that the iterated discrete Legendre Galerkin approximation exhibits superconvergence. 5. Numerical Example In this section, we present the numerical results. For that we take the Legendre polynomials as the basis functions for the subspace X n. We present the errors of the approximate solutions under the discrete Legendre Galerkin and iterated discrete Legendre Galerkin methods in both infinity and L - norm. For computations we use Gauss quadrature rule and the Newton-Raphson method to solve the nonlinear systems. The numerical algorithms are compiled by using atlab. In the following tables, n represents the highest degree of the Legendre polynomials employed in the computation. Example. We consider the following integral equation x(t) k(t, s)ψ (s, x(s))ds = f(t), t, with the kernel function k(t, s) = 3π 6 ψ(s, x(s)) = x(s) and the function f(t) = exact solution is given by (t) = cos ( πt ). Table cos (π s t ), ). The cos (πt n x 0 x n L x 0 x n x 0 x n L x 0 x n e-0 0.0e e e e e-0 0.7e-0 0.0e e e e-0 0.7e e e e- 0.e e e e e- Example We consider the following integral equation x(t) k(t, s)ψ (s, x(s))ds = f(t), t, with the kernel function k(t, s) = cos(πt) sin(πs), ψ(s, x(s)) = x(s) 3 and the function f(t) = sin (πt). The exact solution is given by x(t) = sin(πt) + (0 39)cos (πt). 3 From Tables and, we see that the iterated discrete Legendre Galerkin approximation converges much faster than the discrete Legendre Galerkin solution. 0 Table n x 0 x n L x 0 x n x 0 x n L x 0 x n e-0 0.e e e-07 0.e-09 References 0.556e e-0 0.9e e e e e e e e- 0.9e e e e- 0.78e-3 [] K. E. Atkinson, F. A. Potra, F. A, The discrete Galerkin method for nonlinear integral equations, Journal of Integral Equations and Applications, vol., no, pp. 7-5, 988. [] K. E. Atkinson, A survey of numerical methods for solving nonlinear integral equations, Journal of Integral Equations and Applications, vol., no., pp. 5-6, 99. [3] Z. Chen, G. Long, G. Nelakanti, Y. Zhang, Iterated Fast Collocation ethods for Integral Equations of the Second Kind, Journal of Scientific Computing, vol. 57, no. 3, pp , 03. [] P. Das, G. Nelakanti, Convergence Analysis of Discrete Legendre Spectral Projection ethods for Hammerstein Integral Equations of ixed Type, Applied athematics and Computation, vol. 65, pp , 05. [5] P. Das,.. Sahani, G. Nelakanti. Convergence analysis of Legendre spectral projection methods for Hammerstein integral equations of mixed type, Journal of Applied athematics and Computing, vol. 9, no. -, pp , 0. [6] G. Ben-yu, Spectral methods and their applications World Scientific, 998. [7] H. Kaneko, Y. Xu, Superconvergence of the iterated Galerkin methods for Hammerstein equations, SIA Journal on Numerical Analysis, vol. 33, no. 3, pp , 996. [8] S. Kumar, I. H. Sloan, A new collocation-type method for Hammerstein equations, athematics of Computation, vol. 8, no. 78, pp , 987. [9] S. Kumar, Superconvergence of a collocation-type method for Hammerstein equations, IA Journal of Numerical Analysis, vol. 7, no. 3, pp , 987. [0] S. Kumar, A discrete collocation-type method for Hammerstein equations, SIA Journal on Numerical Analysis, vol. 5, no., pp. 38-3, 988. [] G. Long, G. Nelakanti, X. Zhang, Iterated fast multiscale Galerkin methods for Fredholm integral equations of second kind with weakly singular kernels, Applied Numerical athematics, vol. 6, no.3, pp. 0-, 0. [] G. ichael, Improved convergence rates for some discrete Galerkin methods, Journal of Integral Equations and Applications, vol. 8, no. 3, pp , 996.
9 Payel Das, Gnaneshwar Nelakanti 83 [3]. A. Golberg, C. S. Chen, Discrete projection methods for integral equations, Computational echanics Publications, Southampton, 997. []. Golberg, H. Bowman., Optimal convergence rates for some discrete projection methods, Applied mathematics and computation, vol. 96, no., pp. 37-7, 998. [5] K. Atkinson, J. Flores, The discrete collocation method for nonlinear integral equations, IA Journal of Numerical Analysis, vol. 3, no.., pp. 95-3, 993. [6] K. E. Atkinson, A. Bogomolny. The discrete Galerkin method for integral equations, athematics of computation, vol. 8, no. 78, pp , 987. [7] P. Das,.. Sahani, G. Nelakanti, Legendre spectral projection methods for Urysohn integral equations, Journal of Computational and Applied athematics, vol. 63, pp. 88-0, 0. [8] P. Das, G. Nelakanti, G. Long., Discrete Legendre spectral projection methods for Fredholm-Hammerstein integral equations, Journal of Computational and Applied athematics, vol. 78, pp , 05. [9] H. Kaneko, R. D. Noren, Y. Xu Regularity of the solution of Hammerstein equations with weakly singular kernels, Integral Equations and Operator Theory, vol. 3, no. 5, pp , 990. [0] I. H. Sloan, Polynomial interpolation and hyperinterpolation over general regions, Journal of Approximation Theory, vol. 83, no., pp. 38-5, 995. [] L. L. Schumaker, Spline functions: basis theory, John Wiley and Sons, New York, 98. []. Ahues, A. Largillier, B. Limaye, Spectral computations for bounded operators, CRC Press, 00. [3] G.. Vainikko, Galerkin's perturbation method and the general theory of approximate methods for nonlinear equations, USSR Computational athematics and athematical Physics, vol. 7, no., pp. -, 967.
Wavelets 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 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 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 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 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 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 Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
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 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 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 informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationFunctional Analysis. Martin Brokate. 1 Normed Spaces 2. 2 Hilbert Spaces The Principle of Uniform Boundedness 32
Functional Analysis Martin Brokate Contents 1 Normed Spaces 2 2 Hilbert Spaces 2 3 The Principle of Uniform Boundedness 32 4 Extension, Reflexivity, Separation 37 5 Compact subsets of C and L p 46 6 Weak
More informationINTEGRAL equations are often matched with mathematical
On Taylor Expansion Methods for Multivariate Integral Equations of the Second Kind Boriboon Novaprateep, Khomsan Neamprem, and Hideaki Kaneko Abstract A new Taylor series method that the authors originally
More informationSolving Fredholm integral equations with application of the four Chebyshev polynomials
Journal of Approximation Theory and Applied Mathematics, 204 Vol. 4 Solving Fredholm integral equations with application of the four Chebyshev polynomials Mostefa NADIR Department of Mathematics University
More informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationA 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 informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationMTH 309Y 37. Inner product spaces. = a 1 b 1 + a 2 b a n b n
MTH 39Y 37. Inner product spaces Recall: ) The dot product in R n : a. a n b. b n = a b + a 2 b 2 +...a n b n 2) Properties of the dot product: a) u v = v u b) (u + v) w = u w + v w c) (cu) v = c(u v)
More informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
More informationResearch Article Approximation of Solutions of Nonlinear Integral Equations of Hammerstein Type with Lipschitz and Bounded Nonlinear Operators
International Scholarly Research Network ISRN Applied Mathematics Volume 2012, Article ID 963802, 15 pages doi:10.5402/2012/963802 Research Article Approximation of Solutions of Nonlinear Integral Equations
More informationl(y j ) = 0 for all y j (1)
Problem 1. The closed linear span of a subset {y j } of a normed vector space is defined as the intersection of all closed subspaces containing all y j and thus the smallest such subspace. 1 Show that
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 information16 1 Basic Facts from Functional Analysis and Banach Lattices
16 1 Basic Facts from Functional Analysis and Banach Lattices 1.2.3 Banach Steinhaus Theorem Another fundamental theorem of functional analysis is the Banach Steinhaus theorem, or the Uniform Boundedness
More informationAnalysis Preliminary Exam Workshop: Hilbert Spaces
Analysis Preliminary Exam Workshop: Hilbert Spaces 1. Hilbert spaces A Hilbert space H is a complete real or complex inner product space. Consider complex Hilbert spaces for definiteness. If (, ) : H H
More informationBIBLIOGRAPHY KENDALL E ATKINSON
September 28, 2017 BIBLIOGRAPHY KENDALL E ATKINSON I. BOOKS 1. K. Atkinson, D. Chien, and O. Hansen. Spectral Methods Using Multivariate Polynomials On The Unit Ball, submitted for publication. 2. K. Atkinson
More information2. Dual space is essential for the concept of gradient which, in turn, leads to the variational analysis of Lagrange multipliers.
Chapter 3 Duality in Banach Space Modern optimization theory largely centers around the interplay of a normed vector space and its corresponding dual. The notion of duality is important for the following
More informationTHE SKOROKHOD OBLIQUE REFLECTION PROBLEM IN A CONVEX POLYHEDRON
GEORGIAN MATHEMATICAL JOURNAL: Vol. 3, No. 2, 1996, 153-176 THE SKOROKHOD OBLIQUE REFLECTION PROBLEM IN A CONVEX POLYHEDRON M. SHASHIASHVILI Abstract. The Skorokhod oblique reflection problem is studied
More informationSPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS
SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS G. RAMESH Contents Introduction 1 1. Bounded Operators 1 1.3. Examples 3 2. Compact Operators 5 2.1. Properties 6 3. The Spectral Theorem 9 3.3. Self-adjoint
More information10 The Finite Element Method for a Parabolic Problem
1 The Finite Element Method for a Parabolic Problem In this chapter we consider the approximation of solutions of the model heat equation in two space dimensions by means of Galerkin s method, using piecewise
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationMATH 5640: Fourier Series
MATH 564: Fourier Series Hung Phan, UMass Lowell September, 8 Power Series A power series in the variable x is a series of the form a + a x + a x + = where the coefficients a, a,... are real or complex
More informationFredholm Theory. April 25, 2018
Fredholm Theory April 25, 208 Roughly speaking, Fredholm theory consists of the study of operators of the form I + A where A is compact. From this point on, we will also refer to I + A as Fredholm operators.
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More informationBrownian Motion and Conditional Probability
Math 561: Theory of Probability (Spring 2018) Week 10 Brownian Motion and Conditional Probability 10.1 Standard Brownian Motion (SBM) Brownian motion is a stochastic process with both practical and theoretical
More informationBohl exponent for time-varying linear differential-algebraic equations
Bohl exponent for time-varying linear differential-algebraic equations Thomas Berger Institute of Mathematics, Ilmenau University of Technology, Weimarer Straße 25, 98693 Ilmenau, Germany, thomas.berger@tu-ilmenau.de.
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 informationMA5206 Homework 4. Group 4. April 26, ϕ 1 = 1, ϕ n (x) = 1 n 2 ϕ 1(n 2 x). = 1 and h n C 0. For any ξ ( 1 n, 2 n 2 ), n 3, h n (t) ξ t dt
MA526 Homework 4 Group 4 April 26, 26 Qn 6.2 Show that H is not bounded as a map: L L. Deduce from this that H is not bounded as a map L L. Let {ϕ n } be an approximation of the identity s.t. ϕ C, sptϕ
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 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 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 informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
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 informationORTHOGONAL SERIES REGRESSION ESTIMATORS FOR AN IRREGULARLY SPACED DESIGN
APPLICATIONES MATHEMATICAE 7,3(000), pp. 309 318 W.POPIŃSKI(Warszawa) ORTHOGONAL SERIES REGRESSION ESTIMATORS FOR AN IRREGULARLY SPACED DESIGN Abstract. Nonparametric orthogonal series regression function
More informationEcon Lecture 14. Outline
Econ 204 2010 Lecture 14 Outline 1. Differential Equations and Solutions 2. Existence and Uniqueness of Solutions 3. Autonomous Differential Equations 4. Complex Exponentials 5. Linear Differential Equations
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 informationCOMPACT OPERATORS. 1. Definitions
COMPACT OPERATORS. Definitions S:defi An operator M : X Y, X, Y Banach, is compact if M(B X (0, )) is relatively compact, i.e. it has compact closure. We denote { E:kk (.) K(X, Y ) = M L(X, Y ), M compact
More informationContinuous Functions on Metric Spaces
Continuous Functions on Metric Spaces Math 201A, Fall 2016 1 Continuous functions Definition 1. Let (X, d X ) and (Y, d Y ) be metric spaces. A function f : X Y is continuous at a X if for every ɛ > 0
More informationINTEGRABLE SOLUTIONS OF A FUNCTIONAL-INTEGRAL EQUATION. In this paper we consider the following functional-integral equation. k(t, s)g(s, y(s)) ds
JOURNAL OF INTEGRAL EQUATIONS AN APPLICATIONS Volume 4, Number 1, Winter 1992 INTEGRABLE SOLUTIONS OF A FUNCTIONAL-INTEGRAL EQUATION G. EMMANUELE ABSTRACT. We consider a very general functional-integral
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 informationHilbert space methods for quantum mechanics. S. Richard
Hilbert space methods for quantum mechanics S. Richard Spring Semester 2016 2 Contents 1 Hilbert space and bounded linear operators 5 1.1 Hilbert space................................ 5 1.2 Vector-valued
More information4 Hilbert spaces. The proof of the Hilbert basis theorem is not mathematics, it is theology. Camille Jordan
The proof of the Hilbert basis theorem is not mathematics, it is theology. Camille Jordan Wir müssen wissen, wir werden wissen. David Hilbert We now continue to study a special class of Banach spaces,
More informationShape Preserving Approximation: the Final Frontier???
Shape Preserving Approximation: the Final Frontier??? Kirill Kopotun kopotunk@ccumanitobaca Department of Mathematics and the Institute of Industrial Mathematical Sciences, University of Manitoba, Winnipeg,
More informationI. The space C(K) Let K be a compact metric space, with metric d K. Let B(K) be the space of real valued bounded functions on K with the sup-norm
I. The space C(K) Let K be a compact metric space, with metric d K. Let B(K) be the space of real valued bounded functions on K with the sup-norm Proposition : B(K) is complete. f = sup f(x) x K Proof.
More informationNUMERICAL METHOD FOR THE MIXED VOLTERRA-FREDHOLM INTEGRAL EQUATIONS USING HYBRID LEGENDRE FUNCTIONS
Conference Applications of Mathematics 215, in honor of the birthday anniversaries of Ivo Babuška (9), Milan Práger (85), and Emil Vitásek (85) J. Brandts, S. Korotov, M. Křížek, K. Segeth, J. Šístek,
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationA G Ramm, Implicit Function Theorem via the DSM, Nonlinear Analysis: Theory, Methods and Appl., 72, N3-4, (2010),
A G Ramm, Implicit Function Theorem via the DSM, Nonlinear Analysis: Theory, Methods and Appl., 72, N3-4, (21), 1916-1921. 1 Implicit Function Theorem via the DSM A G Ramm Department of Mathematics Kansas
More informationReview and problem list for Applied Math I
Review and problem list for Applied Math I (This is a first version of a serious review sheet; it may contain errors and it certainly omits a number of topic which were covered in the course. Let me know
More informationFinal Year M.Sc., Degree Examinations
QP CODE 569 Page No Final Year MSc, Degree Examinations September / October 5 (Directorate of Distance Education) MATHEMATICS Paper PM 5: DPB 5: COMPLEX ANALYSIS Time: 3hrs] [Max Marks: 7/8 Instructions
More informationMain topics and some repetition exercises for the course MMG511/MVE161 ODE and mathematical modeling in year Main topics in the course:
Main topics and some repetition exercises for the course MMG5/MVE6 ODE and mathematical modeling in year 04. Main topics in the course:. Banach fixed point principle. Picard- Lindelöf theorem. Lipschitz
More informationA Banach Gelfand Triple Framework for Regularization and App
A Banach Gelfand Triple Framework for Regularization and Hans G. Feichtinger 1 hans.feichtinger@univie.ac.at December 5, 2008 1 Work partially supported by EUCETIFA and MOHAWI Hans G. Feichtinger hans.feichtinger@univie.ac.at
More informationQuadrature approaches to the solution of two point boundary value problems
Quadrature approaches to the solution of two point boundary value problems Seth F. Oppenheimer Mohsen Razzaghi Department of Mathematics and Statistics Mississippi State University Drawer MA MSU, MS 39763
More informationThe Navier-Stokes Equations with Time Delay. Werner Varnhorn. Faculty of Mathematics University of Kassel, Germany
The Navier-Stokes Equations with Time Delay Werner Varnhorn Faculty of Mathematics University of Kassel, Germany AMS: 35 (A 35, D 5, K 55, Q 1), 65 M 1, 76 D 5 Abstract In the present paper we use a time
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationNonlinear Systems Theory
Nonlinear Systems Theory Matthew M. Peet Arizona State University Lecture 2: Nonlinear Systems Theory Overview Our next goal is to extend LMI s and optimization to nonlinear systems analysis. Today we
More information5 Compact linear operators
5 Compact linear operators One of the most important results of Linear Algebra is that for every selfadjoint linear map A on a finite-dimensional space, there exists a basis consisting of eigenvectors.
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 informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationRegularity and Decay Estimates of the Navier-Stokes Equations
Regularity and Decay Estimates of the Navier-Stokes Equations Hantaek Bae Ulsan National Institute of Science and Technology (UNIST), Korea Recent Advances in Hydrodynamics, 216.6.9 Joint work with Eitan
More informationOn the Ψ - Exponential Asymptotic Stability of Nonlinear Lyapunov Matrix Differential Equations
DOI: 1.2478/awutm-213-12 Analele Universităţii de Vest, Timişoara Seria Matematică Informatică LI, 2, (213), 7 28 On the Ψ - Exponential Asymptotic Stability of Nonlinear Lyapunov Matrix Differential Equations
More informationS chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1.
Sep. 1 9 Intuitively, the solution u to the Poisson equation S chauder Theory u = f 1 should have better regularity than the right hand side f. In particular one expects u to be twice more differentiable
More informationHilbert spaces. Let H be a complex vector space. Scalar product, on H : For all x, y, z H and α, β C α x + β y, z = α x, z + β y, z x, y = y, x
Hilbert spaces Let H be a complex vector space. Scalar product, on H : For all x, y, z H and α, β C α x + β y, z = α x, z + β y, z x, y = y, x x, x 0 x, x = 0 x = 0 Orthogonality: x y x, y = 0 A scalar
More informationAn Approximate Solution for Volterra Integral Equations of the Second Kind in Space with Weight Function
International Journal of Mathematical Analysis Vol. 11 17 no. 18 849-861 HIKARI Ltd www.m-hikari.com https://doi.org/1.1988/ijma.17.771 An Approximate Solution for Volterra Integral Equations of the Second
More informationAn introduction to some aspects of functional analysis
An introduction to some aspects of functional analysis Stephen Semmes Rice University Abstract These informal notes deal with some very basic objects in functional analysis, including norms and seminorms
More information2) Let X be a compact space. Prove that the space C(X) of continuous real-valued functions is a complete metric space.
University of Bergen General Functional Analysis Problems with solutions 6 ) Prove that is unique in any normed space. Solution of ) Let us suppose that there are 2 zeros and 2. Then = + 2 = 2 + = 2. 2)
More informationNonlinear equations. Norms for R n. Convergence orders for iterative methods
Nonlinear equations Norms for R n Assume that X is a vector space. A norm is a mapping X R with x such that for all x, y X, α R x = = x = αx = α x x + y x + y We define the following norms on the vector
More informationElements of Positive Definite Kernel and Reproducing Kernel Hilbert Space
Elements of Positive Definite Kernel and Reproducing Kernel Hilbert Space Statistical Inference with Reproducing Kernel Hilbert Space Kenji Fukumizu Institute of Statistical Mathematics, ROIS Department
More informationReal Analysis Notes. Thomas Goller
Real Analysis Notes Thomas Goller September 4, 2011 Contents 1 Abstract Measure Spaces 2 1.1 Basic Definitions........................... 2 1.2 Measurable Functions........................ 2 1.3 Integration..............................
More informationFUNCTIONAL ANALYSIS LECTURE NOTES: WEAK AND WEAK* CONVERGENCE
FUNCTIONAL ANALYSIS LECTURE NOTES: WEAK AND WEAK* CONVERGENCE CHRISTOPHER HEIL 1. Weak and Weak* Convergence of Vectors Definition 1.1. Let X be a normed linear space, and let x n, x X. a. We say that
More informationChapter 7: Bounded Operators in Hilbert Spaces
Chapter 7: Bounded Operators in Hilbert Spaces I-Liang Chern Department of Applied Mathematics National Chiao Tung University and Department of Mathematics National Taiwan University Fall, 2013 1 / 84
More informationViscosity Iterative Approximating the Common Fixed Points of Non-expansive Semigroups in Banach Spaces
Viscosity Iterative Approximating the Common Fixed Points of Non-expansive Semigroups in Banach Spaces YUAN-HENG WANG Zhejiang Normal University Department of Mathematics Yingbing Road 688, 321004 Jinhua
More informationA Spectral Method for Elliptic Equations: The Neumann Problem
A Spectral Method for Elliptic Equations: The Neumann Problem Kendall Atkinson Departments of Mathematics & Computer Science The University of Iowa Olaf Hansen, David Chien Department of Mathematics California
More 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 informationAnalysis Comprehensive Exam Questions Fall 2008
Analysis Comprehensive xam Questions Fall 28. (a) Let R be measurable with finite Lebesgue measure. Suppose that {f n } n N is a bounded sequence in L 2 () and there exists a function f such that f n (x)
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationA brief introduction to ordinary differential equations
Chapter 1 A brief introduction to ordinary differential equations 1.1 Introduction An ordinary differential equation (ode) is an equation that relates a function of one variable, y(t), with its derivative(s)
More informationFunctional Analysis Exercise Class
Functional Analysis Exercise Class Week 2 November 6 November Deadline to hand in the homeworks: your exercise class on week 9 November 13 November Exercises (1) Let X be the following space of piecewise
More informationAnalysis in weighted spaces : preliminary version
Analysis in weighted spaces : preliminary version Frank Pacard To cite this version: Frank Pacard. Analysis in weighted spaces : preliminary version. 3rd cycle. Téhéran (Iran, 2006, pp.75.
More informationMAT 578 FUNCTIONAL ANALYSIS EXERCISES
MAT 578 FUNCTIONAL ANALYSIS EXERCISES JOHN QUIGG Exercise 1. Prove that if A is bounded in a topological vector space, then for every neighborhood V of 0 there exists c > 0 such that tv A for all t > c.
More informationExistence And Uniqueness Of Mild Solutions Of Second Order Volterra Integrodifferential Equations With Nonlocal Conditions
Applied Mathematics E-Notes, 9(29), 11-18 c ISSN 167-251 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Existence And Uniqueness Of Mild Solutions Of Second Order Volterra Integrodifferential
More informationExercise Solutions to Functional Analysis
Exercise Solutions to Functional Analysis Note: References refer to M. Schechter, Principles of Functional Analysis Exersize that. Let φ,..., φ n be an orthonormal set in a Hilbert space H. Show n f n
More informationVolterra Integral Equations of the First Kind with Jump Discontinuous Kernels
Volterra Integral Equations of the First Kind with Jump Discontinuous Kernels Denis Sidorov Energy Systems Institute, Russian Academy of Sciences e-mail: contact.dns@gmail.com INV Follow-up Meeting Isaac
More informationSpectral theory for compact operators on Banach spaces
68 Chapter 9 Spectral theory for compact operators on Banach spaces Recall that a subset S of a metric space X is precompact if its closure is compact, or equivalently every sequence contains a Cauchy
More informationMTH 503: Functional Analysis
MTH 53: Functional Analysis Semester 1, 215-216 Dr. Prahlad Vaidyanathan Contents I. Normed Linear Spaces 4 1. Review of Linear Algebra........................... 4 2. Definition and Examples...........................
More informationShort note on compact operators - Monday 24 th March, Sylvester Eriksson-Bique
Short note on compact operators - Monday 24 th March, 2014 Sylvester Eriksson-Bique 1 Introduction In this note I will give a short outline about the structure theory of compact operators. I restrict attention
More informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
More informationTopics in Harmonic Analysis Lecture 1: The Fourier transform
Topics in Harmonic Analysis Lecture 1: The Fourier transform Po-Lam Yung The Chinese University of Hong Kong Outline Fourier series on T: L 2 theory Convolutions The Dirichlet and Fejer kernels Pointwise
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 informationIdentification of a memory kernel in a nonlinear parabolic integro-differential problem
FACULTY OF ENGINEERING AND ARCHITECTURE Identification of a memory kernel in a nonlinear parabolic integro-differential problem K. Van Bockstal, M. Slodička and F. Gistelinck Ghent University Department
More information