In order to examine nerve impulses by means of simulation methodology, the models of all parts of nervous system, well suited for numerical modeling, are needed. In this paper the problem of setting up such a model, namely, that of a presynaptic bouton, is addressed. Simulation of the neurotransmitter flow inside the presynaptic bouton is performed. The transport is modeled with a partial differential equation with an additional nonlinear term. Two ways of modeling the bouton are applied. One of them let reflect a complex shape of the bouton and of some inner organelles. The influence of the generated mesh quality on the accuracy of numerical simulations is studied by comparing the released amount of neurotransmitter. The only mesh that produced diminished output was the worst one. The conclusion is that even slightly optimized tetrahedral mesh is suitable for calculations.