(eng): Discrete Racah polynomials (DRPs) are highly efficient orthogonal polynomials widely used in various scientific fields for signal representation. They find applications in disciplines like image processing and computer vision. Racah polynomials were originally introduced by Wilson and later modified by Zhu to be orthogonal on a discrete set of samples. However, when the degree of the polynomial is high, it encounters numerical instability issues. In this paper, we propose a novel algorithm called Improved Stabilization (ImSt) for computing DRP coefficients. The algorithm partitions the DRP plane into asymmetric parts based on the polynomial size and DRP parameters. We have optimized the use of stabilizing conditions in these partitions. To compute the initial values, we employ the logarithmic gamma function along with a new formula. This combination enables us to compute the initial values efficiently for a wide range of DRP parameter values and large polynomial sizes. Additionally, we have derived a symmetry relation for the case when the Racah polynomial parameters are zero ($a=0$, $\alpha=0$, $\beta=0$). This symmetry makes the Racah polynomials symmetric, and we present a different algorithm for this specific scenario. We have demonstrated that the ImSt algorithm works for a broader range of parameters and higher degrees compared to existing algorithms. A comprehensive comparison between ImSt and the existing algorithms has been conducted, considering the maximum polynomial degree, computation time, restriction error analysis, and reconstruction error. The results of the comparison indicate that ImSt outperforms the existing algorithms for various values of Racah polynomial parameters.