%function P=customprock(m0,n0,k)
function P=customprock(m0,n0)
global m n
m = m0;
n = n0;
%ww = [9/10 1/20 1/20];
kk = 1:n-1;
%ww = 2.^(-kk);
%ww = 1./(-kk);
ww = kk;
ww = ww/sum(ww);
%ww = ww(end:-1:1);
%ww = ones(1,n-1)/(n-1);

PP = zeros(m^n,1);
for k=1:n-1

  v = 1/(n-k);
  pk = v/2 + 1/2;
  
  P = zeros(m^n,1);
  for xi=1:m^n
    x = i2c(xi);
    b = (x(k)==x(n));
    p = (1/2)^(n-1)*(pk*b + (1-pk)*(1-b));
    P(xi) = p;
  end
  
  %PP = PP + P/(n-1);
  PP = PP + P*ww(k);
end
P = PP;
return