function I = Iin(t)
global i1 A T phi 
I = i1*A*sin(2*pi*t./T+phi); % i1 -- spike threshold, |A|>1 iff activation spike.