function mu=muestim(N,u,v)
    mu=zeros(1,2*N+1)
    Xg=zeros(1,5000);
    Xg(1)=N;
    for k=1:4999
    Xg(k+1)=grand(1,'bin',2*N,Xg(k)/(2*N)*(1-u)+(1-Xg(k)/(2*N))*(1-v));
    end
    for k=1:2*N+1
        mu(k)=sum(Xg==k-1)/5000
    end
endfunction