RNG in DiscreteProblem solver

[sdwfrost]

It’s basically this model (see code at end)

Tips on (sort of) creating a new Problem type in SciML, derived from an existing type

Is this code cleaner? Is composition the best way here, rather than use traits e.g. via SimpleTraits.jl? I mostly need the types for the time and the state, so the typing should be sufficient. What are the scenarios for using the isinplace type? # TODO: # Add MarkovSystem # Generalize to integer steps # Write integrator struct # Write init and step methods import SciMLBase: DiscreteProblem, ODESolution import OrdinaryDiffEq: solve, FunctionMap import Random: AbstractRNG import Plots: plot str…