Multi-wavelength observations of mature supernova remnants (SNRs), especially with recent advances in gamma-ray astronomy, make it possible to constrain energy distribution of energetic particles within these remnants. In consideration of the SNR origin of Galactic cosmic rays and physics related to particle acceleration and radiative processes, we use a simple one-zone model to fit the nonthermal emission spectra of three shell-type SNRs located within 2 degrees on the sky: RX J1713.7-3946, CTB 37B, and CTB 37A. Although radio images of these three sources all show a shell (or half-shell) structure, their radio, X-ray, and gamma-ray spectra are quite different, offering an ideal case to explore evolution of energetic particle distribution in SNRs. Our spectral fitting shows that (1) the particle distribution becomes harder with aging of these SNRs, implying a continuous acceleration process, and the particle distributions of CTB 37A and CTB 37B in the GeV range are harder than the hardest distribution that can be produced at a shock via the linear diffusive shock particle acceleration process, so spatial transport may play a role; (2) the energy loss timescale of electrons at the high-energy cutoff due to synchrotron radiation appears to be always a bit (within a factor of a few) shorter than the age of the corresponding remnant, which also requires continuous particle acceleration; (3) double power-law distributions are needed to fit the spectra of CTB 37B and CTB 37A, which may be attributed to shock interaction with molecular clouds.