reci.record(&soma.ik(0.5))
recg.record(&soma.g_k_3st(0.5))
rect.record(&t)
Where k_3st is the mechanism name that has been inserted into the soma. The model kinetics are as follows where I'm providing values for gbar as well as all the a and z variables in the hoc file.
Code: Select all
STATE { c2 c1 o }
BREAKPOINT {
SOLVE states METHOD sparse
g = gbar*o
ik = (1e-4)*g*(v - ek)
}
INITIAL { SOLVE states STEADYSTATE sparse}
KINETIC states {
rates(v)
~c1 <-> o (k12,k21)
~c2 <-> c1 (k23,k32)
CONSERVE c2+c1+o=1
}
PROCEDURE rates(v(millivolt)) {
k12 = a12*exp(z12*v)
k21 = a21*exp(-z21*v)
k23 = a23*exp(z23*v)
k32 = a32*exp(-z32*v)
}