Tools for estimate (joint) cumulants and products of (joint) cumulants of a given population distribution using (multivariate) k-statistics and (multivariate) polykays, unbiased estimators with minimum variance. Tools for generating univariate and multivariate Faa di Bruno's formula and related polynomials, such as Bell polynomials, generalized complete Bell polynomials, partition polynomials, generalized partition polynomials. For more details see Di Nardo E., Guarino G., Senato D. (2009) , .