function y=c(q,z)
p=1-q;
lambda=z.*q./p;
sum=1;
for k=0:max(z)
poisson=exp(-lambda); for i=1:k poisson=poisson.*lambda./i; end sum=sum-poisson.*(1-((q/p).^(z-k)));
end
y=array2table(sum);
end
poisson
poisson=exp(-lambda); for i=1:k poisson=poisson.*lambda./i; end sum=sum-poisson.*(1-((q/p).^(z-k)));
Best Answer