Improving the accuracy of points coordinates determination from geodetic networks with the use of modern space geodesy techniques (GPS, SLR, VLBI, DORIS) identified factors that were earlier considered as insignificant. Such factors include the movement of tectonic plates. This made possible the transition from the static coordinate system model to the kinematic one. For the effect of the tectonic plates movement on the points coordinates you must know kinematic parameters of the tectonic plates movement – the Euler vector parameters – angular velocity of the tectonic plates rotation and the poles rotation coordinates (latitude and longitude) of tectonic plates. Many modern kinematic models of tectonic plates are used observations, performed by geophysical methods, to estimate kinematic parameters of the tectonic plates movement. The article presents the algorithm that allows to estimate the parameters of the tectonic plates movement from mathematical processing of geodetic measurements carried out at points distributed over the surface of tectonic plates. Based on the presented algorithm the EULER program was developed in the algorithmic FORTRAN language, which was experimentally tested. The analysis results showed that the parameters of the Euler vector obtained by the EULER program are in good agreement with the results obtained by using ITRF 2005 and ITRF 2008 models.