I'm trying to build the "Naundorf model" for a single compartment with the following kinetics equations
Code: Select all
BREAKPOINT {
SOLVE scheme METHOD sparse
gHHO = gHHObar*O
ina = gHHO*(v-ena)
}
KINETIC scheme{
~C <-> O (alpha_a(v+K*J*O),beta_a(v+K*J*O))
~I <-> C (alpha_ci(v),beta_ci(v))
~O <-> I (1/tau_i,0)
CONSERVE C + O + I = 1
}
Is the use of state var O in the rate function (alpha/beta) permited? Does it raise integration method issues?
Does the following use of TMP help or make it worse?
Code: Select all
BREAKPOINT {
TMP = O
SOLVE scheme METHOD sparse
gHHO = gHHObar*O
ina = gHHO*(v-ena)
}
KINETIC scheme{
~C <-> O (alpha_a(v+K*J*TMP),beta_a(v+K*J*TMP))
~I <-> C (alpha_ci(v),beta_ci(v))
~O <-> I (1/tau_i,0)
CONSERVE C + O + I = 1
}
Albert.