a=fopen('momentKrawtchoukBD.txt','r')

A=fscanf(a,'%g',[6 Inf])

A=A'

taux_reconnaissance=0

b=fopen('momentKrawtchoukBDTest.txt','r')

B=fscanf(b,'%g',[6 Inf])

B=B'

 d=zeros(407,627);

 for i=1:627

for j=1:6

    for k=1:407
        for l=1:6
          d(k,i)=d(k,i)+((A(i,j)-B(k,l))^2)
        end
    end
end

end