  Stable numerical simulations for a hyperbolic system of conservation laws of
relaxation type but not in divergence form are obtained by incorporating the
physical entropy into the simulations. The entropy balance is utilized as an
additional equation to eliminate the numerically critical terms with simple
substitutions. The method has potential for a wider applicability than the
particular example presented here.
