`species=Bpen;
ii=1;
    for i=2:50
        for j=i:50
            y=find(distance(:,1)==species(i,1));
            x=find(distance(1,:)==species(1,j+1));
        Species_list(ii,1)=species(i,j+1);
        Distance(ii,1)=distance(y,x);
    ii=ii+1;
    end
end`