Added functional rate laws to the model and fixed some errors in the
rate constants and initial conditions. The model now reproduces the
results from Tashima2006 for DDS=0.0 and 0.005.