Solving a homogeneous Fredholm equation of the 2nd kind whose kernel has simple poles in the domain of interest Consider the Fredholm equation of the 2nd kind
$$ f(s) = \lambda \int_{-\infty}^{\infty} f(s') \Big(\sum_{n=1}^{N} g_n(s) h_n(s') \Big)  ds' , $$
with $f(s)$ an unknown function, $\lambda$ a constant, {$h_1(s') , . . . , h_N(s')$} and {$g_1(s) , . . . , g_N(s)$} two sets of known functions. 
Please note that this integral must be interpreted in the Principal Value sense. If the sets {$g_n(s)$}, {$h_n(s')$} contain only functions which decay sufficiently fast and contain no poles on the interval $(-\infty,\infty)$, then methods to solve this equation are covered in most introductory discussions of integral equations. However, what can be made of the case where at least one element of {$g_n(s)$} has at least one simple pole on the real axis? 
Questions:
(1) How can this equation be classified? Is it a singular integral equation even though it does not involve poles 'along the diagonal'?
(2) Is there an accepted reference for integral equations of this type?
Note: If it helps, I am interested, to start, in a specific example:
$$ f(s) = \frac{1}{a \beta} \int_{-\infty}^{\infty} f(s')  \frac{s \sinh(\frac{\pi s}{2})}{s \sinh(\frac{\pi s}{2}) - 1} \frac{1}{\cosh(\frac{\pi (s-s')}{2})}  ds' , $$
which has simple poles at $ s = \pm s_0 \approx \pm .7202 $. I know, from the solution of a transformed version of this equation that in $\lim_{a->\pm \infty} f(s)$ is proportional to a sum of Dirac deltas. So, clearly there is something pathological going on because the solution is not a function at all. It is not even clear to me that we may legitimately expand the kernel so that it follows the form that I have listed above.
I may need to solve several equations of this form, so I am interested in a general solution method.

Have you seen/considered homotopy perturbation techniques given here - http://www.sciencedirect.com/science/article/pii/S0960077905010362?
They may be applicable to your problem. 
