function dvdt = fdvdt103(nv, vv, iext);% FDVDT103 - computes value of dV/dt for reduced, 2-variable HH model:% dV/dt = {iext - Gna*minf(vv)^3*nvtohv(nv)*(V-Ena) - Gk*ninf^4*(V-Ek) - Gleak*(V-Eleak)}/C% NOTE sethhparms must be called first to set up HH parameters.% Parameters (mS/cm^2, mV, uFd/cm^2)% hps = [gnamx=120; Ena=55; gkmx=36; Ek=-72; gleak=0.3; El=-49.4; C=1]hps = gethhparms;dvdt = (iext - hps(1)*minfhh(vv).^3.*nvtohv(nv).*(vv-hps(2)) - ...          hps(3)*nv.^4.*(vv-hps(4)) - hps(5)*(vv-hps(6)))/hps(7);return