  We study fast approximation of integrals with respect to stationary
probability measures associated to iterated functions systems on the unit
interval. We provide an algorithm for approximating the integrals under certain
conditions on the iterated function system and on the function that is being
integrated. We apply this technique to estimate Hausdorff moments, Wasserstein
distances and Lyapunov exponents of stationary probability measures.
