

%Noble equations

m = W(1,:);
n = W(3,:);
h = W(2,:);
      am = -0.1*(V+48.)/(exp(-(V+48.)/15.)-1.);
      bm = 0.12*(V+8.)/(exp((V+8.)/5.)-1.);
      an = -0.0001*(V+50.)/(exp(-(V+50.)/10.)-1.);
      bn = 0.002* exp(-(V+90.)/80.);
      ah = 0.17*exp(-(V+90.)/20.);
      bh = 1./(1.+exp(-(V+42.)/10.));

      gna = gnabar*m.^3*h ;
      gk1 = 1.2*exp(-(V+90.)/50.) + 0.015*exp((V+90.)/60.);
      gk2 = 1.2*n.^4;
      cna = gna*(V-Vna);
      ck = (gk1+gk2)*(V-Vk);
      can = gan*(V-Van);
      Fv = -(cna+ck+can)/cm;
      Fw(1) = am*(1.-m)-bm*m;
      Fw(2) = ah*(1.-h)-bh*h;
      Fw(3) = an*(1.-n)-bn*n;
	  
 
