RHSBZ:

--------------------------------------------------------------------------
   Model of a Belousov-Zhabotinskii reactor.
--------------------------------------------------------------------------
   Form:
   stateDot = RHSBZ( t, state, d )
--------------------------------------------------------------------------

   ------
   Inputs
   ------
   t           (1,1) Time
   state       (2,1) State [x;z]
   d           (1,1) data structure
                     .eps  (1,1) kC*B/k5*A
                     .q    (1,1) 2*k3*k4/k2*k5
                     .f    (1,1) Stoichiometric factor
   -------
   Outputs
   -------
   stateDot    (2,1) State derivative

--------------------------------------------------------------------------
   Reference: Scott, S.K. (1994.) Oscillations, Waves and Chaos
              in Chemical Reactions. pp. 27-39.
--------------------------------------------------------------------------