% This file is licensed under the Allen Institute Terms of Use, which can be found at:
%http://www.alleninstitute.org/Media/policies/terms_of_use_content.html
%------------------------------------------------------------------------------------

%%%Script to generate plots for Supplement Fig. S11 (matched delta and bimodal
%synaptic distributions differing in 0.1% of their synapses)

clear variables;

%Plot settings
fnt = 36;
fntlab = 36;
xtc = [0;50;100]*10^-3;  %x-axis ticks for time in (ms)
ytc = [0;5;10;15;20]*10^-3; %y-axis ticks for membrane potential in (mV)
xtl = [0;50;100]; %x-axis labels for time in (ms)
ytl = [0;5;10;15;20]; %y-axis labels for membrane potential in (mV)

set(0,'DefaultFigureWindowStyle','docked');
set(0,'defaultAxesFontName', 'Arial');

scale=10^3;   % to convert from V^-1 to mV^-1
scale1 = 10^0; 

load SuppFigS11_frates.mat;

%%%%%FIRING RATE PLOTS
figure;
plot(t_vec_del,fout_del,t_vec_bidel,fout_bidel,'linewidth',4);
legend('Delta','Bimodal','Location','Southeast');
set(gca,'FontSize',fnt);
set(gca,'XTick',xtc);
set(gca,'XTickLabel',xtl);
set(gca,'FontSize',fnt);
ylabel('Firing rate (Hz)');xlabel('time (ms)');
axis square;

print2eps Fig_S11.eps
