Solution of second kind Fredholm integral equations by means of Gauss and anti-Gauss quadrature rules

This paper is concerned with the numerical approximation of Fredholm integral equations of the second kind. A Nyström method based on the anti-Gauss quadrature formula is developed and investigated in terms of stability and convergence in appropriate weighted spaces. The Nyström interpolants corresponding to the Gauss and the anti-Gauss quadrature rules are proved to furnish upper and lower bounds for the solution of the equation, under suitable assumptions which are easily verified for a particular weight function. Hence, an error estimate is available, and the accuracy of the solution can be improved by approximating it by an averaged Nyström interpolant. The effectiveness of the proposed approach is illustrated through different numerical tests.

Mode Stresses for the Interaction between Two Inclined Cracks in Bonded Two Half Planes

The various mode of stresses for the interaction between two inclined cracks in the upper part of bonded two half planes which are normal stress (Mode I), shear stress (Mode II), tearing stress (Mode III) and mixed stress was studied. For this problem, the modified complex potentials (MCPs) method was used to develop the new system of hypersingular integral equations (HSIEs) by applying the conditions for continuity of resultant force and displacement functions with the unknown variable of crack opening displacement (COD) function and the right hand terms are the tractions along the crack. The curve length coordinate method and Gauss quadrature rules were used to solve numerically the obtained HSIEs to compute the stress intensity factors (SIFs) in order to determine the strength of the materials containing cracks. Numerical solutions presented the characteristic of nondimensional SIFs at the cracks tips. It is obtained that the various stresses and the elastic constants ratio are influences to the value of nondimensional SIFs at the crack tips.

Rational averaged gauss quadrature rules

It is important to be able to estimate the quadrature error in Gauss rules. Several approaches have been developed, including the evaluation of associated Gauss-Kronrod rules (if they exist), or the associated averaged Gauss and generalized averaged Gauss rules. Integrals with certain integrands can be approximated more accurately by rational Gauss rules than by Gauss rules. This paper introduces associated rational averaged Gauss rules and rational generalized averaged Gauss rules, which can be used to estimate the error in rational Gauss rules. Also rational Gauss-Kronrod rules are discussed. Computed examples illustrate the accuracy of the error estimates determined by these quadrature rules.

Stable Numerical Solutions Preserving Qualitative Properties of Nonlocal Biological Dynamic Problems

This paper deals with solving numerically partial integrodifferential equations appearing in biological dynamics models when nonlocal interaction phenomenon is considered. An explicit finite difference scheme is proposed to get a numerical solution preserving qualitative properties of the solution. Gauss quadrature rules are used for the computation of the integral part of the equation taking advantage of its accuracy and low computational cost. Numerical analysis including consistency, stability, and positivity is included as well as numerical examples illustrating the efficiency of the proposed method.

Error estimates for certain cubature formulae

We estimate the errors of selected cubature formulae constructed by the product of Gauss quadrature rules. The cases of multiple and (hyper-)surface integrals over n-dimensional cube, simplex, sphere and ball are considered. The error estimates are obtained as the absolute value of the difference between cubature formula constructed by the product of Gauss quadrature rules and cubature formula constructed by the product of corresponding Gauss-Kronrod or corresponding generalized averaged Gaussian quadrature rules. Generalized averaged Gaussian quadrature rule ?2l+1 is (2l + 1)-point quadrature formula. It has 2l + 1 nodes and the nodes of the corresponding Gauss rule Gl with l nodes form a subset, similar to the situation for the (2l + 1)-point Gauss-Kronrod rule H2l+1 associated with Gl. The advantages of bG2l+1 are that it exists also when H2l+1 does not, and that the numerical construction of ?2l+1, based on recently proposed effective numerical procedure, is simpler than the construction of H2l+1.