# comment!
init v=-60

# declare other stuff to plot
aux ica=gca*minf(v)*(v-vca)

# declare the parameters
param vl=-60,vca=120
param i=0,gl=2,gca=4,c=20
param v1=-1.2,v2=18

# user defined functions 
minf(v) .5*(1+tanh((v-v1)/v2))

# right-hand sides 
v'=(i+gl*(vl-v)+gca*minf(v)*(vca-v))/c
done