%%%%%%%%%%%%%%%%%%%%% Reflection in a geodesic  %%%%%%%%%%%%%%%%%%%%%
%
% REFLECT
% Usage: REFLECT(seg,z)
% Inputs:  seg = [x y] -  x,y are points inside or on the unit circle  
%
% Outputs:  each entry of the matrix z is reflected in the  
%           geodesic connecting x and y and the result is returned 
%           in the matrix w.

function w = reflect(seg,z)             
    L = reflmat(seg);
    w = conj(lftact(L,z));