global vna vk vl gna gk gl c I v 

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

for i=1:130 
v = -80 + i + 0.01; 
vv(i) = v; 
nullv(i) = fzero('rhsV',0.5); 
nulln(i) = fzero('rhsn',0.5); 
end 

figure(2) 
hold on; 
plot(vv,nullv,'b'); 
plot(vv,nulln,'r'); 
axis([-80 60 0 1]); 
xlabel('V'); ylabel('n'); 
