We considered solely the \(^{26}Al\)  isotope when calculating the heat from the radiogenic decay.
The radiogenic heat production rate (\(Q_{rad}\) ) is obtained by the following equations [1]}:
\(\begin{split}\lambda & = \frac{ln(2)}{t_{1/2}} \\q_{rad} & = TE\times \lambda \\Q_{rad} & = q_{rad}\, e^{(-t\times \lambda )} \,m_{rock}\, dt\end{split}\) 
