“…On L 1 the total field is a sum of the scattered field and the incident wave: U (ρ, z, t) = U s (ρ, z, t) + U i (ρ, z, t), and one should remember that (25) is formulated for the scattered field and the problem (45) is being solved for the total field. Approximating the derivatives in the left side of the auxiliary initial-boundary value problem in (25) with the central differences, and using the secondorder one-sided differences to approximate the derivatives in the righthand side, an explicit three-step scheme that updates the grid function…”