
function value=star_discrepancy(z)

T=linspace(0,1,1000);

N=size(z,1);

value=[];
for t=T
    iok=find(max(z)<=t);
    area_t=t^2;
    value(end+1)=abs(length(iok)/N-area_t);
end

value=max(value);