# MATLAB: How to code a function of the Probability mass function of the Poisson binomial distribution

functionsum

Hello,
I want to use a function that computes the Probability mass function of the Poisson binomial distribution \Pr(K=k)=\sum \limits {{A\in F{k}}}\prod \limits {{i\in A}}p{i}\prod \limits {{j\in A^{c}}}(1-p{j}) https://en.wikipedia.org/wiki/Poisson_binomial_distribution
But how do I implement this sum over subsets of {1,..,n} that contain k elements
Thank you, Daniel