Via the CTMC.simulate function.
CTMC.simulate.piecewise(rate.fns, jumpLim.fns, T.times, init.state)a LIST of rate functions corresponding to jumpLim.fns and T.times. Length is number of homogeneous pieces, we'll call it M.
a LIST of 'jumpLim' functions of length M like the list rate.fns. See the documentation for CTMC.simulate for an explanation of what each is.
Of length M+1 so that there are M intervals corresponding to rate.fns.
A starting state for the simulated chain.
An object of type CTMC.