Interpolation method for solving Volterra integral equations with weakly singular kernel using an advanced barycentric Lagrange formula

Abstract

We presented an interpolation method for solving weakly singular Volterra integral equations of the second kind (SVK2). The method based on the barycentric Lagrange interpolation.. For the chosen nodes of the two singular kernel variables, we created two rules that confirm that the denominator of the kernel will never become zero or have an imaginary value. The product of four matrices, one of which is the square coefficients matrix of the barycentric functions, expresses both the unknown and data functions. The singular kernel is interpolated twice with to its two variables. Furthermore, the kernel is represented by five matrices: two of which are the monomial basis of the kernel's variables, and one matrix expresses the kernel's values at the two sets of nodes for the kernel's variables. Without using the collocation approach, we get the matrix of the unknown coefficients by substituting the interpolated solution into the integral equation to obtain an equivalent algebraic system. By solving the algebraic system, we get the interpolated solution. We solved five examples, two for non-singular equations and three for weakly singular equations. By adopting lower interpolation degrees, the interpolated solutions for non-singular equations are shown to be equal to the exact ones, whereas, for singular equations, they are found to strongly converge to the exact ones and are superior to those obtained by other cited methods. This ensures the proposed method's uniqueness and high accuracy.