In Fricke-agarose gels, an accurate determination of the spatial dose distribution is hindered by the diffusion of ferric ions. In this work, a model was developed to describe the diffusion process within gel samples of finite length and, thus, permit the reconstruction of the initial spatial distribution of the ferric ions. The temporal evolution of the ion concentration as a function of the initial concentration is derived by solving Fick's second law of diffusion in two dimensions with boundary reflections. The model was applied to magnetic resonance imaging data acquired at high spatial resolution (0.3 mm) and was found to describe accurately the observed diffusion effects.