function plot_simulated_waveform_withstagger(stagger)
%called from command line, plots simulated waveforms for illustrative purposes (including making illustrations)
%writing this rather diminutive routine raised the issue of whether that portion of the PSC due to glutamate release from a second vesicle would also begin with low slope
%if the second vesicle's contents are released onto receptors untouched by the first's, then the current passing through those receptors should also begin to increase slowly as the first few receptors are bound
%if the second vesicle's contents, however, are released directly onto the the same receptor cluster as exposed to the cloud of glu released by the first receptor, would this second cloud be expected to have any slow portion?

n=10;
linewidth=2

t=1:100;
A=5;
tau_rise=3;
tau_decay=15;

pure_psc=A*(1-exp(-t/tau_rise)).^n.*exp(-t/tau_decay);
pure_psc=-1*[0 pure_psc];

alpha=[zeros(1,stagger) pure_psc]; if ~stagger,alpha=alpha-.05*min(alpha);end
beta=[pure_psc zeros(1,stagger)];

figure(3),clf,hold on
plot(alpha,'r','linewidth',linewidth)
plot(beta,'b','linewidth',linewidth)
plot(alpha+beta,'k','linewidth',linewidth)
axis off
title(strcat('Stagger:',num2str(stagger)))