, we have the following auxiliary function: After substituting (2) and (7) into (6), and drop optimization irrelevant terms in KL regularization, we have in a rational function form: are all positive polynomials of .
we can follow the two steps of EBW to derive the GT formulas for .