load hampshire_proximities.dat

inpermone = [2 6 12 19 22 1 11 4 5 23 21 18 3 27 20 24 17 10 25 9 26 13 8 16 7 15 14];
inpermtwo = [18 4 19 16 3 11 15 17 5 14 27 7 8 13 12 21 20 6 2 22 9 10 25 23 1 26 24];

[outpermone,outpermtwo,coordone,coordtwo,fitone,fittwo,addconone,addcontwo,vaf] = ...
    biscalqa_tied(hampshire_proximities,targlin(27),targlin(27),inpermone,inpermtwo,1:27,1:27,3,0)

tiedcoordone = [1 1 1 2 2 2 3 3 3 4 4 4 5 5 5 6 6 6 7 7 7 8 8 8 9 9 9];
tiedcoordtwo = [1 1 1 2 2 2 3 3 3 4 4 4 5 5 5 6 6 6 7 7 7 8 8 8 9 9 9];

[outpermone,outpermtwo,coordone,coordtwo,fitone,fittwo,addconone,addcontwo,vaf] = ...
    biscalqa_tied(hampshire_proximities,targlin(27),targlin(27),inpermone,inpermtwo,...
    tiedcoordone,tiedcoordtwo,3,0)


axes = zeros(27,2);
        
        for i = 1:27
            
            axes(outpermone(i),1) = coordone(i);
            axes(outpermtwo(i),2) = coordtwo(i);
        end
        
        plot(axes(1:27,1),axes(1:27,2),'ko')
        
        hold on
        
        for i = 1:27
            
            objectlabels{i,1} = int2str(i);
            
        end
        
        text(axes(1:27,1),axes(1:27,2),objectlabels,'fontsize',10,'verticalalignment','bottom')
        
        