Abstract: 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 results, 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. \\