function out = pihcm(center,data,U)
%PIHCM finds the total variation to the size of the data set in
%Hard C-means

%Carson KAO, May 29, 1997
mf = U.^2;	
dist = disthcm(center,data);
obj_fun = sum(sum((dist.^2).*mf));
out = obj_fun/size(data,1);
