for (int c0 = 0; c0 <= M; c0 += 1)
  S1(c0);