I´m trying to adding some specifics modeling channels in neuron,
but when i run my .hoc, neuron return this erro:
Code: Select all
loading membrane mechanisms from nrnmech.dll
Additional mechanisms from files
k_mfv_s.mod na_mfv_s.mod soma_k_s.mod soma_na_s.mod
1
oc>at line 31 in file k_mfv_s.mod:
SOLVE states
Error at section location soma(0.5)
The increment in the independent variable is less than machine roundoff error
nrniv: scopmath library error
near line 37
{run()}
^
fadvance( )
advance( )
step( )
continuerun(300 )
and others
Code: Select all
TITLE MVF potassium current to soma MN
NEURON {
SUFFIX soma_k_s
USEION k READ ek WRITE ik
RANGE kgmax, gk, ik
}
UNITS {
(S) = (siemens)
(mV) = (millivolt)
(mA) = (milliamp)
}
PARAMETER {
kgmax = 35 (milisiemens/cm2)
}
ASSIGNED {
v (mV)
ek (mV)
ik (mA/cm2)
}
STATE { n }
BREAKPOINT {
SOLVE states
ik = kgmax*n*n*n*n*(v-ek)
}
INITIAL {
: Assume v has been constant for a long time
n = alpha_n(v)/( alpha_n(v) + beta_n(v) )
}
DERIVATIVE states {
: Computes state variable n at present v & t
n' = ((alpha_n(v)*(1-n)) - (beta_n(v)*n))
}
FUNCTION alpha_n(Vm (mV)) (/ms) {
UNITSOFF
alpha_n = (0.1+(-0.1*Vm))/(exp((Vm-10)/-10) -1)
UNITSON
}
FUNCTION beta_n(Vm (mV)) (/ms) {
UNITSOFF
beta_n = 0.2/( exp((Vm - 35)/70) - 0.04 )
UNITSON
}