IndexError: index 31 is out of bounds for axis 1 with size 31

Question

I am using this code to solve a nonlinear PDE. Since I am using an implicit (backwards-difference) method I am required to solve a set of nonlinear algebraic equations. In m