def compute_expectation_analytic(qgy):
Tk = qgy.Tk
ans = [1]
for k in range(1, len(Tk)):
EtT = qgy.E_tT_simple(0, k, qgy.phi_y(k), qgy.psi_y(k))
XtT = 1
ans.append(EtT * XtT)
return ans
def gen_psi_matx_list(qgy):
ans = []
for i in range(qgy.n):
ans.append(qgy.psi_y(i))
return ans