  Elliptic polylogarithms can be defined as iterated integrals on a genus-one
Riemann surface of a set of integration kernels whose generating series was
already considered by Kronecker in the 19th century. In this article, we employ
the Schottky parametrization of a Riemann surface to construct higher-genus
analogues of Kronecker's generating series, which we refer to as
Schottky-Kronecker forms. Our explicit construction generalizes ideas from
Bernard's higher-genus construction of the Knizhnik-Zamolodchikov connection.
Integration kernels generated from the Schottky-Kronecker forms are defined as
Poincar\'e series. Under technical assumptions, related to the convergence of
these Poincar\'e series on the underlying Riemann surface, we argue that these
integration kernels coincide with a set of differentials defined by Enriquez,
whose iterated integrals constitute higher-genus analogues of polylogarithms.
Enriquez' original definition is not well-suited for numerical evaluation of
higher-genus polylogarithms. In contrast, the Poincar\'e series defining our
integration kernels can be evaluated numerically for real hyperelliptic curves,
for which the above-mentioned convergence assumptions can be verified. We
numerically evaluate several examples of genus-two polylogarithms, thereby
paving the way for numerical evaluation of hyperelliptic analogues of
polylogarithms.
