“…To deal with the above high-dimensionality problem, several numerical methods have been proposed in the literature to reduce the infinite system of differential equations into a loworder system. These can be broadly classified into kinetic lumping methods [19,20], global orthogonal collocation [21,22], method of moments [23][24][25][26][27][28], numerical fractionation methods [29][30][31], discrete weighted Galerkin [32][33][34], orthogonal collocation on finite elements [35,36], and sectional grid methods [15,18,37,38]. The above numerical methods are computationally complex and require special mathematical skills.…”