Numerical methods for hypersingular integrals on the real line

In the present paper the authors propose two numerical methods to approximate Hadamard transforms of the type H_p( f w_eta , t) =int_R rac{f(x)}{(x-t)^p}w_eta8x)dx where p is a nonnegative integer and w_eta (x) = e^{|x|^eta} ,eta > 1, is a Freud weight. One of the procedures employed here is based on a simple tool like the “truncated” Gaussian rule conveniently modified to remove numerical cancellation and overflow phenomena. The second approach is a process of simultaneous approximation of the functions This strategy can be useful in the numerical treatment of hypersingular integral equations. The methods are shown to be numerically stable and convergent and some error estimates in suitable Zygmund-type spaces are proved. Numerical tests confirming the theoretical estimates are given. Comparisons of our methods among them and with other ones available in literature are shown.