function m_inf=mcalc(V)
      z= (0.1 *(25-V)./(exp(-0.1*(V-25))-1))./( 4*exp(-0.0556*V));
      m_inf =z./(1+z);