Using 0.01 ms for the time step, implement the Hodgkin-Huxley equations in MATLAB. Use the parameter values given and apply a current pulse of amplitude Im =30 A/cm2 for a duration of 0.5 ms.

The HH parameters and initial values for the squid giant axon are given as follows:

1) Plot the action potential generated in response to the current pulse.

2) Plot the gating variables n, m and h as a function of time during an action potential.

Vm equilibrium = -64 mv Cm= 1uF/cm2 g Na = 120.0 mS/cm2 Ena = 50 mV Sk = 36.0 mS/cm2 g = 0.3 mS/cm Ek = -77 mv Ei= -55 mV m, = 0.0529 no = 0.31679 h, = 0.5963 The rate constants an, n, Am, m, An, Bn are calculated from the following equations: -65-V. m an 0.01(-55-V.) -55 - V exp -1 10 Br = 0.125 exp (5) 80 m a 0.1(-40 Vm) -40-V exp -1 10 Bm = 4 exp -65-V, 18 (6) m m on = 0.07 exp -65 Vm 20 20") 1 Bn= (-35-VM exp 10 (7) +1 Vm equilibrium = -64 mv Cm= 1uF/cm2 g Na = 120.0 mS/cm2 Ena = 50 mV Sk = 36.0 mS/cm2 g = 0.3 mS/cm Ek = -77 mv Ei= -55 mV m, = 0.0529 no = 0.31679 h, = 0.5963 The rate constants an, n, Am, m, An, Bn are calculated from the following equations: -65-V. m an 0.01(-55-V.) -55 - V exp -1 10 Br = 0.125 exp (5) 80 m a 0.1(-40 Vm) -40-V exp -1 10 Bm = 4 exp -65-V, 18 (6) m m on = 0.07 exp -65 Vm 20 20") 1 Bn= (-35-VM exp 10 (7) +1