It looks simulateDSM operates correctly for nlev=2 and nlev=3 but not for nlev > 3:
for even nlev it gives same results as for nlev=2 and for odd nlev it gives same results as for nlev=3.
Matlab delsig toolbox shows the same same behaviour.
That is because simulateDSM was tested with the same input sin wave amplitude=1 while
the DSM model input range is increased with the quantizer nlev.
It looks simulateDSM operates correctly for nlev=2 and nlev=3 but not for nlev > 3: for even nlev it gives same results as for nlev=2 and for odd nlev it gives same results as for nlev=3.