The main results on the elaboration of the algorithms and programs for the search of the ‘hidden periodicities’ in the ‘unevenly spaced’ data, are discussed. They are mainly based on the application of the cubic spline-functions practically for all purposed of the variable star research — period determination, approximation of the phase curve, search for the possible period changes and/or the shape of the phase curve, detection of the periodic components in the case of the multi-harmonic oscillations, restoration of the smoothed function by removing the ‘apparat function’, numerical integration etc. Contrary to the commonly used splines with the arguments of the basic points coinciding with that of the real observations, we use the smoothing by the spline with the number of the basic points m, which is smaller than the number of the observations n. The main expressions for the corresponding fitting curves and their statistical properties were published by Andronov (1988, 1989), as well as the FORTRAN programs for the fits and the integration (Andronov, 1986, 1987, respectively).