This allows us to solve numerically for ρB(r, t + Δt) around r = 0 while using an analytical extension for long distances (long distance was defined as r > rP, rP such that PAB (r) = 10-6). Equation (18) is then split into a numerical and an analytical part: