There are many methods to solve second order differential equations
numerically. Methods based on variational principles have the
advantage that the energy is conserved but they might not be as
efficient as difference methods. Standard imbedded difference methods
are very efficient but, for general domains, there is always the risk
of instabilities. We will characterize these instabilities and discuss
methods to avoid them.