ABSTRACT
In Section 2 we adapt the techniques used in [8, 9] to derive a numerical
discretization of (1.1)-(1.6) in a polygonal domain and in a bounded interval in
one space dimension, using a fixed spatial mesh and a uniform time step. We
obtain stability estimates in Section 3. We conclude with Section 4 in which we
display some computed results.