The solution to a linear propagation equation of the form L = 0 at each point in a computational lattice may be approximated by the superposition of the field generated by N hypothetical plane waves of strength c and with unit propagation in direction v . The solution at a discrete point x at the centre of a small computational stencil, as depicted in Fig. 1,is