working van der Pol oscillator code, everything up until now with 1D and 2D population units