  We introduce an algorithm to compute the functions belonging to a suitable
set ${\mathscr F}$ defined as follows: $f\in {\mathscr F}$ means that $f(s,x)$,
$s\in A\subset {\mathbb R}$ being fixed and $x>0$, has a power series expansion
centred at $x_0=1$ with convergence radius greater or equal than $1$; moreover,
it satisfies a functional equation of step $1$ and the Euler-Maclaurin
summation formula can be applied to $f$. Denoting the Euler gamma-function as
$\Gamma$, we will show that, for $x>0$, $\log \Gamma(x)$, the digamma function
$\psi(x)$, the polygamma functions $\psi^{(w)}(x)$, $w\in {\mathbb N}$,
$w\ge1$, and, for $s>1$ being fixed, the Hurwitz $\zeta(s,x)$-function and its
first partial derivative $\frac{\partial\zeta}{\partial s}(s,x)$ are in
${\mathscr F}$. In all these cases the coefficients of the involved power
series will depend on the values of $\zeta(u)$, $u>1$, where $\zeta$ is the
Riemann zeta-function. As a by-product, we will also show how to compute the
Dirichlet $L$-functions $L(s,\chi)$ and $L^\prime(s,\chi)$, $s> 1$, $\chi$
being a primitive Dirichlet character, by inserting the reflection formulae of
$\zeta(s,x)$ and $\frac{\partial\zeta}{\partial s}(s,x)$ into the first step of
the Fast Fourier Transform algorithm. Moreover, we will obtain some new
formulae and algorithms for the Dirichlet $\beta$-function and for the Catalan
constant $G$. Finally, we will study the case of the Bateman $G$-function and
of the alternating Hurwitz zeta-function, also known as the $\eta$-function; we
will show that, even if they are not in ${\mathscr F}$, our approach can be
adapted to handle them too. In the last section we will also describe some
tests that show a performance gain with respect to a standard multiprecision
implementation of $\zeta(s,x)$ and $\frac{\partial\zeta}{\partial s}(s,x)$,
$s>1$, $x>0$.
