  A multidimensional generalization of Melvin's solution for an arbitrary
simple Lie algebra $\cal G$ is presented. The gravitational model contains n
2-forms and $l \geq n$ scalar fields, wheren is the rank of $\cal G$. The
solution is governed by a set of n functions obeying n ordinary differential
equations with certain boundary conditions. It was conjectured earlier that
these functions should be polynomials (the so-called fluxbrane polynomials). A
program (in Maple) for calculating of these polynomials for classical series of
Lie algebras is suggested (see Appendix). The polynomials corresponding to the
Lie algebra D_4 are obtained. It is conjectured that the polynomials for A_n-,
B_n- and C_n-series may be obtained from polynomials for D_{n+1}-series by
using certain reduction formulas.
