code: def _EStep(self, K, x, z_h1, pi_h, p_h):
          
          # E-Step:
          for i in range(self.n):
              for k in range(K):
                  denom_sum = 0
                  for k2 in range(K):
                      denom_sum += pi_h[k2] * EMMMixPLAggregator.f(x[i], p_h[k2])
                  z_h1[i][k] = (pi_h[k] * EMMMixPLAggregator.f(x[i], p_h[k])) / denom_sum
docs: Description:
              Internal function for computing the E-Step of the EMM algorithm.