  The large-n expansion is developed for the study of critical behaviour of
d-dimensional systems at m-axial Lifshitz points with an arbitrary number m of
modulation axes. The leading non-trivial contributions of O(1/n) are derived
for the two independent correlation exponents \eta_{L2} and \eta_{L4}, and the
related anisotropy index \theta. The series coefficients of these 1/n
corrections are given for general values of m and d with 0<m<d and
2+m/2<d<4+m/2 in the form of integrals. For special values of m and d such as
(m,d)=(1,4), they can be computed analytically, but in general their evaluation
requires numerical means. The 1/n corrections are shown to reduce in the
appropriate limits to those of known large-n expansions for the case of
d-dimensional isotropic Lifshitz points and critical points, respectively, and
to be in conformity with available dimensionality expansions about the upper
and lower critical dimensions. Numerical results for the 1/n coefficients of
\eta_{L2}, \eta_{L4} and \theta are presented for the physically interesting
case of a uniaxial Lifshitz point in three dimensions, as well as for some
other choices of m and d. A universal coefficient associated with the
energy-density pair correlation function is calculated to leading order in 1/n
for general values of m and d.
