Как векторизовать зависимые петли в Matlab

Question

Мне было интересно, сможет ли кто-нибудь помочь мне векторизовать эту часть моего кода. Здесь bin_pdf - функция биномиального коэффициента. pb и pd - скалярные параметры. Спасибо!

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 

Answer 1

Это Векторизованных Choosenk. если вы не назовете его где-то еще, вы можете интегрировать его в bin_pdf.

function C=Choosenk(n,k) 
if k>n/2 
    k=n-k; 
end; 
C=prod((n-k+1:n)./(1:k)); 
end