function [fm,fs,fmax,fmin] = RD_test_init(num_trials,num_cells)
%

f = zeros(num_trials,1);

for t = 1:num_trials
   [ta,tb] = RD_linear(num_cells,0.0001,1200,4,4,...
      0.25,0.0625,0.05,1);
   f(t) = RD_freq(tb);
end
fm = mean(f);
fs = std(f);
fmax = max(f);
fmin = min(f);
