-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathkslowtab.mod
64 lines (56 loc) · 1.08 KB
/
kslowtab.mod
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
TITLE HH slow potassium channel with FUNCTION_TABLEs (.inf and .tau files)
: Implemented in Rubin and Cleland (2006) J Neurophysiology
: Parameters from Bhalla and Bower (1993) J Neurophysiology
: Adapted from /usr/local/neuron/demo/release/nachan.mod - squid
: by Andrew Davison, The Babraham Institute [Brain Res Bulletin, 2000]
NEURON {
SUFFIX kslowtab
USEION k READ ek WRITE ik
RANGE gkbar, ik
GLOBAL ninf, kinf, ntau, ktau
}
UNITS {
(mA) = (milliamp)
(mV) = (millivolt)
}
INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}
PARAMETER {
v (mV)
dt (ms)
gkbar= 0.120 (mho/cm2) <0,1e9>
ek = -70 (mV)
}
STATE {
n k
}
ASSIGNED {
ik (mA/cm2)
ninf
kinf
ntau (ms)
ktau (ms)
}
INITIAL {
rates(v)
n = ninf
k = kinf
}
BREAKPOINT {
SOLVE states METHOD cnexp
ik = gkbar*n*n*k*(v - ek)
}
DERIVATIVE states {
rates(v)
n' = (ninf - n)/ntau
k' = (kinf - k)/ktau
}
FUNCTION_TABLE tabninf(v(mV))
FUNCTION_TABLE tabntau(v(mV)) (ms)
FUNCTION_TABLE tabkinf(v(mV))
FUNCTION_TABLE tabktau(v(mV)) (ms)
PROCEDURE rates(v(mV)) {
ninf = tabninf(v)
ntau = tabntau(v)
kinf = tabkinf(v)
ktau = tabktau(v)
}