  We present a publicly available software for exponential integrators that
computes the $\varphi_l(z)$ functions using polynomial interpolation. The
interpolation method at Leja points have recently been shown to be competitive
with the traditionally-used Krylov subspace method. The developed framework
facilitates easy adaptation into any Python software package for time
integration.
