On determining the nonnegative relaxation spectrum of viscoelastic materials using complementary error functions