The recent formula for decomposing the expectation value <S(2)> of the total spin operator for general (correlated) wave functions has been rewritten in terms of the cumulant and is realized numerically for the first time. The results confirm its conformity with the physical expectations.