add println code to jacobian2by2 to see output
Jacobi growth process (hasn't been tested this rigorously, but seems to work, up to the answer being transposed and in reverse order in Julia compared to Matlab.)
draft of README for book demos
add code for chapter 21