function y=LR(t,x)

v=x(1);  % potencial de membrana
n=x(2);  % Subunidad n
Ca=x(3); % calcio

% Parametros
gKCa=0.02;
gK=3;
gCa=3.2;
Kd=1;
VK=-75;
VCa=100;
VL=-40;
Cm=.5;
gL=0.012;
va=30;
vp=50;
f=0.007;
kc=0.02;
k1=0.0275;

An=0.01*((10-(v+va))/(exp((10-(v+va))*0.1)-1));
Bn=0.125*exp(-(v+va)/80);
Am=0.1*((25-(v+vp))/(exp((25-(v+vp))*0.1)-1));
Bm=4*exp(-(v+vp)/18);
Ah=0.07*exp(-(v+vp)/20);
Bh=1/(1+exp((30-(v+vp))*0.1));

m_inf=Am/(Am+Bm);
h_inf=Ah/(Ah+Bh);

%ecuaiones
y(1)=-( gCa*m_inf^3*h_inf*(v-VCa) + (gK*n^4+gKCa*Ca/(Ca+Kd))*(v-VK) + gL*(v-VL)) / Cm;  
y(2)=An*(1-n)-Bn*n; 
y(3)=f*(-k1*gCa*m_inf^3*h_inf*(v-VCa)-kc*Ca);

y = y'*1000;