function Y = Ehren(X,d )

V=rand(1);
    Y=X +(V>(X/d))-(V<(X/d));


endfunction