We develop an 1/n expansion method for the systematic approxi- mation of the solutions of classical spin systems (tri-dimensional vec- tors) which can describe general spin-spin interactions and phases with complex magnetization pattern. This method is applied to study the phase transitions observed in pyrochlore crystals (e.g. Gd₂ Ti₂ O₇ ) and the results have qualitative agreement with the experimental ones.