-
Notifications
You must be signed in to change notification settings - Fork 0
/
ka.mod
executable file
·63 lines (54 loc) · 1.13 KB
/
ka.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
TITLE Potasium Type A current for RD Traub, J Neurophysiol 89:909-921, 2003
COMMENT
Implemented by Maciej Lazarewicz 2003 ([email protected])
ENDCOMMENT
INDEPENDENT { t FROM 0 TO 1 WITH 1 (ms) }
UNITS {
(mV) = (millivolt)
(mA) = (milliamp)
}
NEURON {
SUFFIX ka
USEION k READ ek WRITE ik
RANGE gbar, ik
}
PARAMETER {
gbar = 0.0 (mho/cm2)
v ek (mV)
}
ASSIGNED {
ik (mA/cm2)
minf hinf (1)
mtau htau (ms)
}
STATE {
m h
}
BREAKPOINT {
SOLVE states METHOD cnexp
ik = gbar * m * m * m * m * h * ( v - ek )
}
INITIAL {
settables(v)
m = minf
m = 0
h = hinf
}
DERIVATIVE states {
settables(v)
m' = ( minf - m ) / mtau
h' = ( hinf - h ) / htau
}
UNITSOFF
PROCEDURE settables(v) {
TABLE minf, hinf, mtau, htau FROM -120 TO 40 WITH 641
minf = 1 / ( 1 + exp( ( - v - 60 ) / 8.5 ) )
mtau = 0.185 + 0.5 / ( exp( ( v + 35.8 ) / 19.7 ) + exp( ( - v - 79.7 ) / 12.7 ) )
hinf = 1 / ( 1 + exp( ( v + 78 ) / 6 ) )
if( v < -63 ) {
htau = 0.5 / ( exp( ( v + 46 ) / 5 ) + exp( ( - v - 238 ) / 37.5 ) )
}else{
htau = 9.5
}
}
UNITSON