Numerical computation of Neumann controllers for the heat equation on a finite interval