
h0 = rn0 - n0;
Ina = gnabar*aminf.*h0.*(Vn-Vna) ;
Ik = gkbar*n0.^4.*(Vn-Vk);
Icl = gl*(Vn-Vl);
      
Fvn = -Ina-Ik-Icl;