function  [J]= kmeans_evaluation3(y, class_label, center)

%
%
%
[nr, nc] = size(y);
[k, c] = size(center);

for i = 1:1:k
    t = 1;
    for j = 1:1:nr
        if i == class_label(j)
            x(t,1:c) = y(j,1:c);
            t = t+1;
        end
    end
    
    [xr,xc] = size(x);
    d(i) = 0;
    
    for j = 1:1:xr
        d(i) = d(i) + (x(j,:)-center(i,:))*(x(j,:)-center(i,:))';
    end
    clear x
end
J = sum(d)/nr;