the part of the code is
for k= 1:NumWin
X=mixmatr_ICA_st(:,:,k);
Y=mixmatr_ICA_st(:,:,k+1);
%Dimension= zeros(NumChan,NumChan,NumWin);
% D(NumChan,NumChan) = dist((mixmatr_ICA_st(:,:,k+1) -mixmatr_ICA_st(:,:,k))); % euclidean distance
D = pdist2(X,Y,'euclidean');
end