Convergence of the piecewise linear approximation and collocation method for a hypersingular integral equation on a closed surface

We study the numerical solution of a linear hypersingular integral equation arising when solving the Neumann boundary value problem for the Laplace equation by the boundary integral equation method with the solution represented in the form of a double layer potential. The integral in this equation is understood in the sense of Hadamard finite value. We construct quadrature formulas for the integral occurring in this equation based on a triangulation of the surface and an application of the linear approximation to the unknown function on each of the triangles approximating the surface. We prove the uniform convergence of the quadrature formulas at the interpolation nodes as the triangulation size tends to zero. A numerical solution scheme for this integral equation based on the suggested quadrature formulas and the collocation method is constructed. Under additional assumptions about the shape of the surface, we prove a uniform estimate for the error in the numerical solution at the interpolation nodes.     

A remark on the application of interpolatory quadrature rules to the numerical solution of singular integral equations

A Cauchy type singular integral equation of the first or the second kind can be numerically solved either directly or after its reduction (by the usual regularization procedure) to an equivalent Fredholm integral equation of the second kind. The equivalence of these two methods (that is, the equivalence both of the systems of linear algebraic equations to which the singular integral equation is reduced and of the natural interpolation formulae) is proved in this paper for a class of Cauchy type singular integral equations of the first kind and of the second kind (but with constant coefficients) for general interpolatory quadrature rules under sufficiently mild assumptions. The present results constitute an extension of a series of previous results concerning only Gaussian quadrature rules, based on the corresponding orthogonal polynomials and their properties.     

Approximate solution to integral equation with logarithmic kernel of special form

Based on the quadrature formula with non-negative coefficients for integral with a special logarithmic kernel, we construct and substantiate a computational pattern for solving integral equation derived from the boundary-value problem for a function, which is harmonic in the unit disk under the boundary condition of the third kind. We obtain uniform estimates of deviations of the quadrature formula and the approximate solution to integral equation.     

Numerical Quadrature of Periodic Singular Integral Equations

This paper presents quadrature formulae for the numerical integrationof a singular integral equation with Hilbert kernel. The formulaeare based on trigonometric interpolation. By integration a quadratureformula for an integral with a logarithmic singularity is obtained.Finally it is demonstrated how a singular integral equationwith infinite support can be solved by use of the precedingformulae.     

An efficient pseudo‐spectral Legendre–Galerkin method for solving a nonlinear partial integro‐differential equation arising in population dynamics

The pseudo‐spectral Legendre–Galerkin method (PS‐LGM) is applied to solve a nonlinear partial integro‐differential equation arising in population dynamics. This equation is a competition model in which similar individuals are competing for the same resources. It is a kind of reaction–diffusion equation with integral term corresponding to nonlocal consumption of resources. The proposed method is based on the Legendre–Galerkin formulation for the linear terms and interpolation operator at the Chebyshev–Gauss–Lobatto (CGL) points for the nonlinear terms. Also, the integral term, which is a kind of convolution, is directly computed by a fast and accurate method based on CGL interpolation operator, and thus, the use of any quadrature formula in its computation is avoided. The main difference of the PS‐LGM presented in the current paper with the classic LGM is in treating the nonlinear terms and imposing boundary conditions. Indeed, in the PS‐LGM, the nonlinear terms are efficiently handled using the CGL points, and also the boundary conditions are imposed strongly as collocation methods. Combination of the PS‐LGM with a semi‐implicit time integration method such as second‐order backward differentiation formula and Adams‐Bashforth method leads to reducing the complexity of computations and obtaining a linear algebraic system of equations with banded coefficient matrix. The desired equation is considered on one and two‐dimensional spatial domains. Efficiency, accuracy, and convergence of the proposed method are demonstrated numerically in both cases. Copyright © 2012 John Wiley & Sons, Ltd.     

一类含中介值定积分等式证明题的构造

利用多项式插值理论,结合数值求积公式代数精度的概念,给出了一种构造和证明一类含中介值定积分等式证明题的新方法,并给出多个应用实例.     

A High-Order Piecewise Quartic Spline Rule for Hadamard Integral and Its Application of the Cavity Scattering

We develop a fourth-order piecewise quartic spline rule for Hadamard integral. The quadrature formula of Hadamard integral is obtained by replacing the integrand function with the piecewise quartic spline interpolation function. We establish corresponding error estimates and analyze the numerical stability. The rule can achieve fourth-order convergence at any point in the interval, even when the singular point coincides with the grid point. Since the derivative information of the integrand is not required, the rule can be easily applied to solve many practical problems. Finally, the quadrature formula is applied to solve the electromagnetic scattering from cavities with different wave numbers, which improves the whole accuracy of the solution. Numerical experiments are presented to show the efficiency and accuracy of the theoretical analysis.     

偏微分方程边值反问题的数值方法研究

本文研究了奇异积分方程在反边值问题中的应用问题.利用圆周上的自然积分方程及其反演公式,把Laplace方程的边值反问题转化为一对超奇异积分方程和弱奇异积分方程的组合,通过选取三角插值近似奇异积分的计算并构造相应的配置格式,并使用Tikhonov正则化方法求解所得到的线性方程组.数值实验表明了该方法的有效性.     

An exponentially convergent algorithm for nonlinear differential equations in Banach spaces

An exponentially convergent approximation to the solution of a nonlinear first order differential equation with an operator coefficient in Banach space is proposed. The algorithm is based on an equivalent Volterra integral equation including the operator exponential generated by the operator coefficient. The operator exponential is represented by a Dunford-Cauchy integral along a hyperbola enveloping the spectrum of the operator coefficient, and then the integrals involved are approximated using the Chebyshev interpolation and an appropriate Sinc quadrature. Numerical examples are given which confirm theoretical results.

     

Rational interpolation and approximate solution of integral equations

In the space of continuous periodic functions, we construct interpolation rational operators, use them to obtain quadrature formulas with positive coefficients which are exact on rational trigonometric functions of order 2n, and suggest an algorithm for an approximate solution of integral equations of the second kind. We estimate the accuracy of the approximate solution via the best trigonometric rational approximations to the kernel and the right-hand side of the integral equation.     

AN EXTREMAL APPROACH TO BIRKHOFF QUADRATURE FORMULAS

1. Introduction and Main ResultsIn tfor paPer we shaJl use the ddstions and notations of [3l. Let E = (e'k)7t' kt. be anincidence matrir with entries consisting of zeros and ones and satisfying lEl:= Z.,* ei* = n + 1(here we allow a zero row ). Furthermore, in wha follOws we assume that(A) E satisfies the P6lya condition(B) all sequences of E in the interior rows, 0 < i < m + 1, are even.Let Sm denote the set of poiats X = (xo, z1 l "') xm, x.+1) fOr whichand Sm its clOusure. If some O…     

A random integral quadrature method for numerical analysis of the second kind of Volterra integral equations

In this paper, a novel meshless technique termed the random integral quadrature (RIQ) method is developed for the numerical solution of the second kind of the Volterra integral equations. The RIQ method is based on the generalized integral quadrature (GIQ) technique, and associated with the Kriging interpolation function, such that it is regarded as an extension of the GIQ technique. In the GIQ method, the regular computational domain is required, in which the field nodes are scattered along straight lines. In the RIQ method however, the field nodes can be distributed either uniformly or randomly. This is achieved by discretizing the governing integral equation with the GIQ method over a set of virtual nodes that lies along straight lines, and then interpolating the function values at the virtual nodes over all the field nodes which are scattered either randomly or uniformly. In such a way, the governing integral equation is converted approximately into a system of linear algebraic equations, which can be easily solved.     

The best quadrature formula based on Hermite information for the class <Emphasis Type="Italic">KW</Emphasis><Superscript>2</Superscript> [<Emphasis Type="Italic">a,b</Emphasis>]

The best quadrature formula has been found in the following sense: for a function whose norm of the second derivative is bounded by a given constant and the best quadrature formula for the approximate evaluation of integration of that function can minimize the worst possible error if the values of the function and its derivative at certain nodes are known. The best interpolation formula used to get the quadrature formula above is also found. Moreover, we compare the best quadrature formula with the open compound corrected trapezoidal formula by theoretical analysis and stochastic experiments.     

The best quadrature formula based on Hermite information for the class KW~2[a,b]

The best quadrature formula has been found in the following sense:for afunction whose norm of the second derivative is bounded by a given constant and thebest quadrature formula for the approximate evaluation of integration of that function canminimize the worst possible error if the values of the function and its derivative at certainnodes are known.The best interpolation formula used to get the quadrature formula aboveis also found.Moreover,we compare the best quadrature formula with the open compoundcorrected trapezoidal formula by theoretical analysis and stochastic experiments.     

A generalization of discrete Gronwall inequality and its application to weakly singular Volterra integral equation of the second kind

This paper presents a new discrete Gronwall inequality. Using the inequality, we prove convergence and error estimate of the numerical solutions of the second weakly singular Volterra integral equation, where discrete equation is derived by Novot's quadrature formula.     

Solving multivariable mathematical models by the quadrature and cubature methods

The utilization and generalization of quadrature and cubature approximations for numerical solution of mathematical models of multivariable transport processes involving integral, differential, and integro-differential operators, and for numerical interpolation and extrapolation, are presented. The methodology for determination of the quadrature and cubature weights for composite operators is developed to accommodate for general functional representations. Application of these methods is demonstrated by solving two-dimensional steady-state and one-dimensional transient-state problems. The solutions are compared with exact-analytical solutions to evaluate the performance of these methods. It is demonstrated that the quadrature and cubature approximations are simple and universal; i.e., the same formula is applicable irrespective of the order of accuracy of the numerical approximation, the type of linear operator, and the number of temporal and/or spatial variables. Since the quadrature and cubature methods can produce solutions with sufficient accuracy even when using fewer discrete points, both the programming task and computational effort are reduced considerably. Therefore, the quadrature and cubature methods appear to be very practical in solving the mathematical models of a variety of transport processes. © 1994 John Wiley & Sons, Inc.     

Numerical solution of modified Black–Scholes equation pricing stock options with discrete dividend

This paper deals with the numerical solution of the modified Black–Scholes equation modelling the valuation of stock options with discrete dividend payments. By using a delta-defining sequence of the involved generalized Dirac delta function and applying the Mellin transform, an integral formula for the solution is obtained. Then, numerical quadrature approximations and illustrative examples are given.     

Error bounds and estimates for eigenvalues of integral equations

Summary Approximate solutions of the linear integral equation eigenvalue problem can be obtained by the replacement of the integral by a numerical quadrature formula and then collocation to obtain a linear algebraic eigenvalue problem. This method is often called the Nyström method and its convergence was discussed in [7]. In this paper computable error bounds and dominant error terms are derived for the approximation of simple eigenvalues of nonsymmetric kernels.     

Uniform convergence of the rectangle method for singular integral equation with the Hölder density

We study an approximate method for solving singular integral equations. It implies an approximation of a singular operator by means of a compound quadrature formula similar to the rectangle one. The corresponding systems of linear algebraic equations are solvable if so is the integral equation, while its coefficients satisfy the strong ellipticity condition. Under these restrictions we obtain a bound for the rate of convergence of solutions of systems of linear equations to the solution of the considered integral equation in the uniform vector norm.     

On the convergence of the Nyström method for the integral equation eigenvalue problem

Summary Approximate solutions of the linear integral equation eigenvalue problem can be obtained by the replacement of the integral by a numerical quadrature formula and then collocation to obtain a linear algebraic eigenvalue problem. This method is often called the Nyström method and a framework for its error analysis was introduced by Noble [15]. In this paper the convergence of the method is considered when the integral operator is a compact operator from a Banach spaceX intoX.