Уравнение теплопроводности в Матлаб: Вычислить распространение температуры по треугольной поверхност

Задание: Вычислить распространение температуры по треугольной поверхности АВС (координаты точек: А(0;0), В(10;0), С(5;6)) с круглым отверстием с центром в точке О и радиусом R=0.04. Если Уравнение функции температуры имеет вид: dT/dt=d2T/dx2+d2T/dy2+q(T-50), где q=12 и граничные условия:
АВ:
dT/dy=0.
ВС: температура Т растет линейно от Т=0 в точке В до Т=50 в точке С.
СА: температура Т падает линейно от Т=50 в точке С к Т=15 в точке А.
На окружности отверстия Т=10.
Указания: запускать m-файл с названием "aT_In_Triangle".

Реализована задача теплопроводности в программной среде Матлаб