#voltclamp.ode
param V=-60,Vo=-25,So=5,A=10,g=2,Vrev=-65
#initial condition
nO(0)=0.0
#activation functions
ninf=1/(exp(-(V-Vo)/So)+1)
tau=A/cosh((V-Vo)/(2*So))
#equation for activation
dnO/dt=-(nO-ninf)/tau
#auxiliary variables
#current
aux i=g*nO*(V-Vrev)
#maximum current
aux imax=g*ninf*(V-Vrev)
#logarithm of relative current
aux lndi=ln((g*ninf*(V-Vrev)-g*nO*(V-Vrev))/(g*(V-Vrev)))
done