Eq. (4) was solved by 4th-order Runge–Kutta numerical method with initial condition: [TG] = [TG]0 at t = 0 and was used to calculate Xkinetic.
A computer program was developed and applied for the minimization of OF (Eq. (6)).
The obtained reaction rate constant (k) in various reaction temperatures and their coefficients of determinations with comparison to RSM data are presented in Table 7.