The Green function appropriate to a thermal resistance boundary condition is described.
The method of recursive images is presented to solve heat diffusion with a thermal contact resistance boundary condition.
Its actual numerical implementation, using MATLAB, is shown.
Four different heat diffusion problems are solved, illustrating how the method works.