  An algorithm for integration of polynomial functions with variable weight is
considered. It provides extension of the Gaussian integration, with appropriate
scaling of the abscissas and weights. Method is a good alternative to usually
adopted interval splitting.
