I can’t write this computer simulation for you (at least not based on the data provided) but will instead explain a few basic relationships that govern the heating and cooling of objects. I hope this helps.
Consider a building an object that is composed of $n$ objects of masses $m_i$ with specific heat capacities $c_{p,i}$, then the building has an overall heat capacity given by:
$cM=\displaystyle\sum_{i=1}^{n} c_{p,i} m_i$.
1. Heating:
Understand $cM$ to be the total heat energy $\Delta Q$ needed to raise the building’s temperature by 1 degree Celsius:
$\Delta Q=cM\Delta T$ for $\Delta T = 1 C$.
If we add an amount of heat per unit of time $t$, that is $\frac{\Delta Q}{\Delta t}$, then we can write:
$\frac{\Delta Q}{\Delta t}=cM \frac{\Delta T}{\Delta t}$.....(Eq.1).
With $\frac{\Delta T}{\Delta t}$ the rate of temperature change in time.
2. Cooling:
Assuming that the building’s temperature $T$ is higher than that of its surroundings $T_O$, the building is constantly losing heat in accordance with Newton’s cooling law:
$\frac{\Delta Q}{\Delta t}=hA(T(t)-T_O)$.....Eq.2
Where $h$ the heat transfer coefficient, $A$ the total outside surface of the building and $T(t)$ the building’s temperature at any time $t$.
3. Thermostatic heating:
Most buildings require temperature to be kept near a desired set point, say $T_s$ (typically 18 to 20 degrees Celsius).
It should now be apparent that when $T(t) \geq T_s$ and the heating power $\frac{\Delta Q}{\Delta t}$ is still ‘on’, then in accordance with Eq.1, $T(t)$ will continue rising. But that is not a desired outcome. Instead it would be better to switch the heating power to ‘off’ when $T(t) \geq T_s$. This is of course the principle of thermostatic heating. When the building has cooled back down so a that $T(t) \leq T_s$, the thermostatic controller switches the heating power back on. This system requires no timed heating power program and also adapts well to changing values of $T_O$, the building's surrounding temperature.
At this set temperature the average power consumption is given by Eq.2:
$\frac{\Delta Q}{\Delta t}=hA(T_s-T_O)$.