function [pop,Calpop,t1,t2,t3,t4,t5,t6]=dis(pop)
M1=1;M2=1;M3=3;M4=1;M5=1;M6=1;
t=M1+M2+M3+M4+M5+M6+1;
[s,t]=size(pop);
for i=1:1:s
t1=zeros(1,M1);t2=zeros(1,M2);t3=zeros(1,M3);t4=zeros(1,M4);t5=zeros(1,M5);t6=zeros(1,M6);
j1=0;j2=0;j3=0;j4=0;j5=0;j6=0;d11=0;d12=0;d22=0;d24=0;d44=0;d26=0;d46=0;d55=0;d56=0;d66=0;
for j=1:1:t-1;
if find(pop(i,j)<=M1);
j1=j1+1;
t1(1,j1)=j;
Calpop(i,j)=1;
elseif find(pop(i,j)<=M1+M2);
j2=j2+1;
t2(1,j2)=j;
Calpop(i,j)=2;
elseif find(pop(i,j)<=M1+M2+M3);
j3=j3+1;
t3(1,j3)=j;
Calpop(i,j)=3;
elseif find(pop(i,j)<=M1+M2+M3+M4);
j4=j4+1;
t4(1,j4)=j;
Calpop(i,j)=4;
elseif find(pop(i,j)<=M1+M2+M3+M4+M5);
j5=j5+1;
t5(1,j5)=j;
Calpop(i,j)=5;
elseif find(pop(i,j)<=M1+M2+M3+M4+M5+M6)
j6=j6+1;
t6(1,j6)=j;
Calpop(i,j)=6;
end
end
% R1
for i1=1:1:j1-1;
if abs(t1(1,i1)-t1(1,i1+1))==1
d11=d11+1;
end
end
% R2
for i2=1:1:j2-1;
if abs(t2(1,i2)-t2(1,i2+1))<3
d22=d22+3-abs(t2(1,i2)-t2(1,i2+1));
end
end
% R4
for i4=1:1:j4-1;
if abs(t4(1,i4)-t4(1,i4+1))<3
d44=d44+3-abs(t4(1,i4)-t4(1,i4+1));
end
end
% R5
for i5=1:1:j5-1;
if abs(t5(1,i5)-t5(1,i5+1))<21
d55=d55+21-abs(t5(1,i5)-t5(1,i5+1));
end
end
% R6
for i6=1:1:j6-1;
if abs(t6(1,i6)-t6(1,i6+1))<10
d66=d66+10-abs(t6(1,i6)-t6(1,i6+1));
end
end
% R1 & R2
for i12=1:1:j1;
for j12=1:1:j2
if abs(t1(1,i12)-t2(1,j12))<2
d12=d12+1;
end
end
end
% R2 & R4
for i24=1:1:j2;
for j24=1:1:j4
if abs(t2(1,i24)-t4(1,j24))<3
d24=d24+3-abs(t2(1,i24)-t4(1,j24));
end
end
end
% R2 & R6
for i26=1:1:j2;
for j26=1:1:j6
if abs(t2(1,i26)-t6(1,j26))<3
d26=d26+3-abs(t2(1,i26)-t6(1,j26));
end
end
end
% R4 & R6
for i46=1:1:j4;
for j46=1:1:j6
if abs(t4(1,i46)-t6(1,j46))<3
d46=d46+3-abs(t4(1,i46)-t6(1,j46));
end
end
end
% R5 & R6
for i56=1:1:j5;
for j56=1:1:j6
if abs(t5(1,i56)-t6(1,j56))<6
d56=d56+6-abs(t5(1,i56)-t6(1,j56));
end
end
end
pop(i,t)=d11+d12+d22+d24+d44+d55+d56+d26+d46+d66;
Calpop(i,t)=d11+d12+d22+d24+d44+d55+d56+d26+d46+d66;
end