A family of numerical schemes based on finite difference operators is
introduced for the computation of elastic waves. We use a displacement
-stress formulation of the model. After deriving some stability result
s, we give an analysis of the computational cost. Imposing an accuracy
criterion on the phase velocity, we derive the numerical parameters.
We also show that an optimum order of approximation exists for a given
precision.