Запропоновано спосіб спільного використання методів розщеплення, граничних елементів, покрокової часової схеми та ітераційної FD (Finite-Discrete) процедури для побудови інтегрального зображення розв’язку нестаціонарної задачі теплопровідності для замкненої області з заданою на межі умовою Діріхле, всередині якої знаходиться локально-неоднорідна підобласть з залежними від координат фізичними характеристиками. Виконано різносторонній чисельний аналіз цього підходу з урахуванням впливу на теплове поле залежностей від координат коефіцієнта теплопровідності та питомої теплоємності матеріалу.