The shape parameter is a significant factor in determining the effectiveness of the interpolation of the Inverse Multi-quadratic Radial Basis Function (IMQ-RBF). This paper utilizes a simulated annealing algorithm to get the best form parameter for various interpolation challenges. The algorithm is evaluated based on interpolation precision, running time, and iteration count. Numerical experimental results show that, compared with other commonly used methods, the algorithm significantly improves interpolation accuracy while ensuring fast parameter selection.