A M\"untz-Collocation spectral method for weakly singular volterra integral equations

In this paper we propose and analyze a fractional Jacobi-collocation spectral method for the second kind Volterra integral equations (VIEs) with weakly singular kernel $(x-s)^{-\mu},0<\mu<1$. First we develop a family of fractional Jacobi polynomials, along with basic approximation results for some weighted projection and interpolation operators defined in suitable weighted Sobolev spaces. Then we construct an efficient fractional Jacobi-collocation spectral method for the VIEs using the zeros of the new developed fractional Jacobi polynomial. A detailed convergence analysis is carried out to derive error estimates of the numerical solution in both $L^{\infty}$- and weighted $L^{2}$-norms. The main novelty of the paper is that the proposed method is highly efficient for typical solutions that VIEs usually possess. Precisely, it is proved that the exponential convergence rate can be achieved for solutions which are smooth after the variable change $x\rightarrow x^{1/\lambda}$ for a suitable real number $\lambda$. Finally a series of numerical examples are presented to demonstrate the efficiency of the method

Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel

The Jacobi pseudo-spectral Galerkin method for the Volterra integro-differential equations of the second kind with a weakly singular kernel is proposed in this paper. We provide a rigorous error analysis for the proposed method, which indicates that the numerical errors (in the Lωα,β2-norm and the L∞-norm) will decay exponentially provided that the source function is sufficiently smooth. Numerical examples are given to illustrate the theoretical results

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

An efficient spectral method for solving third-kind Volterra integral equations with non-smooth solutions

This paper is concerned with the numerical solution of the third kind Volterra integral equations with non-smooth solutions based on the recursive approach of the spectral Tau method. To this end, a new set of the fractional version of canonical basis polynomials (called FC-polynomials) is introduced. The approximate polynomial solution (called Tau-solution) is expressed in terms of FC-polynomials. The fractional structure of Tau-solution allows recovering the standard degree of accuracy of spectral methods even in the case of non-smooth solutions. The convergence analysis of the method is studied. The obtained numerical results show the accuracy and efficiency of the method compared to other existing methods