%%%%%%%%%%%%%%%%%%%%% Intersection of Hyperbolic Lines %%%%%%%%%%%%%%%%%%%%%
%
% INTLINES
% Usage: INTLINES(zL1,zL2)
% Inputs:  zL1,zL2 are the euclidean centers of two lines L1 and L2
%
% Output: the intersection point of the two lines

function w = intlines(zL1,zL2)             

I = sqrt(-1);
w = zL1-zL2;
w = I*w/abs(w);
t = roots([1,-2*real(zL1*w'),1]);
if abs(t(1)) < 1
   t=t(1)
else t=t(2);
end
   w=t*w;