function dist = soft_delta_distrib(i, n, p)
% SOFT_DELTA_DISTRIB Make a discrete prob. distr. of length n with p of its mass on a single point i.
% dist = delta_dist(i, n, p)

dist = ((1-p)/(n-1))*ones(n,1);
%dist(i) = p;
rest = setdiff(1:n, i);
dist(i) = 1-sum(dist(rest));
