% This is just the card shuffling map with small delta, on the domain 0≤x≤1. So it is just
% a small perturbation to the Dyadic Map. 
%
% My goal is to produce random deviates 0<x<1. Empirically, I find that delta less than 
% a few thousand epsilon leads to obvious, persistent artifacts in the histogram.

function x1 = delta_eps(x0)

persistent delta = eps*3000; 

x1 = mod(delta + 2*(1-delta)*x0, 1);
