function [Ifish, Imut] = fig2_AB(N, fTau)
tic
stderr = 2e-2;
maxiter = 1e5;
tau = 1.0;
F = fTau .* tau;
alpha = 0.5;
fmax = 50.0;
fbg = 10.0;
sigma = 30.0;
nrns = [-180 : 360/N : 180-360/N];
stim = StimulusEnsemble('circular', 360, 360);
popNrns = CircGaussNeurons(nrns, sigma, fmax, fbg, tau, 'Gaussian-independent', [F alpha]);
Ifish = popNrns.Ifisher(stim);
[Imut, ImutSEM, ImutSamps] = popNrns.mi('randMC', stim, stderr, maxiter);
fprintf('fig2_AB.m\n')
fprintf('Parameters: N = %d neurons, F/tau = %g spikes/s^2\n', N, fTau)
fprintf('I_Fisher = %g bits\n', Ifish)
fprintf('I_mut = %g bits with StdErr %g bits\n', Imut, ImutSEM)
fprintf('I_mut - I_Fisher = %g bits\n', Imut - Ifish)