function dh=sodium_inactivation_hh(v,h)

alpha_h=.07*exp(-.05*(v+65));
beta_h=1/(1+exp(-.1*(v+35)));

dh=alpha_j*1(1-h) - beta_h*h;

[v,h] = ode15s(@sodium_inactivation_hh,[0 3000],[2 0]);