n=length(p);
q=p;
m=zeros(n-1,n);
for i=1:n-1
[q,l]=sort(q);
m(i,:)=[l(1:n-i+1),zeros(1,i-1)];
q=[q(1)+q(2),q(3:n),1];
end

这是程序的第一个FOR循环.确定了概率大小值的排列得到M矩正.但是具体每一行表示什么不是很清楚.高手指点下..