fitMultiBin function

Fitting the Multiplicative Binomial Distribution when binomial random variable, frequency, probability of success and theta parameter are given