//Created by libAntimony v2.4 model yanagihara_model_1980__environment(time_) // Variable initializations: time_ = ; end model yanagihara_model_1980__membrane(V, time_, i_Na, i_K, i_Leak, i_s, i_h) // Rate Rules: V' = -(i_Na + i_K + i_Leak + i_s + i_h) / Cm; // Variable initializations: V = -62.83; Cm = 1; time_ = ; i_Na = ; i_K = ; i_Leak = ; i_s = ; i_h = ; end model yanagihara_model_1980__sodium_channel_m_gate(m, V, time_) // Assignment Rules: alpha_m := (1 * (V + 37)) / (-exp((V + 37) / -(10)) + 1); beta_m := 40 * exp((V + 62) / -(17.8)); // Rate Rules: m' = alpha_m * (1 - m) - beta_m * m; // Variable initializations: m = 0.047938; V = ; time_ = ; end model yanagihara_model_1980__sodium_channel_h_gate(h, V, time_) // Assignment Rules: alpha_h := 1.209E-3 * exp((V + 20) / -(6.534)); beta_h := 1 / (1 + exp((V + 30) / -(10))); // Rate Rules: h' = alpha_h * (1 - h) - beta_h * h; // Variable initializations: h = 0.95994; V = ; time_ = ; end model yanagihara_model_1980__sodium_channel(i_Na, time_, V, m, h) // Sub-modules, and any changes to those submodules: sodium_channel_m_gate: yanagihara_model_1980__sodium_channel_m_gate(m, V, time_); sodium_channel_h_gate: yanagihara_model_1980__sodium_channel_h_gate(h, V, time_); // Assignment Rules: i_Na := 0.5 * power(m, 3) * h * (V - 30); end model yanagihara_model_1980__potassium_channel_n_gate(n, V, time_) // Assignment Rules: alpha_n := (0.009 * 1) / (1 + exp(-(V + 3.8) / 9.71)) + 0.0006; beta_n := (-(2.25E-4) * (V + 40)) / (1 - exp((V + 40) / 13.3)); // Rate Rules: n' = alpha_n * (1 - n) - beta_n * n; // Variable initializations: n = 0.509494; V = ; time_ = ; end model yanagihara_model_1980__potassium_channel(i_K, time_, V, n) // Sub-modules, and any changes to those submodules: potassium_channel_n_gate: yanagihara_model_1980__potassium_channel_n_gate(n, V, time_); // Assignment Rules: i_K := (0.7 * n * (exp(0.0277 * (V + 90)) - 1)) / exp(0.0277 * (V + 40)); end model yanagihara_model_1980__leakage_current(i_Leak, V) // Assignment Rules: i_Leak := 0.8 * (1 - exp(-(V + 60) / 20)); // Variable initializations: V = ; end model yanagihara_model_1980__calcium_channel_d_gate(d, V, time_) // Assignment Rules: alpha_d := (0.0145 * (V + 35)) / (1 - exp(-(V + 35) / 2.5)) + (0.03125 * V) / (1 - exp(- V / 4.8)); beta_d := (-(4.21E-3) * (V - 5)) / (-exp((V - 5) / 2.5) + 1); // Rate Rules: d' = alpha_d * (1 - d) - beta_d * d; // Variable initializations: d = 0.000032; V = ; time_ = ; end model yanagihara_model_1980__calcium_channel_f_gate(f, V, time_) // Assignment Rules: alpha_f := (-(3.55E-4) * (V + 20)) / (-exp((V + 20) / 5.633) + 1); beta_f := (0.000944 * (V + 60)) / (1 + exp(-(V + 29.5) / 4.16)); // Rate Rules: f' = alpha_f * (1 - f) - beta_f * f; // Variable initializations: f = 1; V = ; time_ = ; end model yanagihara_model_1980__calcium_channel(i_s, time_, V, d, f) // Sub-modules, and any changes to those submodules: calcium_channel_d_gate: yanagihara_model_1980__calcium_channel_d_gate(d, V, time_); calcium_channel_f_gate: yanagihara_model_1980__calcium_channel_f_gate(f, V, time_); // Assignment Rules: i_s := 12.5 * (0.95 * d + 0.05) * (0.95 * f + 0.05) * (exp((V - 10) / 15) - 1); end model yanagihara_model_1980__hyperpolarization_activated_channel_q_gate(q, V, time_) // Assignment Rules: alpha_q := (0.00034 * (V + 100)) / (exp((V + 100) / 4.4) - 1) + 0.0000495; beta_q := (0.0005 * (V + 40)) / (1 - exp(-(V + 40) / 6)) + 0.0000845; // Rate Rules: q' = alpha_q * (1 - q) - beta_q * q; // Variable initializations: q = 0.010759; V = ; time_ = ; end model yanagihara_model_1980__hyperpolarization_activated_channel(i_h, time_, V, q) // Sub-modules, and any changes to those submodules: hyperpolarization_activated_channel_q_gate: yanagihara_model_1980__hyperpolarization_activated_channel_q_gate(q, V, time_); // Assignment Rules: i_h := 0.4 * q * (V + 45); end model *yanagihara_model_1980____main() // Sub-modules, and any changes to those submodules: environment: yanagihara_model_1980__environment(time_); membrane: yanagihara_model_1980__membrane(V, time_, i_Na, i_K, i_Leak, i_s, i_h); sodium_channel: yanagihara_model_1980__sodium_channel(i_Na, time_, V, m, h); potassium_channel: yanagihara_model_1980__potassium_channel(i_K, time_, V, n); leakage_current: yanagihara_model_1980__leakage_current(i_Leak, V); calcium_channel: yanagihara_model_1980__calcium_channel(i_s, time_, V, d, f); hyperpolarization_activated_channel: yanagihara_model_1980__hyperpolarization_activated_channel(i_h, time_, V, q); end