mlzipf function

Zipf distribution maximum likelihood estimation