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

  // Variable initializations:
  time_ = ;
end

model martinov_2000__Ado(Ado)

  // Variable initializations:
  Ado = 1;
end

model martinov_2000__Met(Met, time_)

  // Assignment Rules:
  Met := piecewise( 45 , ( geq(time_, 0)) && (time_ < 5 ), 52 , ( geq(time_, 5)) && (time_ < 15 ), 55 , ( geq(time_, 15)) && (time_ < 60 ), 52 , ( geq(time_, 60)) && (time_ < 75 ), 45 , geq(time_, 75) );

  // Variable initializations:
  time_ = ;
end

model martinov_2000__Hcy(Hcy, AdoHcy, Ado, K_AHC)

  // Assignment Rules:
  Hcy := (AdoHcy * K_AHC) / Ado;

  // Variable initializations:
  AdoHcy = ;
  Ado = ;
  K_AHC = ;
end

model martinov_2000__K_AHC(K_AHC)

  // Variable initializations:
  K_AHC = 0.1;
end

model martinov_2000__AdoMet(AdoMet, V_MET, V_GNMT, V_MATI, V_MATIII, time_)

  // Rate Rules:
  AdoMet' = (V_MATI + V_MATIII) - (V_MET + V_GNMT);

  // Variable initializations:
  AdoMet = 60;
  V_MET = ;
  V_GNMT = ;
  V_MATI = ;
  V_MATIII = ;
  time_ = ;
end

model martinov_2000__AdoHcy(AdoHcy, V_MET, V_GNMT, V_D, Ado, K_AHC, time_)

  // Rate Rules:
  AdoHcy' = ((V_MET + V_GNMT) - V_D) / (1 + K_AHC / Ado);

  // Variable initializations:
  AdoHcy = 3;
  V_MET = ;
  V_GNMT = ;
  V_D = ;
  Ado = ;
  K_AHC = ;
  time_ = ;
end

model martinov_2000__V_MATI(V_MATI, Met, AdoMet)

  // Assignment Rules:
  V_MATI := V_MATImax / (1 + (Km_MATI / Met) * (1 + AdoMet / Ki_MATI));

  // Variable initializations:
  Met = ;
  AdoMet = ;
  V_MATImax = 561;
  Km_MATI = 41;
  Ki_MATI = 50;
end

model martinov_2000__V_MATIII(V_MATIII, AdoMet, Met)

  // Assignment Rules:
  V_MATIII := V_MATIIImax / (1 + (Km1_MATIII * Km2_MATIII) / (power(Met, 2) + Met * Km2_MATIII));
  Km1_MATIII := 20000 / (1 + 5.7 * power(AdoMet / (AdoMet + 600), 2));

  // Variable initializations:
  AdoMet = ;
  Met = ;
  V_MATIIImax = 22870;
  Km2_MATIII = 21.1;
end

model martinov_2000__V_MET(V_MET, AdoMet, AdoHcy)

  // Assignment Rules:
  V_MET := V_METmax / (1 + Km1_MET / AdoMet + Km2_MET_A + (Km2_MET_A * Km1_MET) / AdoMet);
  Km1_MET := 10 * (1 + AdoHcy / 4);

  // Variable initializations:
  AdoMet = ;
  AdoHcy = ;
  V_METmax = 4544;
  Km2_MET_A = 10;
end

model martinov_2000__V_GNMT(V_GNMT, AdoMet, AdoHcy)

  // Assignment Rules:
  V_GNMT := ((V_GNMTmax / (1 + power(Km_GNMT / AdoMet, 2.3))) * 1) / (1 + AdoHcy / Ki_GNMT);

  // Variable initializations:
  AdoMet = ;
  AdoHcy = ;
  V_GNMTmax = 10600;
  Km_GNMT = 4500;
  Ki_GNMT = 20;
end

model martinov_2000__V_D(V_D, Hcy)

  // Assignment Rules:
  V_D := alpha_d * Hcy;

  // Variable initializations:
  alpha_d = 1333;
  Hcy = ;
end

model *martinov_2000____main()

  // Sub-modules, and any changes to those submodules:
  environment: martinov_2000__environment(time_);
  Ado: martinov_2000__Ado(Ado0);
  Met: martinov_2000__Met(Met0, time_);
  Hcy: martinov_2000__Hcy(Hcy0, AdoHcy0, Ado0, K_AHC0);
  K_AHC: martinov_2000__K_AHC(K_AHC0);
  AdoMet: martinov_2000__AdoMet(AdoMet0, V_MET0, V_GNMT0, V_MATI0, V_MATIII0, time_);
  AdoHcy: martinov_2000__AdoHcy(AdoHcy0, V_MET0, V_GNMT0, V_D0, Ado0, K_AHC0, time_);
  V_MATI: martinov_2000__V_MATI(V_MATI0, Met0, AdoMet0);
  V_MATIII: martinov_2000__V_MATIII(V_MATIII0, AdoMet0, Met0);
  V_MET: martinov_2000__V_MET(V_MET0, AdoMet0, AdoHcy0);
  V_GNMT: martinov_2000__V_GNMT(V_GNMT0, AdoMet0, AdoHcy0);
  V_D: martinov_2000__V_D(V_D0, Hcy0);
end