// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;

// SBML property definitions
property sbmlRole=string;
property sbmlName=string;
property sbmlCompartment=string;

// SBML reactions
// TC: EmptySet <=> M 
// TL: EmptySet <=> P 
// Md: M <=> EmptySet 
// Pd: P <=> EmptySet 

math main {
  realDomain time second;
  time.min=0;
  extern time.max;
  extern time.delta;

  // variable definitions
  real compartment_0000004 = 1E-15 L;
  real EmptySet = 0 M;
  real M(time) M;
  real P(time) M;
  real TC(time) mole/hr;
  real rM = 1;
  real k = 1;
  real n = 2;
  real TL(time) mole/hr;
  real rP = 1;
  real m = 3;
  real parameter_0000009 = 4;
  real Md(time) mole/hr;
  real qM = .21;
  real Pd(time) mole/hr;
  real qP = .21;

  // equations
  when (time=time.min) M = 1E-15/compartment_0000004;
  M:time = (TC + -1*Md)/compartment_0000004;
  when (time=time.min) P = 1E-15/compartment_0000004;
  P:time = (TL + -1*Pd)/compartment_0000004;
  TC = compartment_0000004*(rM/(1+(P/k)^n));
  TL = compartment_0000004*rP*M(if (time<(time.min+parameter_0000009)) time.min else time-(parameter_0000009))^m;
  Md = compartment_0000004*qM*M;
  Pd = compartment_0000004*qP*P;

  // variable properties
  compartment_0000004.sbmlRole="compartment";
  compartment_0000004.sbmlName="cytoplasm";
  EmptySet.sbmlRole="species";
  EmptySet.sbmlCompartment="compartment_0000004";
  M.sbmlRole="species";
  M.sbmlCompartment="compartment_0000004";
  P.sbmlRole="species";
  P.sbmlCompartment="compartment_0000004";
  TC.sbmlRole="rate";
  rM.sbmlRole="parameter";
  k.sbmlRole="parameter";
  n.sbmlRole="parameter";
  TL.sbmlRole="rate";
  rP.sbmlRole="parameter";
  m.sbmlRole="parameter";
  parameter_0000009.sbmlRole="parameter";
  parameter_0000009.sbmlName="tau";
  Md.sbmlRole="rate";
  qM.sbmlRole="parameter";
  Pd.sbmlRole="rate";
  qP.sbmlRole="parameter";
}