function dm=sodium_activation_hh(v,m)

alpha_m=.1*(v+40)/(1-exp(-.1*(v+40)));
beta_m=4*exp(-.0556(v+65));

dm=alpha_m*(1-m) - beta_m*m;