AbstractThis paper discusses three numerical integration schemes connected with a finite element solution to the transient heat‐conduction equation. The main emphasis is on the achievable accuracy. A simple recurrence relation derived from the Galerkin process is suggested for the treatment of fast‐varying boundary conditions as being significantly better than the usual Crank‐Nicholson scheme in matters of short‐time a