A concise operator form of the Fokker-Planck equation agreeing with that proposed by Weizenecker [Phys. Med. Biol. 63, 035004 (2018)10.1088/1361-6560/aaa186] for the joint orientational distribution of the coupled physical and magnetodynamic rotational diffusion of a single-domain ferromagnetic nanoparticle suspended in a liquid is written from the postulated Langevin equations for the stochastic dynamics. Series expansion of its solution in a complete set yields, using the theory of angular momentum, differential-recurrence equations for statistical moments for coupled motion with uniaxial symmetry of the internal anisotropy-Zeeman energy of a nanoparticle. The numerical results via the matrix iteration method suggest that the susceptibility is adequately approximated by a single Lorentzian with peak frequency given by the inverse integral relaxation time and are discussed in relation to those of the well-known "egg model".