  We propose numerical integration methods for Choquet integrals where the
capacities are given by distortion functions of an underlying probability
measure. It relies on the explicit representation of the integrals for step
functions and can be seen as quasi-Monte Carlo methods in this framework. We
give bounds on the approximation errors in terms of the modulus of continuity
of the integrand and the star discrepancy.
