function glu_in_cleft_over_time=single_syn_stsd_for_montecarlo(seed_values,times_over_which_to_solve,number_of_iterations)
%this mfile should be called from synaptic_dynamics_model_optimizer
%the destiny of this m-file is making predictions about [Ca]bouton, exocytosis of transmitter, [glutamate]cleft, and AMPAR conductance(s), among other things

% %finding maxima
% maxima_ind=[];maxima_val=[];
% [maxima_val(1) maxima_ind(1)]=max(glu_in_cleft_over_time_hx(AP_times(1):AP_times(1)+45));
% [maxima_val(2) maxima_ind(2)]=max(glu_in_cleft_over_time_hx(AP_times(2):AP_times(2)+45));
% maxima_ind(1)=maxima_ind(1)+AP_times(1);
% maxima_ind(2)=maxima_ind(2)+AP_times(2);

%seed_values=[membrane_potential_due_to_AP, cytosolic_calcium_conc, number_ca_channels, ca_channel_conductance, ca_channel_close_tau, cytosolic_calcium_decay_tau, RRP_size, ...
%        RP_to_RRP_tau, RP_size,PEP_size, PEP_to_RP_tau, PEP_to_RRP_tau, glu_per_vesicle, glutamate_cleft_tau, cleft_glutamate_conc];

options=odeset('RelTol',1e-4,'AbsTol',repmat([1e-4],1,length(seed_values)));

glu_in_cleft_over_time=[];
for i=1:number_of_iterations
    seeeds=seed_values;
    %the randomization here represents the small differences in the state of the synapse at the beginning of each sweep
    
    seeeds(1)=seeeds(1)*(rand*.4+.8);
    seeeds(2)=seeeds(2)*(rand*.4+.8);
% 	  seeeds(3)=round(seeeds(3)*(rand*.4+.8));
    seeeds(4)=seeeds(4)*(rand*2+.9);
%     seeeds(5)=
%     seeeds(6)=
    seeeds(7)=round(seeeds(7)*rand*.8+.6);
%     seeeds(8)=
%     seeeds(9)=
%     seeeds(10)=
%     seeeds(11)=
%     seeeds(12)=
    seeeds(13)=seeeds(13)*(rand*.2+.9);
%     seeeds(14)=
%     seeeds(15)=
    
    state_of_synapse=seeeds;
    cleft_glutamate_conc_record=[0 0];
    for j=1:length(times_over_which_to_solve)
        [T,Y] = ode45(@single_syn_coupled_equations,[0 times_over_which_to_solve(j)],[state_of_synapse],options);
        plot(T,Y)
        break,break
        state_of_synapse=Y(end,:);
        %calcium_conc_addon=rand*3; %due to AP
        %state_of_synapse(1,2)=state_of_synapse(1,2)+calcium_conc_addon;
        cleft_glutamate_conc_record=[cleft_glutamate_conc_record; [T Y(:,xxxx)]];
        state_of_synapse(1)=seed_values(1)*(rand*.4+.8);    %this represents the action potential
        %also need to round the discrete values (number of vesicles in each pool, number of calcium channels open...) before starting the next round -- really should round each timestep
    end
   glu_in_cleft_over_time=[glu_in_cleft_over_time cleft_glutamate_conc_record]; 
end