In this paper, we present an implementation of the harmonic polylogarithm of Remiddi and Vermaseren for Mathematica. It contains an implementation of the product algebra, the derivative properties, series expansion and numerical evaluation. The analytic continuation has been treated carefully, allowing the user to keep the control over the definition of the sign of the imaginary parts. Many options enables the user to adapt the behavior of the package to his specific problem.