This paper proposes a method for calculating the Landau damping of a low-energy collective mode in a harmonically trapped Bose-Einstein condensate. Based on the divergence-free analytical solutions for ground-state wavefunction of the condensate and eigenvalues and eigenfunctions for thermally excited quasiparticles, obtained beyond Thomas-Fermi approximation, this paper calculates the coupling matrix elements describing the interaction between the collective mode and the quasiparticles. With these analytical results this paper evaluates the Landau damping rate of a monopole mode in a spherical trap and discusses its dependence on temperature, particle number and trapping frequency of the system.