1. A method to estimate at least one function of a probability distribution of at least one variable, where the variables consist of a primary function and arguments of said primary function, and where at least one said argument is a datum, and the data are separated into at least two groups and processed sequentially, comprising the steps of: (a) assigning prior bounds for said datum, said arguments, and said primary function; (b) assigning prior probability distributions within said prior bounds for all said arguments of said primary function; (c) generating realizations of said primary function from realizations of said arguments using Monte-Carlo technique; (d) determining at least one posterior bound among said realizations of said primary function and said arguments; (e) stopping the sequential processing when there are less than a minimum number of acceptable said realizations of said primary function within its said prior bounds; (f) stopping said sequential processing when at least one said posterior bound is inconsistent with at last one said prior bound; (g) determining at least one final bound after said data are sequentially processed; whereby at least one said function of said probability distribution of at least one said variable is estimated, where said variables consist of said values of said primary function and said arguments of said primary function, and where at least one said argument is said datum, and said data are separated into at least two said groups and processed sequentially.