function y = shift(x)
%
% shift(x) = 2x mod 1
%
y = rem(2*x,1);