Displaying 1 result from an estimated 1 matches for "e_ca".
Did you mean:
e_a
2013 Feb 18
1
attempt to apply non-function
...t;- function(time, init, parms) {
with(as.list(c(init, parms)),{
m_Na <- function(v) 1 / (1 + 15 * exp(-0.769 * (v+34)) * ((1 -
exp(-0.05(v+6)/0.11 * v + 0.66))))
dv <- I_ext - (
+ g_Na * (m_Na(v))^3 * h_Na * (E_Na - v)
+ ((g_CaT * m_CaT * h_CaT) + (g_CaS * m_CaS)) * (E_Ca - v)
+ (g_K * m_K^4) * (E_K - v)
+ (g_KCa * m_KCa * h_KCa * (E_K - v))
+ (g_Af * h_Af + g_As * h_As) * m_A^3 * (E_K - v)
+ (g_L * (E_L - v)))
dh_Na <- k_h * (0.08 * (1 - h_Na) * exp(-0.125 * (v + 39)) - h_Na / (1
+ exp(-(0.2 * (2 + 8)))))
dm_Ca <- -C_Ca * (g_...