Debug output to file is working.
Current version - not yet completely working. For now, assume N is a power of two. The logic for arbitrary N is built into the code, but not yet tested or debugged. For certain values of N, a first order error appears on the right end of the domain, decaying from the right boundary into the domain. Yet for other values of N, the error seems to vanish.