[loalive]

在你这个问题中，我几点让我感到疑惑的地方，可能是未表达清楚：
1、p为一m维向量，那（1-p）^(n-1)是不是应该表达成(1-p).^(n-1)?同样 an=p.*(1-p).^(n-1)
2、Pr是不是应该为Pn？
我便按照我的理解（上面提到的）写了下面的这个函数：

function P=plot_P_p(p,n)
pl=length(p);
a=[];
pm=zeros(1,pl);
P=[];
for m=1:n
am=[a p.*(1-p).^(m-1)];
pm=pm+am;
P=[P;pm];
end
figure;
plot(p,P);