Mathematical models for the description, in a quantitative way, of the damages induced on monuments by the action of specific pollutants are often systems of nonlinear, possibly degenerate, parabolic equations. Although some of the asymptotic properties of the solutions are known, one needs a numerical approximation scheme in order to have a quantitative forecast at any time of interest. In thi...