A method for getting a finite $\alpha$ in the IR region from an all-order beta function

The analytical method of QCD running coupling constant is extended to a model with an all-order beta function which is inspired by the famous Novikov-Shifman-Vai\-n\-s\-htein-Zakharov beta function of N=1 supersymmetric gau\-g\-e theories. In the approach presented here, the running coupling is determined by a transcendental equation with non-elementary integral of the running scale $\mu$. In our approach $\alpha_{an}(0)$, which reads 0.30642, does not rely on any dimensional parameters. This is in accordance with results in the literature on the analytical method of QCD running coupling constant. The new "analytically im\-p\-roved" running coupling constant is also compatible with the property of asymptotic freedom.