For the following parameter values: = 1, = 0.15, 0 = 1 mg/L, 0 = 3 mg/L, 0 = 0, = 5E 5, Numerically solve the above mentioned system of ODEs from = 0 to = 12 using MATLAB's ODE89 function.

The reaction leading to the formation of DBPs can generally be described by a second-order two- component model as follows: aA + bB pP where, A, B and P are the concentrations of Chlorine, NOM, and DBPs, respectively (mg/L); and a, b, and p, are stoichiometric reaction coefficients. The reaction rate is typically assumed to be first order in the concentration of chlorine (A) and a NOM (B), yielding the following system of ODES: dA/dt = -K*A*B = (b/a) * (dA/dt) -Y8*(K A + B) dp/dt -p/a) * (dA/dt) = Yp * (K + A+B) where, K is the rection rate constant (L/mg/sec). dB/dt = = = * The analytical solution for the concentration of A can be described as: A(t) = 1 - -PA exp(- (A. 1) + K + Bo *t] A - B./YB (YgA , where, Ao and Bo are the initial concentrations of Chlorine and NOM (mg/L)