Methods are described for determining the a priori variations of elements of small-eccentricity orbits with a specified probability under the influence of instrument errors in the carrier vehicle control system and perturbing factors. A method is proposed for approximating the probability density functions of the orbit elements by a composition of a normal and a gamma-distribution function.