function P = eraseproc(M,N)
global m n
m = 2;
n = M + N;
E = 0;

P = zeros(m^n,1);
for xi = 1:m^n
  xy = i2c(xi);
  x = xy(1:M);
  y = xy(M+1:M+N);
  P(xi) = (1/2)^M * count_subseq(x,y) / choose(M,N);
end
return