function mu=LoilimWF(N,pAB,pBA)
    mu=zeros(1,N+1)
    Xg=zeros(1,5000);
    Xg(1)=floor(N/2);
    for k=1:4999
        Xg(k+1)=grand(1,'bin',Xg(k),1-pAB) + grand(1,'bin',N-Xg(k),pBA);
    end
    for k=1:N+1
        mu(k)=sum(Xg==k-1)/5000
    end
endfunction