[EM algorithm] added first EM implementation for case M1.L1 (performance seems quite dependent on initial values)