//Created by libAntimony v2.4
model lenbury_2005__environment(time_)

  // Variable initializations:
  time_ = ;
end

model lenbury_2005__x(x, y, beta, delta1, time_)

  // Rate Rules:
  x' = delta1 * exp(beta * (1.0 - power(y, 2.0))) - delta1 * x;

  // Variable initializations:
  x = 1.0;
  y = ;
  beta = ;
  delta1 = ;
  time_ = ;
end

model lenbury_2005__y(y, x, z, beta, delta2, time_)

  // Rate Rules:
  y' = delta2 * x * exp(beta * (1.0 - power(z, 2.0))) - delta2 * y;

  // Variable initializations:
  y = 0.01;
  x = ;
  z = ;
  beta = ;
  delta2 = ;
  time_ = ;
end

model lenbury_2005__z(z, y, delta3, time_)

  // Rate Rules:
  z' = delta3 * y - delta3 * z;

  // Variable initializations:
  z = 0.01;
  y = ;
  delta3 = ;
  time_ = ;
end

model lenbury_2005__model_parameters(delta1, delta2, delta3, beta)

  // Variable initializations:
  delta1 = 0.5;
  delta2 = 0.38;
  delta3 = 0.6;
  beta = 1.091;
end

model *lenbury_2005____main()

  // Sub-modules, and any changes to those submodules:
  environment: lenbury_2005__environment(time_);
  x: lenbury_2005__x(x0, y0, beta, delta1, time_);
  y: lenbury_2005__y(y0, x0, z0, beta, delta2, time_);
  z: lenbury_2005__z(z0, y0, delta3, time_);
  model_parameters: lenbury_2005__model_parameters(delta1, delta2, delta3, beta);
end