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

  // Variable initializations:
  time_ = ;
end

model tran_model_2009__SERCA(time_)

  // Assignment Rules:
  kdatp := k_m1 / k_p1;
  T_ATP := ATP / kdatp;
  T_Cai := Ca_i / kdcai;
  T_Casr := Ca_sr / kdcasr;
  T_H1 := H_i / kdh1;
  T_Hi := power(H_i, n) / kdhi;
  T_Hsr := power(H_i, n) / kdhsr;
  T_H := H_i / kdh;
  a_p1 := (k_p2 * T_ATP * power(T_Cai, 2)) / (T_ATP * power(T_Cai, 2) + T_Hi * (1 + T_ATP * (1 + T_H1 + power(T_Cai, 2))));
  a_p2 := (k_p3 * T_Hsr) / (T_Hsr * (1 + T_H) + T_H * (1 + power(T_Casr, 2)));
  a_m1 := (k_m2 * ADP * power(T_Casr, 2) * T_H) / (T_Hsr * (1 + T_H) + T_H * (1 + power(T_Casr, 2)));
  a_m2 := (k_m3 * P_i * T_Hi) / (T_ATP * power(T_Cai, 2) + T_Hi * (1 + T_ATP * (1 + T_H1 + power(T_Cai, 2))));
  s2 := a_p1 + a_p2 + a_m1 + a_m2;
  v_cycle := (a_p1 * a_p2 - a_m1 * a_m2) / s2;

  // Rate Rules:
  ATP' = 0.1;

  // Variable initializations:
  k_p1 = 25900;
  k_p2 = 2540;
  k_p3 = 20.5;
  k_m1 = 2;
  k_m2 = 67200;
  k_m3 = 149;
  kdcai = 0.9;
  kdcasr = 2.24;
  kdh1 = 1.09e-5;
  kdhi = 3.54e-3;
  kdhsr = 1.05e-8;
  kdh = 7.24e-5;
  n = 2;
  Ca_i = 1e-3;
  Ca_sr = 1;
  H_i = 1e-4;
  ATP = 0.01;
  ADP = 20e-3;
  P_i = 1;
  time_ = ;
end

model *tran_model_2009____main()

  // Sub-modules, and any changes to those submodules:
  environment: tran_model_2009__environment(time_);
  SERCA: tran_model_2009__SERCA(time_);
end