function dn=potassium_gate(n,v)

alpha_n=.01*(v+55)/(1-exp(-.1*(v+55)));
beta_n=.125*exp(-.0125*(v+65));

dn=0;
dn(1)=alpha_n*(1-n) - beta_n*n;