A Hybrid Collocation Method for Volterra Integral Equations with Weakly Singular Kernels

The commonly used graded piecewise polynomial collocation method for weakly singular Volterra integral equations may cause serious round-off error problems due to its use of extremely nonuniform partitions and the sensitivity of such time-dependent equations to round-off errors. The singularity preserving ( nonpolynomial) collocation method is known to have only local convergence. To overcome the shortcoming of these well-known methods, we introduce a hybrid collocation method for solving Volterra integral equations of the second kind with weakly singular kernels. In this hybrid method we combine a singularity preserving ( nonpolynomial) collocation method used near the singular point of the derivative of the solution and a graded piecewise polynomial collocation method used for the rest of the domain. We prove the optimal order of global convergence for this method. The convergence analysis of this method is based on a singularity expansion of the exact solution of the equations. We prove that the solutions of such equations can be decomposed into two parts, with one part being a linear combination of some known singular functions which reflect the singularity of the solutions and the other part being a smooth function. A numerical example is presented to demonstrate the effectiveness of the proposed method and to compare it to the graded collocation method

Spectral collocation method for compact integral operators

We propose and analyze a spectral collocation method for integral equations with compact kernels, e.g. piecewise smooth kernels and weakly singular kernels of the form $\frac{1}{|t-s|^\mu}, \; 0\u3c\mu\u3c1.$ We prove that 1) for integral equations, the convergence rate depends on the smoothness of true solutions $y(t)$. If $y(t)$ satisfies condition (R): $\|y^{(k)}\|_{L^\infty[0,T]}\leq ck!R^{-k}$}, we obtain a geometric rate of convergence; if $y(t)$ satisfies condition (M): $\|y^{(k)}\|_{L^{\infty}[0,T]}\leq cM^k$, we obtain supergeometric rate of convergence for both Volterra equations and Fredholm equations and related integro differential equations; 2) for eigenvalue problems, the convergence rate depends on the smoothness of eigenfunctions. The same convergence rate for the largest modulus eigenvalue approximation can be obtained. Moreover, the convergence rate doubles for positive compact operators. Our numerical experiments confirm our theoretical results

A Product Integration type Method for solving Nonlinear Integral Equations in L

This paper deals with nonlinear Fredholm integral equations of the second kind. We study the case of a weakly singular kernel and we set the problem in the space L 1 ([a, b], C). As numerical method, we extend the product integration scheme from C 0 ([a, b], C) to L 1 ([a, b], C)

A nonpolynomial collocation method for fractional terminal value problems

NOTICE: this is the author’s version of a work that was accepted for publication in Journal of Computational and Applied Mathematics. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Journal of Computational and Applied Mathematics, 275, February 2015, doi: 10.1016/j.cam.2014.06.013In this paper we propose a non-polynomial collocation method for solving a class of terminal (or boundary) value problems for differential equations of fractional order α, 0 < α < 1. The approach used is based on the equivalence between a problem of this type and a Fredholm integral equation of a particular form. Taking into account the asymptotic behaviour of the solution of this problem, we propose a non-polynomial collocation method on a uniform mesh. We study the order of convergence of the proposed algorithm and a result on optimal order of convergence is obtained. In order to illustrate the theoretical results and the performance of the method we present several numerical examples.The work was supported by an International Research Excellence Award funded through the Santander Universities scheme

High-Order Multivariate Spectral Algorithms for High-Dimensional Nonlinear Weakly Singular Integral Equations with Delay

One of the open problems in the numerical analysis of solutions to high-dimensional nonlinear integral equations with memory kernel and proportional delay is how to preserve the high-order accuracy for nonsmooth solutions. It is well-known that the solutions to these equations display a typical weak singularity at the initial time, which causes challenges in developing high-order and efficient numerical algorithms. The key idea of the proposed approach is to adopt a smoothing transformation for the multivariate spectral collocation method to circumvent the curse of singularity at the beginning of time. Therefore, the singularity of the approximate solution can be tailored to that of the exact one, resulting in high-order spectral collocation algorithms. Moreover, we provide a framework for studying the rate of convergence of the proposed algorithm. Finally, we give a numerical test example to show that the approach can preserve the nonsmooth solution to the underlying problems. © 2022 by the authors.King Saud University, KSUM. A. Zaky and A. Aldraiweesh extend their appreciation to Distinguished Scientist Fellowship Program (DSFP) at King Saud University (Saudi Arabia)