The existence of a fast algorithm with an
O
A
(
n
(
log
n
)
2
)
{O_A}(n{(\log n)^2})
time complexity for the multiplication of generalized Hilbert matrices with vectors is shown. These matrices are defined by
(
B
p
)
i
,
j
=
1
/
(
t
i
−
s
j
)
p
{({B_p})_{i,j}} = 1/{({t_i} - {s_j})^p}
,
i
,
j
=
1
,
…
,
n
i,j = 1, \ldots ,n
,
p
=
1
,
…
,
q
p = 1, \ldots ,q
,
q
≪
n
q \ll n
, where
t
i
{t_i}
and
s
i
{s_i}
are distinct points in the complex plane and
t
i
≠
s
j
{t_i} \ne {s_j}
,
i
,
j
=
1
,
…
,
n
i,j = 1, \ldots ,n
. The major contribution of the paper is the stable implementation of the algorithm for two important sets of points, the Chebyshev points and the nth roots of unity. Such points have applications in the numerical approximation of Cauchy singular integral equations. The time complexity of the algorithm, for these special sets of points, reduces to
O
A
(
n
log
n
)
{O_A}(n\log n)
.