如何在Matlab中向量化相关的For循环

    for t=0:min(T,r)
        for n=0:r-t
            pp = (bin_pdf(n,pb,r-t) * bin_pdf(t,pd,min(T,r)));
            pt = pt + t/r * pp;
            pn = pn + n/r * pp;
            pc = pc + (r-t-n)/r * pp;
        end
    end

function p = bin_pdf(x,rho,n) 
if (x > n) || (n < 0) 
    p = 0; 
else 
    p = Choosenk(n,x) * rho^x * (1-rho)^(n-x); 
end

function C=Choosenk(n,k) 
if k>n/2 
    k=n-k; 
end; 
C=1; 
for i=0:k-1 
     C=C*(n-i)/(k-i); 
end 
end