The present program solves the heat transfer equations for spherical foods by a numerical procedure that uses the Crank-Nicolson method. The program is implemented for conditions with as few restrictions as possible to enhance versatility. The results agree with the analytical solutions for two particular cases of boundary conditions. The program also computes the thermal effect in the foodstuff from the temperature data obtained by the numerical program.