Results1.m

function [dis_mat] = Results1( images )
    n = length(images);
    dis_mat = zeros(n, n); % dissimilarity matrix
    torf = zeros(n, n); % 1 means have calculated dissimilarity of the comb
    hists = cell(1,n); % stores histogram of imgs
    %havecaled = zeros(1, n); % have calculated the histogram of this img
    
%%
	for row = 1:n
       dis_mat(row,row) = 0;
       torf(row,row) = 1;
    end
%%
    if (mod(n,2)==0)
        for row = 1:2:n
           col = row+1;
           img1 = images{row};
           img2 = images{col};
           [dis_mat(row,col) hist1 hist2] = ColorDist1(img1, img2);
           %havecaled(row) = 1;  havecaled(col) = 1;
           hists{row} = hist1;  hists{col} = hist2;
           dis_mat(col,row) = dis_mat(row,col);
           torf(row,col) = 1; torf(col,row) = 1;
           %fprintf('processing: %d, %d\n',row,col);
        end
    else
        for row = 1:2:n-1
           col = row+1;
           img1 = images{row};
           img2 = images{col};
           [dis_mat(row,col) hist1 hist2] = ColorDist1(img1, img2);
           %havecaled(row) = 1;  havecaled(col) = 1;
           hists{row} = hist1;  hists{col} = hist2;
           dis_mat(col,row) = dis_mat(row,col);
           torf(row,col) = 1; torf(col,row) = 1;
        end
        img1 = images{1};
        img2 = images{n};
        [dis_mat(row,col) hist1 hist2] = ColorDist1(img1, img2);
        hists{1} = hist1;  hists{n} = hist2;
        dis_mat(col,row) = dis_mat(row,col);
        torf(row,col) = 1; torf(col,row) = 1;
    end
 %%   
    for row = 1:n
       for col = 1:n
               if (torf(col,row)==1)
                   dis_mat(row,col) = dis_mat(col,row);
                   torf(row,col) = 1;
               else
                   hist1 = hists{row};
                   rhist1 = hist1(1,:); 
                   ghist1 = hist1(2,:); 
                   bhist1 = hist1(3,:);
                   hist2 = hists{col};
                   rhist2 = hist2(1,:); 
                   ghist2 = hist2(2,:); 
                   bhist2 = hist2(3,:);
                    a=0;
                    for k=0:255
                       a = a + (rhist1(k+1) - rhist2(k+1))^2;
                    end
                    b=0;
                    for k=0:255
                       b = b + (ghist1(k+1) - ghist2(k+1))^2;
                    end
                    c=0;
                    for k=0:255
                       c = c + (bhist1(k+1) - bhist2(k+1))^2;
                    end
                   dis_mat(row,col) = sqrt(a)+sqrt(b)+sqrt(c); dis_mat(col,row) = dis_mat(row,col);
                   torf(row,col) = 1; torf(col,row) = 1;
               end
           %fprintf('row: %d, col: %d\n',row,col);
       end
    end
end