i = 1:10;
const = sum(1./i);
pmf1 = (1./i)./const;
pmf2 = (1-(1./i))./const;
plot(i,pmf1, 'k');
hold on
plot(i,pmf2, 'k--');