Text this: Uncertain viscoelastic models with fractional order: a new spectral tau method to study the numerical simulations of the solution