Závoti (2002) presented the mathematical description of the interpolation method especially for modeling the orbit of artificial satellites, which is suitable for approaching only certain 9 points. The task in this form stems from Grafarend and Schaffrin's (1993) study. During the time passed since the elaboration of the method, the generalization of the algorithm became necessary in the case when we have an arbitrary amount of measurement points, which must be approached according to a certain principle. The generalized method was successfully applied for modeling geodynamical processes, for describing the motion of the Earth's poles and for analyzing economical time series.