如何求解下面的一个概率分布问题？ 第1页

       Redball <- function(M, m, N) {     ball = c(rep(1,m), rep(0,M-m))     ball = sample(ball, M)     List = list()     order = 1:M     k = 0     for(i in 1:(M/N))     {         B = sample(order,N)         order = order[-B]         List[[i]] = ball[B]         if(max(ball[B])==1)k = k + 1     }     list(组 = List, 含红球组数 = k) }  pro.Redball <- function(M, m, N, n)    #n为实验重复次数 {     g = M/N + 1     b = rep(0,g)     for(i in 1:n)     {         num = Redball(M,m,N)[[2]] + 1         b[num] = b[num] + 1     }     p = b/n     g = g - 1      barplot(p, names.arg = paste(0:g), col = rgb(0.3,0.5,0.7))     p } pro.Redball(24, 5, 6, 1000)