function z=mod1(x,n)

% z=mod1(x,n)
% z=x-n.*floor((x-.5)./n)
% with this definition, mod1(n,n)=n

z=x-n.*floor((x-.5)./n);
