code: def standard_kinetics(target, quantity, prefactor, exponent):
      r
      X = target[quantity]
      A = target[prefactor]
      b = target[exponent]
  
      r = A*(X**b)
      S1 = A*b*(X**(b - 1))
      S2 = A*(1 - b)*(X**b)
      values = {'S1': S1, 'S2': S2, 'rate': r}
      return values
docs: r"""