// Created by libAntimony v2.8.0
model *Mears1997_CRAC_PancreaticBetaCells()

  // Compartments and Species:
  compartment COMpartment_;
  species $V_membrane in COMpartment_, $n in COMpartment_, $jm in COMpartment_;
  species $Ca_er_Ca_equations in COMpartment_, $Ca_i in COMpartment_;

  // Assignment Rules:
  i_Ca := i_Ca_f + i_Ca_s;
  i_K := g_K*n*(V_membrane - V_K);
  i_K_ATP := g_K_ATP*(V_membrane - V_K);
  i_K_Ca := (g_K_Ca*Ca_i^5)/(Ca_i^5 + kdkca^5)*(V_membrane - V_K);
  i_CRAC := g_CRAC*r_infinity*(V_membrane - V_CRAC);
  i_leak := g_leak*(V_membrane - V_CRAC);
  n_infinity := 1/(1 + exp((-15 - V_membrane)/6));
  tau_n := 9.09/(1 + exp((V_membrane + 15)/6));
  jm_infinity := 1 - 1/(1 + exp((V_membrane + 53)/2));
  tau_j := 50000/(exp((V_membrane + 53)/4) + exp((-53 - V_membrane)/4)) + 1500;
  J_er_tot := (J_er_leak + J_er_IP3) - J_er_p;
  J_mem_tot := -f*(gamma*i_Ca + Jmp);
  g_K_ATP := piecewise(2000, time > 60000 && time < 660000, 150);
  i_Ca_f := g_Ca_f*m_f_infinity*(V_membrane - V_Ca);
  m_f_infinity := 1/(1 + exp((-20 - V_membrane)/7.5));
  i_Ca_s := g_Ca_s*m_s_infinity*(1 - jm)*(V_membrane - V_Ca);
  m_s_infinity := 1/(1 + exp((-16 - V_membrane)/10));
  r_infinity := 1/(1 + exp((Ca_er_Ca_equations - Ca_er_bar)/sloper));
  J_er_p := (verp*Ca_i^2)/(Ca_i^2 + kerp^2);
  a_infinity := 1/(1 + dact/Ca_i);
  b_infinity := IP3/(IP3 + dip3);
  h_infinity := 1/(1 + Ca_i/dinh);
  O := a_infinity^3*b_infinity^3*h_infinity^3*1;
  J_er_leak := perl*(Ca_er_Ca_equations - Ca_i);
  J_er_IP3 := O*(Ca_er_Ca_equations - Ca_i);
  Jmp := (vmp*Ca_i^2)/(Ca_i^2 + kmp^2);

  // Rate Rules:
  V_membrane' = -(i_Ca + i_K + i_K_ATP + i_K_Ca + i_CRAC + i_leak)/Cm;
  n' = (lambda_n*(n_infinity - n))/tau_n;
  jm' = (jm_infinity - jm)/tau_j;
  Ca_er_Ca_equations' = -J_er_tot/(lambda_er*sigma_er);
  Ca_i' = J_er_tot/lambda_er + J_mem_tot;

  // Species initializations:
  V_membrane = -61/COMpartment_;
  n = 0.0005/COMpartment_;
  jm = 0.12/COMpartment_;
  Ca_er_Ca_equations = 60/COMpartment_;
  Ca_i = 0.11/COMpartment_;

  // Compartment initializations:
  COMpartment_ = 1;

  // Variable initializations:
  Cm = 6158;
  lambda_n = 1.85;
  lambda_er = 250;
  sigma_er = 1;
  g_K = 3900;
  V_K = -70;
  Vn = -15;
  Sn = 6;
  g_Ca_f = 810;
  V_Ca = 100;
  Vm_f = -20;
  Sm_f = 7.5;
  g_Ca_s = 510;
  Vm_s = -16;
  Sm_s = 10;
  Vj = -53;
  Sj = 2;
  g_K_Ca = 1200;
  kdkca = 0.55;
  g_CRAC = 75;
  V_CRAC = 0;
  Ca_er_bar = 40;
  sloper = 3;
  g_leak = 0;
  verp = 0.24;
  kerp = 0.09;
  IP3 = 0;
  dact = 0.35;
  dinh = 0.4;
  dip3 = 0.2;
  perl = 0.003;
  kmp = 0.35;
  vmp = 0.08;
  gamma = 3.607e-006;
  f = 0.01;

  // Other declarations:
  var i_Ca, i_K, i_K_ATP, i_K_Ca, i_CRAC, i_leak, n_infinity, tau_n, jm_infinity;
  var tau_j, J_er_tot, J_mem_tot, g_K_ATP, i_Ca_f, m_f_infinity, i_Ca_s, m_s_infinity;
  var r_infinity, J_er_p, a_infinity, b_infinity, h_infinity, O, J_er_leak;
  var J_er_IP3, Jmp;
  const COMpartment_, Cm, lambda_n, lambda_er, sigma_er, g_K, V_K, Vn, Sn;
  const g_Ca_f, V_Ca, Vm_f, Sm_f, g_Ca_s, Vm_s, Sm_s, Vj, Sj, g_K_Ca, kdkca;
  const g_CRAC, V_CRAC, Ca_er_bar, sloper, g_leak, verp, kerp, IP3, dact;
  const dinh, dip3, perl, kmp, vmp, gamma, f;
end