An \(hp\)-version spectral collocation method for nonlinear Volterra integro-differential equation with weakly singular kernels (Q1676917)

The numerical solutions of the following nonlinear Volterra integro-differential equations with weakly singular kernels are investigated \[ y'(t)=f(t,y(t))+\int\limits_{0}^{t}(t-s)^{-\mu}K(t,s,y(s))ds, \quad t\in(0,T], \quad y(0)=y_0, \tag{1} \] where \(0<\mu<1\), \(f(t,y)\) and \(K(t,s,y)\) are given functions, and \(y_0\) is the initial data. An efficient \(hp\)-version Legendre-Jacobi spectral collocation method is presented to solve this problem. Numerical results show that the \(p\)-version gives twice the rate of convergence as the \(h\)-version for singular solutions on quasi-uniform meshes. A priori error estimates under the \(H^1\)-norm are established for equation (1) with smooth solutions on arbitrary meshes and singular solutions on quasi-uniform meshes. Numerical results demonstrating the effectiveness of the proposed method are also suggested.