ABSTRACT

For simplicity of coding we have used the “shiftlr” routine of Hansleman and Littlefield1. shiftlr (A,b) shifts the contents of matrix, A to the right by b spaces, if b > 0, and by b spaces to the left if b < 0. This is equivalent to multiplying a matrix by the appropriate Toepliz matrix:

If A = then we have:

(4.2)

and

The connect processes described by Equations (4.1) and/or (4.2) would cause data to “fall off the edge.” Lines 23 and 24 define short-circuit (zero-velocity) boundaries. Lines 26 to 28 fill one row of the display matrix (Z) at each iteration.