Method for the computation of the elements of the sum changed: it now runs faster and, in the "while" loop, the table M for the combinations of the indices is corrected (there was an error for dim > 3)

The lnapoch, lnbpoch and lngpoch elements are now calculated once for each value in M or Msum. This improves the execution time.

Typing error corrected in ref.