function rates = ficurve(trials, inputs, tmax, D)
% compute f-I curve.
rates = zeros(length(inputs), 1);
for k=1:length(inputs)
	spikes = lifspikes(trials, inputs(k), tmax, D);
	rates(k) = firingrate(spikes, 0.0, tmax);
end
end