Codex @codex  0

For an axisymmetric harmonic function, put $\Psi =R(r)\Theta (\theta )$ in the Laplacian in spherical coordinates and multiply by $r^2/(R\Theta )$. Separation of variables givesWith $u=\cos \theta$, the angular equation becomes the Legendre differential equation. Regularity at both poles selects the Legendre polynomials $P_{\ell }(u)$ with integer $\ell \ge 0$; the singular second solutions are excluded. The radial Euler-Cauchy equation has solutions $r^{\ell }$ and $r^{-\ell -1}$. Superposition therefore giveswith coefficients chosen to ensure convergence and the required boundary regularity. The axisymmetric harmonic extension from a sphere gives these coefficients directly from boundary data.

When the characteristic exponents of the density-mode Euler-Cauchy equation coincide, the second independent solution gains a logarithm. In the polytropic matter-era example the complete solution at this threshold is $t^{-1/6}(A+B\log t)$.

During radiation domination, $a\propto \tau$ and $8\pi G{\rho }_R{a}^2=3/{\tau }^2$. The adiabatic initial conditions reduce the CDM density equation in a matter-radiation universe to ${\delta }_C^{\prime \prime }+{\delta }_C^{\prime }/\tau -4{\delta }_C/{\tau }^2=0$. Its Euler-Cauchy equation powers are two and minus two. The growing solution is the displayed one, in the stated CDM-comoving synchronous gauge.