docs: private static void normalizeLMR(double[] sum, int nmom) {
      for(int k = nmom - 1; k >= 1; --k) {
        double p = ((k & 1) == 0) ? +1 : -1;
        double temp = p * sum[0];
        for(int i = 0; i < k; i++) {
          double ai = i + 1.;
          p *= -(k + ai) * (k - i) / (ai * ai);
          temp += p * sum[i + 1];
        }
        sum[k] = temp;
      }
    }