%% general settings for the model neuron:
trials = 10;
tmax = 50.0;

%% f-I curves:
figure()
Ds = [0, 0.001, 0.01, 0.1];
for j = 1:length(Ds)
	D = Ds(j);
	inputs = 0.0:0.5:20.0;
	rates = ficurve(trials, inputs, tmax, D);
	plot(inputs, rates);
	hold on;
end
hold off;

%% spike raster and CVs
input = 12.0;
for j = 1:length(Ds)
	D = Ds(j);
	spikes = lifspikes(trials, input, tmax, D);
	subplot(4, 2, 2*j-1);
	spikeraster(spikes, 0.0, 1.0);
	subplot(4, 2, 2*j);
	isih(spikes, [0:0.001:0.04]);
end