MATLAB R2016b software was employed for simulations and the models were constructed in SIMULINK environment as interlaced subsystems using the provided building blocks. The simulations were run with MATLAB's ode45 solver with variable step size and relative tolerance of 10-3. Both model files are run through a single main script which lets the user select the stimulation modes for examining the effects of changing different