global vna vk vl gna gk gl c I 

vna=50; vk=-77; 
vl=-54.4; 
gna=120; 
gk=36; 
gl=.3; 
c=1; 
I=20; 

[T,Y] = ode23('func_HHapprox',[0,100],[-65,.317]); 

figure(1); subplot(2,1,1) 
plot(T,Y(:,1)); 
xlabel('t'); 
ylabel('V'); 

subplot(2,1,2); 
plot(T,Y(:,2)); 
xlabel('t'); 
ylabel('n'); 
