Работа посвящена возможности применения матричного метода решения уравнения теплопроводности для многослойных сред в случае, когда в каком-либо слое возможен фазовый переход. Рассматриваются только стационарные процессы при отсутствии в среде внутренних источников тепла. Предложен общий метод решения для систем слоев, обладающих сдвиговой, осевой или центральной симметрией за счет использования аппарата обобщенных степеней Берса. По указанному методу проведены вычисления для одного вещества, когда в результате фазового перехода система становится двухслойной. Рассмотрена зависимость координаты точки фазового перехода от внешней температуры, проведено сравнение для сред с вышеперечисленными видами симметрии. Построено поле температур для многослойных сред с различными видами симметрии, когда в каком-либо слое произошел фазовый переход.