Wrapper around `scipy.integrate.odeint` that solves any chain given initial conditions and decay constants.