\[
  P^{(1)} =
    \left(T \omega  \dot{\omega}-\frac{p_1 (T-1)^2}{T} \right)
    + 2 T \omega  \dot{\omega} a
    + \frac{2 T \omega  \dot{\omega}}{1-T} xy,
\]