Fixed bug which would cause non-convergence of the solver in the
second timestep and reshaped/resized vectors and matrices for
Von Neumann boundary conditions.