The time moment of appearance of zero relative velocity during displacements of contacting bodies in dynamical systems with dry friction can be precisely computed. In this work, there is proposed an event driven numerical scheme allowing for integration of discontinuous differential equations. The study shows, that a direct application of standard integration methods may produce some inaccurate estimations of stick-slip transitions observed during frictional relative displacements. Results of numerical computations obtained by means of the proposed method of approximation of one degree-of-freedom dynamical system with dry friction confirm good accuracy of the determination of sticking phases. By comparison of two apparently similar solutions there has been shown and explained a significant difference in their smoothness and accuracy of estimation of the sticking phase's boundaries. The P. Olejnik, J. Awrejcewicz and M. Fe膷kan phase observed in an approximate solution of the analyzed discontinuous counterpart is smooth, and its ends are more precisely estimated. Reasonable accuracy of the solution could be obtained by a decrease of fitting parameter of the proposed approximate method to an exact sticking phase. For a comparison, three numerical procedures performed in Python and Matlab have been included. The usefulness of this method was verified by calculation of elementary work against the friction force to yield an acceptable accuracy.Mathematics Subject Classification: 65D30, 68W25, 34A36, 74M10