  To calculate the intercept of the multigluon system in a symmetric spatial
configuration a variational method is developed based on a complete system of
one-gluon functions. The method is applied to two- and three- gluon cases to
compare with the known results. The convergence turns out rather slow. Ways to
improve results are discussed.
