function y=modpi(x)

% y = modpi(x)

y = rem( x + pi, 2*pi ) - pi;

ii = find( y < -pi );
if ~isempty( ii )
    y(ii) = y(ii) + 2*pi;
end
