function [E] = c_2_p_10_mk_E(n)
%
inds = 1:n;
C = diag( inds .* ( n - inds + 1 ) );
for ii=1:n
  for jj=(ii+1):n
    C(ii,jj) = C(ii,jj-1) - ii;
    C(jj,ii) = C(ii,jj);
  end
end
E = (1/(n+1)) * C;