We develop a numerical technique for solving the one-dimensional heat equation that combine classical and integral boundary conditions. The combined Laplace transform, high-precision quadrature schemes, and Stehfest inversion ...