Solution (source code)

= Solution

Linearizing the constraint gives $\delta\rho=-2\rho_0\delta a/a_0$. Since $\delta p=\delta\rho/3$, the <acceleration> equation gives $\delta\ddot a=2\Lambda\delta a$. Hence perturbations contain $e^{\sqrt{2\Lambda}t}$: the Einstein static universe is unstable.