An efficient method is proposed to determine the deformation function of a viscoelastic material from experimental data. The deformation function is assumed to be an integral operator with Rabotnov's fractional-exponential kernel or a sum of such kernels. This representation enables effective use of the method of operator continued fractions. To illustrate the method, deformation data for polymethylmethacrylate are used. The viscoelastic characteristics of a composite based on this material are obtained using the method of operator continued fractions Introduction. There are different methods to represent creep data using various functions of time. Such functions are needed to solve viscoelastic problems to obtain many characteristics of the material under investigation. Experimental creep data can be represented: (a) as a table if strains ε i = ε(t i ) at t i , i = 0, 1, …, n, are known (ε max is the maximum strain corresponding to the time t max ) and (b) graphically if strain ε = ε(t) at each t is known. This representation implies either a spline running through experimental points or already obtained standard creep curves.Many hereditary-elasticity problems are solved using the following analytic time-dependence of strain: