I finished up the code that handles the number of neurotransmitter vesicles that fuse with the plasma membrane (and release into the synapse) in relation to calcium concentration in terminals. The code will autoselect the appropriate method of computation depending on the probability and total number of vesicles. For a large number of vesicles and low probability, the Poisson model is used. For a large number of vesicles and high probability, the binomial estimation method is used. For a low number of vesicles, the full binomial method is used.

Love the Poisson graph, as it has enough vesicles to get a good curve, but low enough probability to show the stochastic nature of the process.