SciComp@UNC Courseware: MATH564

Interactive demo: Hodgkin-Huxley experimental m,n,h slopes

Initialization code

(*Sodium channel (m)*)
fm[V_,Vr_]=0.1 (25-(V-Vr));
am[V_,Vr_]=fm[V,Vr]/(Exp[fm[V,Vr]]-1);
bm[V_,Vr_]=4. Exp[-(V-Vr)/18.];
(*Potassium channel (n)*)
fn[V_,Vr_]=0.01 (10-(V-Vr));
an[V_,Vr_]=fn[V,Vr]/(Exp[fn[V,Vr]]-1);
bn[V_,Vr_]=0.125 Exp[-(V-Vr)/80.];
(*Chloride channel (h)*)
ah[V_,Vr_]=0.07 Exp[-0.05 (V-Vr)];
bh[V_,Vr_]=1./(Exp[0.1 (30-(V-Vr))]+1.);

αβPlot[Vr_]:=Plot[{am[V,Vr]/am[Vr,Vr],bm[V,Vr]/bm[Vr,Vr],an[V,Vr]/an[Vr,vr],bn[V,Vr]/bn[Vr,Vr],ah[V,Vr]/ah[Vr,Vr],bh[V,Vr]/bh[Vr,Vr]} ,
{V,-75,-60},PlotLegends->{"(Subscript[α, m])","(Subscript[β, m])","(Subscript[α, n])","(Subscript[β, n])","(Subscript[α, h])","(Subscript[β, h])"},Frame->True,PlotRange->All,GridLines->Automatic,FrameLabel->{"V [mV]","α/(Subscript[α, r]),β/(Subscript[β, r])"}]

Experimental m,n,h slopes relative to values at resting potential

(* Manipulate is set up for Plot. Modify as needed *)

Manipulate[ αβPlot[Vr],{{Vr,-68},-75,-60},
SaveDefinitions->True]