Ordinary differential equations have recently been extended to fractional equations that are transformed using fractional differential equations. These fractional equations are believed to have high accuracy and low computational cost compared to ordinary differential equations. For the first time, this paper focuses on extending the nonlinear heat equations to a fractional order in a Caputo order. A new perturbation iteration algorithm (PIA) of the fractional order is applied to solve the nonlinear heat equations. Solving numerical problems that involve fractional differential equations can be challenging due to their inherent complexity and high computational cost. To overcome these challenges, there is a need to develop numerical schemes such as the PIA method. This method can provide approximate solutions to problems that involve classical fractional derivatives. The results obtained from this algorithm are compared with those obtained from the perturbation iteration method (PIM), the variational iteration method (VIM), and the Bezier curve method (BCM). All solutions are tested with numerical simulations. The study found that the new PIA algorithm performs better than the PIM, VIM, and BCM, achieving high accuracy and low computational cost. One significant advantage of this algorithm is that the solutions obtained have established that the fractional values of alpha, specifically α, significantly influencing the accuracy of the outcome and the associated computational cost.