:Comment :
:Reference : : Reuveni, Friedman, Amitai, and Gutnick, J.Neurosci. 1993
NEURON {
SUFFIX Ca
USEION ca READ eca WRITE ica
RANGE gCabar, gCa, ica
}
UNITS {
(S) = (siemens)
(mV) = (millivolt)
(mA) = (milliamp)
}
PARAMETER {
gCabar = 0.00001 (S/cm2)
}
ASSIGNED {
v (mV)
eca (mV)
ica (mA/cm2)
gCa (S/cm2)
mInf
mTau
mAlpha
mBeta
hInf
hTau
hAlpha
hBeta
}
STATE {
m
h
}
BREAKPOINT {
SOLVE states METHOD cnexp
gCa = gCabar*m*m*h
ica = gCa*(veca)
}
DERIVATIVE states {
rates()
m' = (mInfm)/mTau
h' = (hInfh)/hTau
}
INITIAL{
rates()
m = mInf
h = hInf
}
PROCEDURE rates(){
UNITSOFF
if((v == 27) ){
v = v+0.0001
}
mAlpha = (0.055*(27v))/(exp((27v)/3.8)  1)
mBeta = (0.94*exp((75v)/17))
mInf = mAlpha/(mAlpha + mBeta)
mTau = 1/(mAlpha + mBeta)
hAlpha = (0.000457*exp((13v)/50))
hBeta = (0.0065/(exp((v15)/28)+1))
hInf = hAlpha/(hAlpha + hBeta)
hTau = 1/(hAlpha + hBeta)
UNITSON
}
