% warm up

for run = 1:100
    dv = sampling(GaussianTest0p1, 'centroid', 'centroid', 100);
end

% Gaussian

fprintf('\n\n')
fprintf('Test set: %s \n', 'Gaussian');

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(GaussianTest0p1, 'centroid', 'centroid', 3);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 3);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(GaussianTest0p1, 'centroid', 'centroid', 5);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 5);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(GaussianTest0p1, 'centroid', 'centroid', 10);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 10);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(GaussianTest0p1, 'centroid', 'centroid', 100);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 100);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(GaussianTest0p1, 'centroid', 'centroid', 1000);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 1000);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(GaussianTest0p1, 'centroid', 'centroid', 10000);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 10000);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:10
    dv = sampling(GaussianTest0p1, 'centroid', 'centroid', 100000);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 100000);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

% Scurve

fprintf('\n\n')
fprintf('Test set: %s \n', 'Scurve');

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(ScurveTest0p1, 'centroid', 'centroid', 3);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 3);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(ScurveTest0p1, 'centroid', 'centroid', 5);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 5);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(ScurveTest0p1, 'centroid', 'centroid', 10);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 10);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(ScurveTest0p1, 'centroid', 'centroid', 100);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 100);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(ScurveTest0p1, 'centroid', 'centroid', 1000);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 1000);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(ScurveTest0p1, 'centroid', 'centroid', 10000);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 10000);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:10
    dv = sampling(ScurveTest0p1, 'centroid', 'centroid', 100000);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 100000);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

% shopping

fprintf('\n\n')
fprintf('Test set: %s \n', 'shopping');

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(shoppingTest0p1, 'centroid', 'centroid', 3);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 3);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(shoppingTest0p1, 'centroid', 'centroid', 5);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 5);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(shoppingTest0p1, 'centroid', 'centroid', 10);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 10);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(shoppingTest0p1, 'centroid', 'centroid', 100);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 100);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(shoppingTest0p1, 'centroid', 'centroid', 1000);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 1000);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(shoppingTest0p1, 'centroid', 'centroid', 10000);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 10000);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:10
    dv = sampling(shoppingTest0p1, 'centroid', 'centroid', 100000);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 100000);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

% powder

fprintf('\n\n')
fprintf('Test set: %s \n', 'powder');

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(powderTest0p1, 'centroid', 'centroid', 3);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 3);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(powderTest0p1, 'centroid', 'centroid', 5);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 5);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(powderTest0p1, 'centroid', 'centroid', 10);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 10);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(powderTest0p1, 'centroid', 'centroid', 100);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 100);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(powderTest0p1, 'centroid', 'centroid', 1000);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 1000);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:100
    dv = sampling(powderTest0p1, 'centroid', 'centroid', 10000);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 10000);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);

fprintf('\n')
DVlist = [];
tic
for run = 1:10
    dv = sampling(powderTest0p1, 'centroid', 'centroid', 100000);
    DVlist = [DVlist dv];
end
toc
meanDV = mean(DVlist);
stdDV = std(DVlist);
fprintf('Sample size: %0.0f \n', 100000);
fprintf('Mean defuzzified value: %0.10f \n', meanDV);
fprintf('Standard deviation: %0.10f \n', stdDV);
