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

  // Variable initializations:
  time_ = ;
end

model komarova_2003__x1(x1, alpha1, beta1, g11, g21, x2, time_)

  // Rate Rules:
  x1' = alpha1 * power(x1, g11) * power(x2, g21) - beta1 * x1;

  // Variable initializations:
  x1 = 10.06066;
  alpha1 = ;
  beta1 = ;
  g11 = ;
  g21 = ;
  x2 = ;
  time_ = ;
end

model komarova_2003__x2(x2, alpha2, beta2, g12, g22, x1, time_)

  // Rate Rules:
  x2' = alpha2 * power(x1, g12) * power(x2, g22) - beta2 * x2;

  // Variable initializations:
  x2 = 212.132;
  alpha2 = ;
  beta2 = ;
  g12 = ;
  g22 = ;
  x1 = ;
  time_ = ;
end

model komarova_2003__z(z, k1, k2, y1, y2, time_)

  // Rate Rules:
  z' = k2 * y2 - k1 * y1;

  // Variable initializations:
  z = 100.0;
  k1 = ;
  k2 = ;
  y1 = ;
  y2 = ;
  time_ = ;
end

model komarova_2003__y1(y1, x1, x1_bar)

  // Assignment Rules:
  y1 := piecewise( x1 - x1_bar , x1 > x1_bar , 0.0 );

  // Variable initializations:
  x1 = ;
  x1_bar = ;
end

model komarova_2003__y2(y2, x2, x2_bar)

  // Assignment Rules:
  y2 := piecewise( x2 - x2_bar , x2 > x2_bar , 0.0 );

  // Variable initializations:
  x2 = ;
  x2_bar = ;
end

model komarova_2003__x1_bar(x1_bar, alpha1, beta1, alpha2, beta2, g22, g21, gamma)

  // Assignment Rules:
  x1_bar := power(beta1 / alpha1, (1.0 - g22) / gamma) * power(beta2 / alpha2, g21 / gamma);

  // Variable initializations:
  alpha1 = ;
  beta1 = ;
  alpha2 = ;
  beta2 = ;
  g22 = ;
  g21 = ;
  gamma = ;
end

model komarova_2003__x2_bar(x2_bar, alpha1, beta1, alpha2, beta2, g11, g12, gamma)

  // Assignment Rules:
  x2_bar := power(beta1 / alpha1, g12 / gamma) * power(beta2 / alpha2, (1.0 - g11) / gamma);

  // Variable initializations:
  alpha1 = ;
  beta1 = ;
  alpha2 = ;
  beta2 = ;
  g11 = ;
  g12 = ;
  gamma = ;
end

model komarova_2003__model_parameters(alpha1, beta1, alpha2, beta2, k1, k2, g11, g21, g12, g22, gamma)

  // Assignment Rules:
  gamma := g12 * g21 - (1.0 - g11) * (1.0 - g22);

  // Variable initializations:
  alpha1 = 3;
  beta1 = 0.2;
  alpha2 = 4;
  beta2 = 0.02;
  k1 = 0.24;
  k2 = 0.0017;
  g11 = 0.5;
  g21 = -0.5;
  g12 = 1;
  g22 = 0;
end

model *komarova_2003____main()

  // Sub-modules, and any changes to those submodules:
  environment: komarova_2003__environment(time_);
  x1: komarova_2003__x1(x10, alpha1, beta1, g11, g21, x20, time_);
  x2: komarova_2003__x2(x20, alpha2, beta2, g12, g22, x10, time_);
  z: komarova_2003__z(z0, k1, k2, y10, y20, time_);
  y1: komarova_2003__y1(y10, x10, x1_bar0);
  y2: komarova_2003__y2(y20, x20, x2_bar0);
  x1_bar: komarova_2003__x1_bar(x1_bar0, alpha1, beta1, alpha2, beta2, g22, g21, gamma);
  x2_bar: komarova_2003__x2_bar(x2_bar0, alpha1, beta1, alpha2, beta2, g11, g12, gamma);
  model_parameters: komarova_2003__model_parameters(alpha1, beta1, alpha2, beta2, k1, k2, g11, g21, g12, g22, gamma);
end