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