%==========================================================================
%Modelo Hodgkin-Huxley
%EDOs
%==========================================================================
function dS=Hodgkin_Huxley_EDOs(t,S)
  
%Variables globales--------------------------------------------------------
global E_N
global E_K
global g_N_bar
global g_K_bar
global g_leak
global E_leak
global C_M

%Obtención de valores------------------------------------------------------ 
%Potencial de Membrana
V=S(1);
%Activadores
m=S(2);
h=S(3);
n=S(4);
 
%Resolución de dm/dt-------------------------------------------------------
m_inf = (0.10*(V+35)/(1.0 - exp(-(V+35)/10.0)))/((0.10*(V+35)/(1.0 - exp(-(V+35)/10.0)))+(4.0*exp( -(V+60)/18.0)));
t_m   = 1.0/((0.10*(V+35)/(1.0 - exp(-(V+35)/10.0)))+(4.0*exp( -(V+60)/18.0)));
dmdt  = (m_inf - m)/t_m;
%Resolución de dh/dt-------------------------------------------------------

%Resolución de dn/dt-------------------------------------------------------


%Cálculo de las corrientes iónicas-----------------------------------------
%Colocar la expresión de las corrientes conforme las conductancias y los
%valores de m,n y h



%Corriente de excitación---------------------------------------------------
%Debe ser un pulso de 1ms de valor 6.85mV en un tiempo a elección dentro 
%del intervalo 



%Derivada del potencial de Acción------------------------------------------   
%Expresión de dV/dt



%Salida   
dS=[dVdt dmdt dhdt dndt]';