Implement pruning at order>3

The algorithm should follow Andreasen et al. (2018) as already implemented in simult_.m for orders 2 and 3. Should be done in the C++ routines of dynare_simul_ as discussed in 1e92e308 (comment 8364)