%%%%%%%%%%%%%%%%%%%%% Linear Fractional Transformation %%%%%%%%%%%%%%%%%%%
%
% LFTACT
% Usage: w = LFTACT(L,z)
% Inputs: L,z  - L in SL_2(C), and z a matrix of complex numbers       
%
% Output: The matrix obtained by finding the LFT-image of z under L, 
%         computed componentwise

function w = lftact(L,z)
  [m,n] = size(z);
  a = L(1,1);
  b = L(1,2);
  c = L(2,1);
  d = L(2,2);
  w = (a*z+b*ones(m,n))./(c*z+d*ones(m,n));