For effective patient treatment in proton therapy, it is therefore important to accurately measure the beam range. For measuring beam range, various researchers determine the beam range by measuring the prompt gammas generated during nuclear reactions of protons with materials. However, the accuracy of the beam range determination can be lowered in heterogeneous phantoms, because of the differences with respect to the prompt gamma production depending on the properties of the material. In this research, to improve the beam range determination in a heterogeneous phantom, we derived a formula to correct the prompt-gamma distribution using the ratio of the prompt gamma production, stopping power, and density obtained for each material. Then, the prompt-gamma distributions were acquired by a multi-slit prompt-gamma camera on various kinds of heterogeneous phantoms using a Geant4 Monte Carlo simulation, and the deduced formula was applied to the prompt-gamma distributions. For the case involving the phantom having bone-equivalent material in the soft tissue-equivalent material, it was confirmed that compared to the actual range, the determined ranges were relatively accurate both before and after correction. In the case of a phantom having the lung-equivalent material in the soft tissue-equivalent material, although the maximum error before correction was 18.7 mm, the difference was very large. However, when the correction method was applied, the accuracy was significantly improved by a maximum error of 4.1 mm. Moreover, for a phantom that was constructed based on CT data, after applying the calibration method, the beam range could be generally determined within an error of 2.5 mm. Simulation results confirmed the potential to determine the beam range with high accuracy in heterogeneous phantoms by applying the proposed correction method. In future, these methods will be verified by performing experiments using a therapeutic proton beam.