function r = rhsV(n) 

global vna vk vl gna gk gl c I v 

r = (I - gna*(m_inf(v))^3*(0.8-n)*(v-vna) - gk*n^4*(v-vk)-gl*(v-vl))/c; 
