ENH: moments for the zipf distribution

Should close gh-1886.
A remaining issue is that  `expect` reports the lack of convergence for small enough values of `a`. Which is a known issue (and the knownfail I'm leaving in this PR can serve as a useful smoke test if/when discrete `expect` is improved).