function Z=LinearM3(X,Y)
%%GZCM3 function
global M

NC1=size(M,1);
NC2=size(M{1},1);



I=1;
for i=1:NC1
    
    AllZero=0;% 0 means all M{i}{j}(5) is zero for all j in i
    
    for j=1:NC2
       if M{i}{j}(5)==1 
           AllZero=1;
           break;
       end
    end
    
    if j==NC2 & AllZero==0
        continue;
    end
    
    TMIN=Linear(X,Y,i,j);  %TEMP MIN
      
    for k=j:NC2
        if M{i}{k}(5)==1 
            TT=Linear(X,Y,i,k);
            TMIN=min(TMIN,TT);
        end
    end
    
    
    MIN(:,:,I)=TMIN;
    I=I+1;
end


MAX=MIN(:,:,1);
if (I-1)>=2
  for i=2:(I-1)
    TMAX=MIN(:,:,i);    %TEMP MAX
    MAX=max(MAX,TMAX);
  end
end


Z=MAX;