The Method of Discretization in Time (MDT) is a hybrid numerical technique intended to alleviate upfront the computational procedure of timedependent partial differential equations of parabolic type upfront. The MDT engenders a sequence of adjoint second order ordinary differential equations, wherein the space coordinate is the independent variable and time becomes an embedded parameter. Essentially, the adjoint second order ordinary differential equations are considered of “quasistationary” nature. In this work, the MDT is used for the analysis of unsteady heat conduction in regular bodies (large wall, long cylinder and sphere) accounting for nearly constant thermophysical properties, uniform initial temperature and surface heat flux. In engineering applications, the surface heat flux is customarily provided by electrical heating, radiative heating and pool fire heating. It is demonstrated that the approximate, semianalytical temperature solutions of the first adjoint “quasistationary” heat conduction equations using the first time jump are easily obtainable for each regular body. For enhanced acccuracy, regression analysis is applied to the deviations of the dimensionless surface temperature as a function of the dimensionless time for each regular body.