where p(δ) represents the probability distribution of the decision variable disturbance.
In general, the distribution takes a normal distribution. Usually, this effective objective
function is not available. Hence, it is approximated, for instance, by a Monte Carlo
integration:
f ′robust(x) =
_N
i=1 f (x + δi)