function [ alpha_r ] = RestrictAngle(alpha)
%Restricts any value of alpha to +/- pi radians

  alpha_r = mod(alpha,2*pi);
  idx = find(alpha_r > pi);
  alpha_r(idx) = alpha_r(idx)-2*pi;
end

