Computer-aided automatic computation of eigenvalues for multi-layer transient heat conduction problems