S = 1   // S initial conc
 E = 1   // E initial conc
 ES = 0  // ES initial conc
 S + E -> ES; k1*S*E - k2*ES // enzyme-substrate binding
 k1 = 1  // forward rate
 k2 = .1 // backward rate