function Rxz=fzcompose(Rxy,Ryz,mode)
% Usage  Rxz=fzcompose(Rxy,Ryz,mode)
% fuzzy relation composition
% mode = 0 (default): Max-min composition
% mode = 1: Max-product composition
% Rxy: L x M,  Ryz: M x N, Rxz: L x N
% (C) 2005 by Yu Hen Hu
% created: 12/4/2005

if nargin<3, mode=0; end % default mode = 0

[L,M]=size(Rxy);
[M2,N]=size(Ryz);
if M~=M2,
   error(' dimension of Rxy and Ryz does not match, abort!');
end
Rxz=zeros(L,N);

for i=1:L, % for each row of Rxy matrix
   for j=1:N, % for each column of Ryz matrix
      if mode==0, % max-min composition
         Rxz(i,j)=max(min([Rxy(i,:); Ryz(:,j)']));
      elseif mode==1,
         Rxz(i,j)=max(Rxy(i,:).*Ryz(:,j)');
      end
   end
end