The numerical treatment of Love's integral equation having very small parameter

作者：Pastore Pietro^*

来源：Journal of Computational and Applied Mathematics, 2011, 236(6): 1267-1281.

DOI：10.1016/j.cam.2011.08.011

• Love's integral equation
• System of Fredholm integral equations
• Nystrom method
• Gaussian quadrature rule
• Condition number

In this paper, we describe a procedure to find the numerical solution of Love's integral equation
f (y) + 1/pi integral(1)(-1) c/(x-y)(2) + c(2) f(x)dx = 1. vertical bar Y vertical bar <= 1. 0 <= c is an element of R.
\A crucial role is played by the parameter c. In fact, the complex poles of the kernel function get closer to the real axis as c decreases, and numerical difficulties appear ill approximating the solution f. We reduce the above integral equation to an equivalent system of Fredholm integral equations that we solve using a stable and convergent method. Moreover, since the matrices of coefficients of the derived linear systems are structured and can become very "large", we focus our attention also on the computational aspects related to the numerical solution of such linear systems. Finally, we give some numerical tests considering different values of the parameter c, and compare our results with the ones obtained with other methods in the literature.

• 出版日期2011-10-15