MLE for computing smoothing constant

MLE for computing smoothing constant - all credits to Wilker