function m = binary(n)
%m = (n.^2./2).*(2.^(-3.*(n-1)));
%m = 1-(3./n);
m = (n.^2-3.*n+log(4.*(n.^2))./log(2))./(n.^2);