clear all
close all

DATA = rand(300000,21);
rows = size(DATA,1);
cols = size(DATA,2);

C = mean(DATA);

rozdily = DATA - repmat(C,[rows 1]);
rozdily = rozdily.^2;
sumactvercu = sum(rozdily,2);

m = mean(sumactvercu);
s = std(sumactvercu);

hornimez = m+4*s;
dolnimez = m-4*s;

suspektni = find(sumactvercu>hornimez | sumactvercu<dolnimez);

anomalnidata = [suspektni DATA(suspektni,:)];