求解：写一个产生多项分布随机数的程序，参数为n,size,prob,其中n为要产生的多项分布随机数个数，故程序返回一个n*length(prob)的矩阵，每一行表示一个多项随机数，并用此程序产生1000组样本，检查估计的频率分布和理论分布是否拟合。