During decay, if we neglect the triple correlation scalars, then the equations reduce to axially symmetric five-dimensional heat equations,
Solutions to these five-dimensional heat equation was solved by Chandrasekhar. The initial conditions can be expressed in terms of Gegenbauer polynomials (without loss of generality),
where  are Gegenbauer polynomials. The required solutions are
where  is the Bessel function of the first kind. 
As  the solutions become independent of 
where