“…On the second day, we substitute the initial values u(0, x, y) = u 0 (x, y) and c(0, x, y) = c 0 (x, y) with u(1, x, y) = u 1 (x, y) and c(1, x, y) = c 1 (x, y), and then proceed to solve ( 5), ( 6) using this revised initial condition on the time interval t ∈ [1,2]. The identical procedure, employing equations such as (7) for the second day, extends to time intervals [2,3], [3,4].…”