nFeats = 39;
nn1 = zeros(1,nFeats);
nn3 = zeros(1,nFeats);
nn5 = zeros(1,nFeats);
CentroiD = zeros(1,nFeats);
averageCent = zeros(1,nFeats);
OurMethod = zeros(1,nFeats);
for ind = 2:40
    fileName = ['numFeatures-' int2str(ind)];
    load(fileName);
    nn1(1,ind-1) = oneNN(1,ind-1);
    nn3(1,ind-1) = threeNN(1,ind-1);
    nn5(1,ind-1) = fiveNN(1,ind-1);
    CentroiD(1,ind-1) = cent(1,ind-1);
    averageCent(1,ind-1) =  avgcent(1,ind-1);
    OurMethod(1,ind-1) = centparam(1,ind-1);
end
qwe = 1:39;
plot(qwe+1,nn1,'-.k',qwe+1,nn3,'+k',qwe+1,nn5,'--k',qwe+1,CentroiD,'*k',qwe+1,averageCent,'^k',qwe+1,OurMethod, '-k');