dcmri.conc_mmcomp#

dcmri.conc_mmcomp(J, Vmax, Km, t=None, dt=1.0, solver='SM')[source]#

Indicator concentration inside a Michaelis-Menten compartment.

See section Michaelis-Menten compartment for more detail.

Parameters:
  • J (array_like) – the indicator flux entering the compartment.

  • Vmax (float) – Limiting rate in the same units as J. Must be non-negative.

  • Km (float) – Michaelis-Menten constant in units of concentration (or flux x time). Must be non-negative.

  • t (array_like, optional) – the time points of the indicator flux J, in the same units as Km/Vmax. If t=None, the time points are assumed to be uniformly spaced with spacing dt. Defaults to None.

  • dt (float, optional) – spacing between time points for uniformly spaced time points, in the same units as Km/Vmax. This parameter is ignored if t is explicity provided. Defaults to 1.0.

  • solver (str, optional) – choose which solver to use. The options are ‘SM’ for the Schnell and Mendoza solution or ‘prop’ for a numerical solution by forward propagation. Defaults to ‘SM’.

Returns:

Concentration as a 1D array.

Return type:

numpy.ndarray

Raises:

ValueError – if one of the parameters is out of bounds.

See also

flux_mmcomp

Example

>>> import dcmri as dc
>>> t = [0,5,15,30,60]
>>> J = [1,2,3,3,2]
>>> Vmax, Km = 1, 12
>>> dc.conc_mmcomp(J, Vmax, Km, t)
array([  0.        ,   7.5       ,  29.26723718,  64.27756059,
114.97656637])