Reversible Series Reactions see equations 2.10

The solution vector sol and its elements (sol [1,1], e.g.) are too long to present here. They are functions of the rate constants (kl, k2, k3, and k4) and the independent variable t. One can obtain the solution from the CD with this book by changing the colons to semicolons at the end of the command line in the worksheet. As in Example 2.1, the solution obtained can be plotted for a particular set of values for the rate constants. Let s set these to be the same as those used on page 282 of Constantinides and Mostoufi. [6] 

Obtain the concentrations as functions of time by substituting these values for the rate constants into the components of the solution vector  

We can find an equation in terms of the rate constants that will provide a means of finding the time at which the maximum in Cb occurs. We do this by using Maple s simplify command and by finding dcbldi by using Maple s diff command. 

we use Maple s solve command to set the equation (eq) equal to zero and solve for tmax  

The concentration of B at tmax can be obtained by substitution of the tmax expression into the solution for the concentration of B as a function of time. 

---