Modelling of the laser heating process is fruitful, since it enhances the understanding of the physical processes involved and minimizes the experimental cost. In the present study, an analytical solution for the temperature distribution inside the solid substrate is obtained using a Laplace transform method. A time exponentially decaying laser pulse profile is introduced in the analysis. The phase change process and recession velocity are accommodated to account for the evaporation at the surface. The closed-form solution obtained is compared with the analytical solution obtained previously for a conduction limited heating case. It is found that the closed-form solution obtained from the present study reduces to a previously obtained analytical solution when the pulse parameter, beta*, is set to zero in the closed-form solution. Temperature predictions from simulations agree well with the results obtained from the closed-form solution.