Proper elements for highly inclined asteroidal orbits

Based on Williams' work and rewritten in action angle variables, a method for the calculation of proper elements is here presented. The averaging over the long periodic terms is performed by the semi numerical method developed by Henrard (1990); no series expansion in eccentricity or inclination of the asteroid is used which allows calculating proper elements for highly inclined orbits. Conversely, the theory is truncated at the first degree in the eccentricity and the inclination of the perturbing planets. A few tests about accuracy and consistency are presented.