integrator: stacking results into numpy arrays
test_2d_uinverse: checking that the total population does not change