In this research article, a model of the temperature distribution in the human skull is considered by using the Caputo fractional derivative. When interfering thermometry is lacking, these models are quite useful in estimating the temporal course of temperatures. The given model is numerically solved with the application of the polynomial least-squares scheme. We study various fractional-order cases to simulate the proposed phenomena more clearly compared to the previously proposed integer-order studies. The motivation behind this study is to involve memory effects in the model by using fractional time derivatives.