M4B function

Model M4B MLE function