A maximum-likelihood algorithm for determining precession corrections and stellar kinematic parameters is presented. It takes into account the ellipsoidal distribution of residual velocities, as well as observational errors. Numerical tests of the algorithm and the least squares method using a set of synthetic data are described. Both the algorithm of maximum-likelihood and the least squares method are applied to the 512 FK4/FK4 sup. stars selected by Fricke (1967).