compute_sup.m

function sup=compute_sup(S,D)
 [n,m]=size(D)
 sup=0;
 for i=1:n
     if all(D(i,:)-S>=0)==true
         sup=sup+1;
     end
 end
 sup=sup/n;
end
  

compute_conf.m

function conf=compute_conf(R,D)
conf=compute_sup(abs(R),D)/compute_sup((abs(R)-R)/2,D);
end