Solution
= Solution
Linear <eccentricity damping> adds $-z/\tau$ to the complex equation:
$$
\dot z=(iA-\tau^{-1})z+i\sum_k\nu_ke^{i(g_kt+\beta_k)}.
$$
Its solution is
$$
\boxed{z(t)=C e^{(iA-1/\tau)t}
+\sum_{k=1}^2
\frac{\nu_k}{g_k-A-i/\tau}e^{i(g_kt+\beta_k)}}.
$$
Unlike the undamped free eccentricity, the homogeneous term decays exponentially. The forced response survives, acquires a phase lag, and has finite amplitude
$$
\frac{|\nu_k|}{\sqrt{(g_k-A)^2+\tau^{-2}}}.
$$