An automatic code is proposed for the perturbation analysis of strongly non-linear two-degree-of-freedom systems with cubic non-linearities. The recently proposed [1] modified Lindstedt-Poincare method is adopted, both because of its excellent performance and of its straightforward implementation, The symbolic software Mathematica is used in order to speed up all the cumbersome algebra which is inherent to every perturbation method. (C) 1998 Academic Press Limited.