the difference between the 0.975 and 0.025 percentiles (divided by 2*1.96) as the accuracy estimate of a given parameter in the bootstrap method. The variability of the semi-regular pulsating star Z UMa is analyzed. The presence of multicomponent variability of an object, including, four periodic oscillations (188.88(3), 197.89(4) days and halves of both) and significant variability of the amplitudes and phases of individual oscillations are shown. Approximation using the parabolic spline is only slightly better than the asymptotic parabola, for our sampling of the complete interval. It is expectedly better for larger subintervals. The use of different complementary methods allows to get a statistically optimal phenomenological approximation.