//Created by libAntimony v2.4 model bertram_2008__environment(time_) // Variable initializations: time_ = ; end model bertram_2008__membrane(V, Ica, Is1, Is2, Il, Ik, time_) // Rate Rules: V' = -(Ica + Is1 + Is2 + Il + Ik) / Cm; // Variable initializations: V = -40.0; Cm = 4525.0; Ica = ; Is1 = ; Is2 = ; Il = ; Ik = ; time_ = ; end model bertram_2008__m(m_infinity, V) // Assignment Rules: m_infinity := 1.0 / (1.0 + exp((vm - V) / sm)); // Variable initializations: vm = -22.0; sm = 7.5; V = ; end model bertram_2008__Ica(Ica, V, m_infinity) // Sub-modules, and any changes to those submodules: m: bertram_2008__m(m_infinity, V); // Assignment Rules: Ica := gCa * m_infinity * (V - VCa); // Variable initializations: gCa = 280.0; VCa = 100.0; end model bertram_2008__s1(s1, time_, V) // Assignment Rules: s1_infinity := 1.0 / (1.0 + exp((vs1 - V) / ss1)); // Rate Rules: s1' = autos1 * ((s1_infinity - s1) / tau_s1) + (1.0 - autos1) * (s1knot - s1); // Variable initializations: s1 = 0.9; autos1 = 1; s1knot = 1; tau_s1 = 1000.0; vs1 = -50.0; ss1 = 5; time_ = ; V = ; end model bertram_2008__Is1(Is1, time_, VK, V, s1) // Sub-modules, and any changes to those submodules: s1_mod: bertram_2008__s1(s1, time_, V); // Assignment Rules: Is1 := gs1 * s1 * (V - VK); // Variable initializations: gs1 = 22.0; VK = ; end model bertram_2008__n(n, time_, V) // Assignment Rules: n_infinity := 1.0 / (1.0 + exp((vn - V) / sn)); tau_n := tau_n_bar / (1.0 + exp((V - vn) / sn)); // Rate Rules: n' = (n_infinity - n) / tau_n; // Variable initializations: n = 0.0; tau_n_bar = 8.25; vn = -9.0; sn = 10.0; time_ = ; V = ; end model bertram_2008__Ik(Ik, VK, time_, V, n) // Sub-modules, and any changes to those submodules: n_mod: bertram_2008__n(n, time_, V); // Assignment Rules: Ik := gK * n * (V - VK); // Variable initializations: VK = -80.0; gK = 1300.0; end model bertram_2008__Il(Il, V, time_) // Assignment Rules: Il := gl * q * (V - Vl); q := (1.0 + p) / 2.0; beta_p := alpha_p * (1.0 / p0 - 1.0); sigma := root((alpha_p * (1.0 - p) + beta_p * p) / (tau_p * nstoc)); nstoc := 1000.0 / power(delNoise, 2.0); // Rate Rules: p' = (alpha_p * (1.0 - p) - beta_p * p) / tau_p + noise * w * sigma; // Variable initializations: gl = 41.0; Vl = -40.0; p = 0.14; alpha_p = 1.0; tau_p = 100.0; p0 = 0.2; noise = 1; w = 1; delNoise = 3; V = ; time_ = ; end model bertram_2008__s2(s2, time_, V) // Assignment Rules: s2_infinity := 1.0 / (1.0 + exp((vs2 - V) / ss2)); // Rate Rules: s2' = autos2 * ((s2_infinity - s2) / tau_s2) + (1.0 - autos2) * (s2knot - s2); // Variable initializations: s2 = 0.5; autos2 = 1; s2knot = 0.47; tau_s2 = 30000.0; vs2 = -40.0; ss2 = 15; time_ = ; V = ; end model bertram_2008__Is2(Is2, time_, VK, V, s2) // Sub-modules, and any changes to those submodules: s2_mod: bertram_2008__s2(s2, time_, V); // Assignment Rules: Is2 := gs2 * s2 * (V - VK); // Variable initializations: gs2 = 16; VK = ; end model *bertram_2008____main() // Sub-modules, and any changes to those submodules: environment: bertram_2008__environment(time_); membrane: bertram_2008__membrane(V, Ica0, Is10, Is20, Il0, Ik0, time_); Ica: bertram_2008__Ica(Ica0, V, m_infinity); Is1: bertram_2008__Is1(Is10, time_, VK, V, s1); Ik: bertram_2008__Ik(Ik0, VK, time_, V, n); Il: bertram_2008__Il(Il0, V, time_); Is2: bertram_2008__Is2(Is20, time_, VK, V, s2); end