Write Fortran k-order simulation routine

This is required as the Dynare++ routines have a heavy overhead. It should be used to replace local_state_space_iteration (#1773 (closed)) and will be used for simulated moments in the context of higher-order welfare (!1866 (closed)).