“…Obviously, it is difficult to find the analytical solution, and it is more suitable to find the numerical solution by a computer. In this paper, the Runge–Kutta fourth-order method is used to calculate the amount of each component in each time step . The right side of eqs – is written as the function F i ( i = 1, 2, 3, 4, 5, 6, and 7). normald C Sa normald t = F 1 ( C Sa , C Ar , C Re , C As , C SO , C G , C W ) normald C Ar normald t = F 2 ( C Sa , C Ar , C Re , C As , C SO , C G , C W ) normald C Re normald t = F 3 …”