R. Hirota and K. Kimura discovered integrable discretizations of the Euler and the Lagrange tops, given by birational maps. Their method is a specialization to the integrable context of a general discretization scheme introduced by W. Kahan and appli
cable to any vector field with a quadratic dependence on phase variables. According to a proposal by T. Ratiu, discretizations of the Hirota-Kimura type can be considered for numerous integrable systems of classical mechanics. Due to a remarkable and not well understood mechanism, such discretizations seem to inherit the integrability for all algebraically completely integrable systems. We introduce an experimental method for a rigorous study of integrability of such discretizations. Application of this method to the Hirota-Kimura type discretization of the Clebsch system leads to the discovery of four functionally independent integrals of motion of this discrete time system, which turn out to be much more complicated than the integrals of the continuous time system. Further, we prove that every orbit of the discrete time Clebsch system lies in an intersection of four quadrics in the six-dimensional phase space. Analogous results hold for the Hirota-Kimura type discretizations for all commuting flows of the Clebsch system, as well as for the $so(4)$ Euler top.
