The study of synchronous or asynchronous in (stochastic) neuronal populations is an important concept both in theory and in practice in neuroscience. The NPDS toolbox provides an interactive simulation platform for exploring such processes in Matlab looking through the lens of nonlinear dynamical systems. NPDS includes two main components: neural population (de)synchronization, and neural dynamics. One can investigate distribution controls on various neural models such as HH, FHN, RH, and Thalamic. Also, it supports many numerical approaches for simulation: finite-difference, pseudo-spectral, radial basis function, and Fourier methods. In addition, this toolbox can be used for population phase shifting and clustering.