Various algorithms have been proposed to solve the interconversion equation of linear viscoelasticity when Prony series are used for the relaxation and creep moduli, G(t) and J(t). With respect to a Prony series for G(t), the key step in recovering the corresponding Prony series for J(t) is the determination of the coefficients {jk} of terms in J(t). Here, the need to solve a poorly conditioned matrix equation for the {jk} is circumvented by deriving elementary and easily evaluated analytic formulae for the {jk} in terms of the derivative dG(s)/ds of the Laplace transform G(s) of G(t).