sum(k*(k-1)*p^r*(1-p)^k*binomial(r+k-1,k),k=0..infinity)