%%%%%%%%%%%%%%%%%%%%% Parametrization of Geodesic Ray %%%%%%%%%%%%%%%%%%%%%
%
% GEORAY
% Usage: GEORAY(seg,d)
% Inputs:  seg = [x y] -  x,y are points inside or on the unit circle  
%
% Outputs: the hyperbolic geodesic arc from x to y is constructed for 
%          a hyperbolic distance of d 

function w = georay(seg,d)             

L = raymat(seg);
r = tanh(d);
z = 0:abs(r)/(narcpnts-1):abs(r);
if (r < 0) z = -z; end; 
w = lftact(L,z);               
clear z;