Clarify / revisit boundary conditions in Step 4

[gforsyth]

There have been a few issues around the periodic boundary conditions in Step 4. At the very least, they need further clarification, especially if we treat u[0] == u[-1] (which we do, for the moment)