Fitting methods for relaxation modulus of viscoelastic materials