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.