% cell, stimulus and grid parameters

Ra = .3;	% k Ohm cm
Cm = 1;		% micro F / cm^2
Gm = 0.3;	% mS / cm^2
Gs = 2*Gm;	% mS / cm^2

%stim = inline('t*exp(-t)/2000');        % soma current stimulus
stim = inline('-(1+sign(t-5)*sign(40-t))/10000');  % soma current stimulus

dt = .1;	% time step, ms
Nt = 600;	% number of time steps
h = 1;		% space step, microns