function xprime = SB(t,x)
global R L C Iin eps
xprime=zeros(3,1);
V_C=x(1);I_L=x(2);I_g=x(3);
%xprime(1)=-(I_g+I_L+Iin(t))/C;
xprime(1)=-(I_g+I_L+Iin)/C; 
xprime(2)=(V_C-R*I_L)/L;
xprime(3)=KNa(V_C,I_g)/eps;