Heat.m program