“…In this paper, a numerical scheme is proposed to solve the reaction diffusion problem: find a nonnegative function u defined on Ω, a smooth domain R d (d ≥ 1) and such that u t − ∆u m = αu This problem and analogous problems have been studied from a theoretical point of view by several authors: Aronson-Crandall-Peletier [2], Berryman-Holland [3], Friedman-Lacey [4], Friedman-McLeod [5], [6], Levine-Sacks [13], Sabinina [15], and Sacks [16], [17], [18]. M.-N.…”