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)=;