We provide a fast iteration method to calculate the ion equilibrium temperature in an ultracold neutral plasma (UNP).The temperature as functions of electron initial temperature and ion density is obtained and compared with the recent UNP experimental data.The theoretical predictions agree with the experimental results very well.The calculated ion equilibrium temperature by this method can be applied to study the UNP expansion process more effectively.