function dy=single_synapse_stsd_solve(t,y,i_vgcc)
%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
%y(5)=
%y(5)=
%y(5)=
dy=zeros(8,1);
dy(1)=;
dy(2)=;
dy(3)=;
dy(4)=;
dy(5)=;
dy(6)=;
dy(7)=;
dy(8)=;