Fixed some bugs in the EM code.
Newton-Raphson takes half-steps when LL increases on a new guess. Improved EM stability and better estimates.
limit number of iterations in an optimization loop.