A special system of canonical variables is considered. An algorithm for expanding the principal functions of Keplerian motion in new elements is presented. The advantage of the proposed system is a relatively small number of terms in the classical expansions of the unperturbed two-body problem. A method for expanding the time derivatives of the rectangular coordinates is proposed. Some estimates of the number of terms in the presented expansions have been obtained through numerical experiments.