\[\zeta_t=0.5\left(\phi_2 + \sqrt{\phi_2^2 + 4 \rho_t\phi_1\phi_2}\right)\]