  Two types of finite series of products of harmonic numbers involving
nonnegative integer powers are evaluated, also yielding two other important
harmonic number identities. The recursion formulas for these sums are derived,
which are easily translated into a computer program.
