######### PARAMETERs

# PARAMETERS - Applied Current
p IApp = 0

# CONDUCTANCEs
p gKDR  = 100
p gK    = 32
p gBK   = 80
p gCaPQ = 5
p gCaL  = 23
p gNaP  = 31
p gL    = 18.1



######### ODEs

# ODE - VOLTAGE
dV/dt = -(IK(mK,V)+IBK(mBK,hBK,V)+IKDR(mKDR,V)+ICaPQ(V)+ICaL(hCaL,V)+INaP(V)+IL(V)-IApp)/C

# ODE - ACTIVATION VARIABLEs
dhCaL/dt  = ( hCaLI(V)  -  hCaL ) / hCaLT
dmK/dt    = ( mKI(V)    -    mK ) / mKT
dmKDR/dt  = ( mKDRI(V)  -  mKDR ) / mKDRT
dmBK/dt   = ( mBKI(V)   -   mBK ) / mBKT
dhBK/dt   = ( hBKI(V)   -   hBK ) / hBKT



######### CURRENTS

IK(m,V)     = gK    *    m^4     * 1 * ( V  -   EK )
IBK(m,h,V)  = gBK   *     m      * h * ( V  -   EK )
IKDR(m,V)   = gKDR  *     m      * 1 * ( V  -   EK )
ICaPQ(V)    = gCaPQ *  mCaPQI(V) * 1 * ( V  -  ECa )
ICaL(h,V)   = gCaL  *  mCaLI(V)  * h * ( V  -  ECa )
INaP(V)     = gNaP  * mNaPI(V)^3 * 1 * ( V  -  ENa )
IL(V)       = gL    *     1      * 1 * ( V  -   EL )



######### ION CHANNEL DYNAMICS

# ICaL ------------------------------------------------
num vmCaL = -38, smCaL = 6
num vhCaL = -42, shCaL = -4
num hCaLT =  20

mCaLI(V)  = 1/(1 + exp(-(V-vmCaL)/smCaL))
hCaLI(V)  = 1/(1 + exp(-(V-vhCaL)/shCaL))


# ICaPQ ------------------------------------------------
num vmCaPQ = -22, smCaPQ = 8

mCaPQI(V) = 1/(1 + exp(-(V-vmCaPQ)/smCaPQ))


# IK -------------------------------------------------
num vmK = -55, smK = 15
num mKT = 2

mKI(V) = 1/(1 + exp(-(V - vmK)/smK))


# IKDR ------------------------------------------------
num vmKDR = -11, smKDR = 8
num mKDRT = 1.3

mKDRI(V) = 1/(1 + exp(-(V-vmKDR)/smKDR))


# IBK -------------------------------------------------
num mBKT = 2
num dhBK = 0.0025, KCa = 0.67, ghBK = 0.002
par hBKT = 5

mBKI(V) = mCaPQI(V)
hBKI(V) = ghBK/(mCaPQI(V)*dhBK*KCa*abs(V-ECa) + ghBK)


# INaP ------------------------------------------------
num vmNaP = -50, smNaP = 15

mNaPI(V) = 1/(1 + exp(-(V - vmNaP)/smNaP))



######### ELECTROPHYSIOLOGY

# CAPACITANCE
p C = 10

# NERNST POTENTIALs
p ECa  =  60
p EK   = -85
p EL   = -80
p ENa  =  50



######### SETTNGS

# VISUALIZATION
@ MULTIWIN=1
@ SIMPLOT=0
@ NPLOT=1
@ XP=t
@ YP=V
@ XLO=2500
@ XHI=3000
@ YLO=-90
@ YHI=10

# AUTO
@ NTST=30      
@ NPR=10000
@ NMAX=10000
@ DS=1e-3
@ DSMAX=1
@ DSMIN=1e-3
@ NORMMIN=0
@ NORMMAX=Inf
@ PARMIN=-300
@ PARMAX=1000
@ AUTOVAR=V
@ AUTOXMIN=-300
@ AUTOYMIN=-90
@ AUTOXMAX=1000
@ AUTOYMAX=20

# INTEGRATION
@ BOUND=1e50
@ TOTAL=5000
@ TOL=1e-6
@ ATOLER=1e-6
@ METHOD=C
done