An efficient and positivity-preserving layer method is introduced to solve the radiation belt diffusion equation and is applied to study the bounce resonance interaction between relativistic electrons and magnetosonic waves. The layer method with linear interpolation, denoted by LM-L (layer method-linear), requires the use of a large number of grid points to ensure accurate solutions. We introduce a monotonicity- and positivity-preserving cubic interpolation method to be used with the Milstein-Tretyakov layer method. The resulting method, called LM-MC (layer method-monotone cubic), can be used to solve the radiation belt diffusion equation with a much smaller number of grid points than LM-L while still being able to preserve the positivity of the solution. We suggest that LM-MC can be used to study long-term dynamics of radiation belts. We then develop a 2-D LM-MC code and use it to investigate the bounce resonance diffusion of radiation belt electrons by magnetosonic waves. Using a previously published magnetosonic wave model, we demonstrate that bounce resonance with magnetosonic waves is as important as gyroresonance; both can cause several orders of magnitude increase of MeV electron fluxes within 1ﾠday. We conclude that bounce resonance with magnetosonic waves should be taken into consideration together with gyroresonance.