若用matlab
设a(k)=C(n,k)*(p^k)*(q^(n-k)), s(k)=a(1)+a(2)+…+a(k)
coef = 1;
ak1 = 1;
ak2 = q^n;
step = k;
s = 0;
for i = 1 : step
coef = (step-i+1)/i * step;
ak1 = ak1 * p;
ak2 = ak2 / q;
ak = coef * ak1 * ak2;
s = s + ak;
end
fprintf(‘a(1)+a(2)+…+a(k)= %fn’, s);