This paper presents some results of analytical work focused on finding the optimal description of the viscoelasticity of solid polymers. Memory functions of creep and relaxation are the basic functions in the constitutive equations. Applying a specific type of this kind of function is the purpose of this work. Based on the known exponential function, a new function called the “root function” has been created. Its very important feature is that the corresponding plots of spectra are of the same character as the curves of molecular weight distribution. The theory is verified by applying uniaxial creep tests. The experimental curves are compared with the numerical simulation based on the newly introduced root function. Apart from the purely phenomenological verification, the molecular weight distribution curves obtained by the traditional GPC (gel permeation chromatographic) method are compared with the spectra calculated analytically on the basis of the parameters calculated from the creep tests. A very good convergence of the curves confirms the selection of the memory function and the applied method.