This paper will develop a computer program based on application of the finite difference method for predicting temperature distribution in fruits and vegetables represented by real axisymmetric shape. Owing to the irregular shape of the fruits and vegetables, the body-fitted coordinate grid generation method is used to transform the original irregular shape into a regular computational domain. The transient heat transfer characteristics of fruits and vegetables during the cooling process will be computed on the computational domain. The numerical results are checked and compared with previous ones done by Pan and Bhowmik(1991)who use the finite element method to calculate the temperature distribution of fruits during the cooling process. Our numerical results show a good agreement with their results.