UMAT nonlinear viscoelastic behavior

I wonder how to write a UMAT of a nonlinear viscoelastic behavior of creep type (unidirectional schapery model) we have the young modulus, the Poisson's ratio, the interpolation functions of the non-linear parameters (g0,gk:two second-degree polynomials) and the couples of the prony series (Di,τi )