code: def C_stock(self):
          
          return self._C_sys * (self._Q_sys / self._Q_stock).to(u.dimensionless)
docs: Return the required concentration of material in the stock given a
          reactor's desired system flow rate, system concentration, and stock
          flow rate.
  
          :return: Concentration of material in the stock
          :rtype: float