Solution (source code)

= Solution

The homogeneous Euler equation removes the background force. Products of perturbations are second order, while differentiating the peculiar velocity inside the expanding flow produces the Hubble drag term. The first-order equation is therefore
$$
\boxed{m\bar n a(\dot{\delta\mathbf v}+H\delta\mathbf v)
=-\nabla\delta P-m\bar n\nabla\delta\Phi}.
$$