function xie = xie_matrix(Z,X,U,m,n,k)
temp1 = zeros(m,m);
temp2 = 0;
for j=1:k
    for i=1:n
        temp1 = temp1 + Z(i,j).*(X(i,:)-U(j,:))'*(X(i,:)-U(j,:));
        temp2 = temp2 + Z(i,j);
    end
    xie(:,:,j) = temp1/temp2;
    temp1 = zeros(m,m);
    temp2 = 0;
end