We numerically solve a diffusion-type equation for the accumulated int
erdiffusion coefficient in an interacting system (nonideal ion-exchang
e reaction) by using a modified Crank-Nicolson algorithm. We obtain ex
act and approximate index profiles in such a way that bump effects fou
nd in experimental ionic concentrations in glass are explained as a re
sult of cation-cation interactions. (C) 1997 Optical Society of Americ
a.