simult_.m: simulate user-requested order with k_order_solver

Also prevent inconsistent simulations; every two orders, the second derivative wrt the perturbation parameter will change the g_* matrices