The analysis of heat conduction through a solid with heat generation leads to a linear matrix differential equation with separated boundary conditions. We present a symmetric second order exponential integrator for the ...