added plots of initial guess for J, and for J from the formal solution
more minor fixes. getting closer, but still not there
fixed index error in calculating J
changed all array initializations from "empty" to "zeros" just to be safe...