      function [reflect] = refl_xz
%     --------------------------------------------
%     reflection operator with respect to xz-plane    
%     --------------------------------------------
%     reflect = sub-set of finite motion matrix
      reflect = eye(3)   ;
      reflect(2,2) = -1 ;
%     end of refl_xz
