Write two MATLAB codes,

a$)$ to calculate and plot $K^{+}$channel activation rate constants ${}_n(t)$ and ${}_n(t),$ and the activation

parameter $n(t),$

b$)$ to calculate and plot $N{a}^{+}$channel activation and deactivation rate constants ${}_m(t),{}_m(t),{}_h(t)$

and ${}_h(t),$ and the activation and deactivation parameters $m(t)$ and $h(t).$

Run the two MATLAB codes of parts $($a$)$ and $($b$)$ for the following three membrane voltages, $V_m(t),$

and comment on the results.

a$)V_m(t)=-90mV$ for $0t<mn>1</mn><mn>0</mn><mi>m</mi><mi>s</mi>$ and $V_m(t)=0$ for $10mst.$

b$)V_m(t)=-90mV$ for for $10mst<mn>2</mn><mn>0</mn><mn>0</mn><mi>m</mi><mi>s</mi>$ and $V_m(t)=-90mV$ for

$200mst.$

c$)V_m(t)=-90mV$ for for $10mst<mn>2</mn><mn>0</mn><mn>0</mn><mi>m</mi><mi>s</mi>$ and $V_m(t)=-90mV$ for

$200mst.$

You need to submit your MATLAB codes along with your plots and comments.

NOTE: You can assume that the resting membrane potential equals $V_r=-90mV.$