Following Cicerone and Bowhill [1970a, b, 1971] a Monte Carlo method is adopted as a basis for the derivation of the photoelectron‐heat input into the ionospheric plasma. Considering that the size of the source elements must be small compared with the range over which photoelectrons dissipate their energy, a great number of Monte Carlo runs are required normally for the computation of the heating rates. This approach is modified in an attempt to minimize the computation time. The heat‐input distributions are computed for arbitrarily small source elements that are spaced apart at distances corresponding to the photoelectron dissipation range. By means of a nonlinear interpolation procedure their individual heating‐rate distributions are utilized to produce synthetic ones that fill the gaps between the Monte Carlo generated distributions. By varying these gaps and the corresponding number of Monte Carlo runs the accuracy of the results is tested to verify the validity of this procedure. It is concluded that this model can reduce the computation time by as much as an order of magnitude, thus improving the feasibility of including Monte Carlo calculations in self‐consistent ionosphere models.