Equations
pi=sigma(-E_pi-G_pi)
s_{t+1}=B_t*s_t
gamma_A=1/beta_A
G_pi=sum_t[o_{pi, t} \dot (ln(o_{pi, t})-C)-diag(A \dot lnA) \dot s^{bar}_{pi, t}]
F_pi=sum_t[s^{bar}{pi, t} \dot (ln{s^{bar}{pi, t})-ln(A) \dot o_t - 0.5 * ln(B_{t-1}*s^{bar}{pi, t-1})-0.5 * ln(B{t+1}*s^{bar}_{pi, t+1})]
pi^{bar}=sigma(-E_pi-G_pi-F_pi)
s^{bar}{pi,t}=sigma(ln(B{pi, t-1}s_{t-1})+ln(A^{bar} \dot o^{bar}t)+ln(B{pi, t}*s_{t+1}))
beta^{bar}_A=beta_A - sum_t(ln(A) \dot (o^{bar}_t*s_t-A^{bar}*s_t))
o^{bar}_t=delta(o_t)