// 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
// Reaction1: EmptySet <=> z 
// Reaction2: u <=> EmptySet 
// Reaction3: z <=> u 

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

  // variable definitions
  real cell = 1 L;
  real kappa = .015;
  real k6 = 1;
  real k4 = 180;
  real k4prime = .018;
  real alpha(time);
  real EmptySet = 1 mol;
  real u(time) mol;
  real z(time) mol;
  real v(time) mol;
  real Reaction1(time) katal;
  real Reaction2(time) katal;
  real Reaction3(time) katal;

  // equations
  alpha = k4prime/k4;
  when (time=time.min) u = 0;
  u:time = k4*(v-u)*(alpha+u^2)-k6*u;
  z = v-u;
  when (time=time.min) v = 0;
  v:time = kappa-k6*u;
  Reaction1 = kappa;
  Reaction2 = k6*u;
  Reaction3 = k4*z*(k4prime/k4+u^2);

  // variable properties
  cell.sbmlRole="compartment";
  kappa.sbmlRole="parameter";
  k6.sbmlRole="parameter";
  k4.sbmlRole="parameter";
  k4prime.sbmlRole="parameter";
  alpha.sbmlRole="parameter";
  EmptySet.sbmlRole="species";
  EmptySet.sbmlCompartment="cell";
  u.sbmlRole="species";
  u.sbmlCompartment="cell";
  z.sbmlRole="species";
  z.sbmlCompartment="cell";
  v.sbmlRole="species";
  v.sbmlCompartment="cell";
  Reaction1.sbmlRole="rate";
  Reaction2.sbmlRole="rate";
  Reaction3.sbmlRole="rate";
}