This contains code applying Crank-Nicolson to solve the heat equation with Neumann boundary conditions, along with an approximating scheme.