%the right hand side for ode simulation:
function s_prime=deRHS(t,s)
global gam alf eps Inp BCL
V = s(1);
W = s(2);

Input = Inp/cosh(mod(t,BCL)-BCL/4);

Fw = eps*(V - gam*W);
Fv = -V.*(V-1).*(V-alf) - W ;
s_prime = [Fv+Input Fw]';
