For rarefied gases, classical hydrodynamic equations such as Euler equations and Navier-Stokes equations cannot accurately describe the gas dynamics. The kinetic models such as the Boltzmann equation and the Enskog equation are usually too expensive to solve numerically. This inspires people to search for models between the hydrodynamic models and kinetic models. We developed a systematic methodology to derive a class of intermediate models called regularized 13-moment models’’, which is suitable for slip and transitional flows. The method is implemented for general inverse power law models, which has never been considered before. In our implementation, the derivation of the models are automated by computer codes, and the numerical solver is also automatically generated by programs. The resulting models are tested for shock structure problems, and the results show its capability to capture the correct flow structure in strong nonequilibirum.