Just to see whether increasing region 1 radius makes the comparison with figure 2 of reference (Minervini et al., 1999) better.
For the doped system, a series of calculations of
and
was
performed at increasing region 1 radius. Region 2a radius was set at