The excitatory synaptic function is subject to a huge amount of researches and fairly all the structural elements of the synapse are investigated to determine their specific contribution to the response. A model of an excitatory (hippocampal) synapse, based on time discretized Langevin equations (time-step = 40 fs), was introduced to describe the Brownian motion of Glutamate molecules (GLUTs) within the synaptic cleft and their binding to postsynaptic receptors. The binding has been computed by the introduction of a binding probability related to the hits of GLUTs on receptor binding sites. This model has been utilized in computer simulations aimed to describe the random dispersion of the synaptic response, evaluated from the dispersion of the peak amplitude of the excitatory post-synaptic current. The results of the simulation, presented here, have been used to find a reliable numerical quantity for the unknown value of the binding probability. Moreover, the same results have shown that the coefficient of variation decreases when the number of postsynaptic receptors increases, all the other parameters of the process being unchanged. Due to its possible relationships with the learning and memory, this last finding seems to furnish an important clue for understanding the basic mechanisms of the brain activity. © 2014 Springer Science+Business Media.