C1=[0.4 0.5
    0.5 0.4
    0.5 0.75
    0.8 0.1]  % positive class

C2=[0.4 0.5
    0.15  0.4
    0.25  0.75
    0.35  0.6
   ] %
   

minD=pointDist(X,M{i}{1}(1:2));
for i=1:NC1
    pDist(:,:,i)=pointDist(X,M{i}{1}(1:2));
    minD=min(minD,pDist(:,:,i));
end

minI=i*(minD==pDist(:,:,i));
for i=1:NC1
    mI(:,:,i)=i*(minD==pDist(:,:,i));
    minI=max(minI,mI(:,:,i)); 
end