function old_fashioned
%y(1)=rate (probability) of release -- should release affect the probability of release? (uniquantal rule)
%y(2)=readily releasible pool size
%y(3)=reserve pool size
%y(4)=rate of transfer of vesicles from RP to RRP
%y(5)=size of post-exocytosis pool (increments each time there is a release, until retrieved) -- should there be kiss-and-run (exocytosis w/o vesicle fusion)
%y(6)=rate of transfer of vesicles from post-exocytosis pool to RP
%y(7)=rate of transfer of vesicles from post-exocytosis pool to RRP
%y(8)=amount of transmitter in the cleft
% y0(1)=0.5; %a probability
% y0(2)=10; %in vesicles
% y0(3)=)100; %in vesicles
% y0(4)=.01; %in vesicles per time step
% y0(5)=0; %in vesicles
% y0(6)=.01; %in vesicles per time step
% y0(7)=0; %in vesicles per time step
% y0(8)=0; %in pA(?)
calcium_conc=0;
RRP_size=10;
AP_times=[100 200]
aa=[]
for i=1:500
if find(i==AP_times)
calcium_conc=1;
else calcium_conc=calcium_conc*.9;
end
aa(i)=calcium_conc;
end
plot(aa)