function aH = h_inf(V)
Vi = -63;     % Sodium inactivation reversal potential (mV)
ki = 6;       % Inactivation slope
aH = 1./(1+exp((V-Vi)/ki));
end