In practice, the probability density function p(y) is often not available, or is very difficult
to estimate. Instead, we might have S different scenarios Y = (y[1], y[2], . . . , y[S]
) that are sampled from the probability distribution or that have been obtained from computer
simulations. Evaluating the auxiliary function F˜β(x, α) using the scenarios Y, we have