Consistency and stability for some nonnegativity-conserving methods (Q1318376)

The paper deals with a one-dimensional diffusion equation \(\partial u/\partial t- \text{div}(d(x) \nabla u)+ q(x)u= f(x,t)\) with \(d>0\) and \(g,f\geq 0\) and such discretizations of it which paraphraze the behavior of the original equation in the sense that the solution, being nonnegative at a current time level, remains nonnegative also at any further time level. The so-called ``diagonally splitted'' and ``diagonally splitted implicit'' methods are of a particular importance, the latter one being essentially identified as a combination of the former one in place of a predictor and the implicit Euler method in place of a corrector. Numerical experiments are included, too.